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PREFACE 


7  This  final  report  describes  work  completed  under  Research  Number 
F33600-80-C-^*23  during  the  period  July  1980  through  December  1980.  This 
effort  was  performed  for  the  United  States  Air  Force  Logistics  Command 
at  Wright-Patterson  Air  Force  Base,  Ohio. 

The  work  was  performed  by  the  School  of  Industrial  Engineering  and 
Management  at  Oklahoma  State  University.  Dr.  M.  Palmer  Terrell  and 
Dr.  Philip  M.  Wolfe  were  co-principle  researchers,  assisted  by  Ph.D. 
students,  Mr.  Umit  Yuceer  and  Mr.  Shawn  Yu. 

This  effort  was  undertaken  to  develop  a  new  model  and  approach  for 
solving  the  War  Readiness  Spares  Kit  (WRSK)  problem.  The  development 
and  evaluation  of  a  new  model  and  computer  program,  referred  to  as  the 

Greedy  Algorithm,  is  described  herein. 

< 

The  final  report  is  divided  into  six  chapters  and  four  appendices.  1 
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PART  I 


INTRODUCTION 

Problem  Statement 

The  management  of  spare  parts  Inventories  is  indispensable  to  the 
logistics  operations  of  the  Air  Force.  Effective  inventory  management 
systems  are  of  great  importance  to  the  Air  Force  since  a  large  amount 
of  money  is  spent  annually  on  purchases  of  spare  parts.  Two  systems 
'METRIC  and  MOD-METRIC*  have  been  implemented  to  assist  in  the  manage¬ 
ment  of  recoverable  spare  parts.  However,  some  of  the  logistics  problems 
need  specialized  and  sophisticated  solution  strategies  and  methods. 

A  special  type  of  inventory  is  the  war  readiness  spares  kit  (WRSK) . 

This  consists  of  a  group  of  spare  parts  that  are  deployed  with  a  flying 
unit.  A  WRSK  system  is  in  operation  and  assists  in  determining  the  spare 
parts  that  should  be  In  a  kit.  The  objective  of  this  system  is  to  determine 
the  contents  of  a  kit  which  will  result  in  a  high  level  of  performance  at  a 
minimum  cost. 

An  important  aspect  of  this  problem  is  the  definition  of  measures  of 
supply  performance  which  are  operationally  meaningful.  A  thorough  investi¬ 
gation  of  the  problem  and  the  underlining  assumptions  suggests  that  the  NORS 
(Not  Operationally  Ready-Supply)  and  SDO  (Stock  Due  Out)  measures  can  be  used 
to  describe  the  performance  of  the  kit.  (See  references  [3]  and  [6].)  The 
NORS  function  is  non-linear  and  nonseparable  over  the  nonnegative  integers. 

This  makes  the  problem  difficult  to  solve.  In  addition,  the  number  of  items 
in  the  kit  is  very  large  in  practice.  Consequently,  a  large  size  non-linear 
integer  programming  needs  to  be  solved.  Dynamic  programming  [1]  or  the  branch- 


and-bound  method  [5]  could  be  employed  to  solve  the  WRSK  problem,  but  the 
computation  time  required  is  impractical. 

Brief  Description  of  Current  D029  WRSK  Model 

The  current  D029  Model  and  its  computer  propram  is  designed  to  solve 
this  problem.  It  has  two  passes.  The  first  pass  tries  to  determine  a  kit 
which  will  yield  a  high  performance  for  the  initial  six  day  period  for  a 
squadron  of  aircrafts.  The  second  pass  attempts  to  find  the  contents  of 
the  whole  kit  for  the  thirty  day  period.  A  detailed  description  is  pro¬ 
vided  by  Wright  Patterson  Air  Force  Base  Logistics  Command  in  [8]. 

This  model  employs  marginal  analysis  to  obtain  a  kit  which  will  perform 
better  than  the  requirement  kit  57-18.  The  tradeoff  between  NORS  and  SDO  is 
handled  by  manipulating  the  weight  between  them  at  each  step  according  to 
some  heuristic  rules.  In  addition,  at  each  step  of  the  algorithm  some  number 
of  units  of  an  item  are  added  to  the  kit.  This  number  is  also  manipulated 
according  to  some  heuristics.  This  is  necessary  because  the  current  D029  Model 
uses  marginal  improvements  which  must  be  recalculated  after  one  unit  of  an 
item  is  added  to  the  kit  due  to  the  nonseparability  of  the  expected  NORS  func¬ 
tion,  E(NORS/X) . 

The  current  D029  computer  program  has  some  special  subroutines  to  handle 
input,  word  packing  and  unpacking,  and  calculation  and  interpolation  of  pro¬ 
babilities.  One  of  the  subroutines  provides  intermediate  kits  at  various 
budget  levels. 


Problems  With  Current  Solution  Methodology 

The  marginal  analysis  as  currently  used  will  not  guarantee  an  optimum 
solution.  Fox  [4]  states  conditions  when  discrete  marginal  analysis  will 


2 


« 

T 

S 

i 


result  in  an  optimum  answer.  One  criterion  is  that  the  functions  must  be 
separable.  Hie  model  used  in  the  WRSK  system  involves  functions  which  are 
not  separable.  Consequently,  the  marginal  analyis  approach  as  applied  to 
the  WRSK  must  be  categorized  as  a  heuristic. 

Another  problem  involves  the  calculation  of  E(NORS/X)  and  E(SDO/X)  functions. 
These  expectation  formulas  are  truncated  to  ease  the  computations.  Truncation 
may  cause  numerical  errors  if  the  upper  bound  of  the  truncation  is  not  chosen 
suitably. 

Relative  to  the  current  WRSK  system,  models  involving  large  flying  units 
(greater  than  24  aircraft)  require  a  large  amount  of  memory.  So  much  so,  that 
portions  of  the  resident  software  must  be  removed  to  provide  the  required  stor¬ 
age.  In  order  to  avoid  the  use  of  excessive  storage,  word  packing  and  unpacking 
is  employed.  This  however  creates  some  numerical  errors  throughout  the  compu¬ 
tations.  Some  of  the  probability  values  are  stored  and  others  are  interpolated 
when  needed.  This  causes  some  probabilities  to  be  underestimated.  More  impor¬ 
tantly,  this  process  is  very  time  consuming.  Approximately  half  of  the  execution 
time  is  spent  on  packing  and  unpacking. 
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PART  II 


OBJECTIVES  OF  THIS  RESEARCH 

General 

The  objective  of  this  research  has  been  to  develop,  validate,  and 
document  an  improved  War  Readiness  Spares  Kit  (WRSK)  model  and  solution 
procedure  which  may  be  used  by  the  United  States  Air  Force  in  the 
management  of  recoverable  spare  parts.  * 

This  research  effort  has  been  achieved  by  accomplishing  a  series 
of  subobjectives.  These  objectives  are  stated  below,  and  reference  is 
made  to  the  pertinent  parts  of  this  report. 

Objectives 

Objective  1:  Develop  and  validate  an  improved  solution  procedure  for 
the  current  WRSK  model. 

Development  and  validation  of  such  a  solution  method  was  accomplished 
in  two  stages.  In  the  first  stage,  measures  of  supply  performance  and 
their  mathematical  properties  were  investigated.  Part  III  of  this 
report  is  devoted  to  the  results  of  this  investigation.  The  second  stage 
was  the  development  of  a  mathematical  model  which  is  reported  in  Part  IV 
where  a  multiobjective  approach  is  proposed  for  the  WRSK  problem.  Poten¬ 
tial  solution  methods  were  then  investigated.  A  small  example  problem 
was  solved  by  these  algorithms  in  order  to  demonstrate  how  they  generate 
solutions  and  how  they  converge  to  an  optimum.  Evaluation  of  these 
methods  led  to  the  selection  of  a  preferred  solution  method  which  was 
investigated  and  tested  using  F4D  data.  The  final  section  of  Part  IV 
discusses  the  current  D029  model  which  has  features  pertinent  to  this  research. 
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Objective  2:  Develop  procedures  which  will  eliminate  the  problems 
caused  by  using  expected  value  approximations. 

Part  III  considers  the  measures  of  supply  performance  which  are 
closely  related  to  the  stocastic  behavior  of  a  WRSK  kit.  Accurate 
expected  value  computations  are  of  vital  importance  to  the  generation 
of  high  quality  solutions.  Truncation  of  an  infinite  sum  is  compared 
to  the  approximation  of  the  sum.  Proper  use  of  the  truncation  procedure 
is  explained  by  some  examples  in  the  last  section  of  Part  III. 

Objective  3:  Implement  the  algorithm  to  solve  the  War  Readiness  Spares 
Kit  problem. 

A  major  aspect  of  this  research  has  been  to  implement  the  proposed 
solution  method  and  use  it  to  solve  a  WRSK  problem  with  real  data.  The 
implementation  and  computational  features  using  F4D  data  is  discussed  in 
Part  V. 

Objective  4:  Compare  the  performance  of  the  algorithm  with  the  current 
system  in  terms  of  time  and  computational  feasibility. 

In  order  to  establish  a  base  for  meaningful  comparison,  the  current 
D029  computer  program  was  modified  for  compatibility  with  IBM  compilers. 
Then,  performance  of  the  two  algorithms  was  observed  in  terms  of  time 
and  computational  feasibility.  These  results  are  discussed  in  Part  V. 

Objective  5:  Document  results  and  computer  programs. 

This  report,  in  its  totality,  documents  the  research  results.  The 
detailed  documentation  and  FORTRAN  listing  of  the  proposed  solution 
method  is  given  in  Appendix  B.l  and  B.2.  In  addition,  the  FORTRAN  listing 
of  the  IBM  modified  D029  program  is  provided  in  Appendix  C.l  and  c.2. 
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PART  III 


MEASURES  OF  SUPPLY  PERFORMANCE 

Introduction 

This  section  focuses  on  measures  of  supply  performance.  The  stochastic 
prediction  of  the  behavior  of  a  kit  during  the  mission  is  of  vital  importance 
in  the  study  of  the  WRSK  problem.  Another  important  aspect,  when  measuring 
the  performance  of  a  kit,  is  the  relationship  of  the  performance  measures  to 
the  operations. 

The  performance  of  a  kit  can  be  measured  in  several  ways;  fill  rate, 
operational  rate,  average  back  orders,  average  number  of  Not  Operationally 
Ready  Supply  Aircraft,  and  kit  cost.  These  can  be  described  briefly  as 
follows.  Fill  rate  is  the  ratio  of  the  number  of  parts  issued  over  a  fixed 
time  period  to  the  number  of  parts  demanded  over  the  same  period.  Operational 
rate  is  the  probability  that  there  will  be  no  due  outs  from  the  base  supply 
at  any  time.  Back  orders  contain  the  number  of  stock  due-outs  (SDO)  from  base 
supply  during  a  fixed  period.  Finally,  NORS  is  the  number  of  grounded  aircraft 
due  to  lack  of  spare  parts.  Experience  has  revealed  that  average  SDO,  average 
NORS,  and  the  cost  of  a  kit  relate  to  operations  better  than  the  other  measures. 

Mathematical  Development  of  Performance  Measures 

The  following  definitions  and  assumptions  are  made  in  order  to  develop 
the  mathematical  expressions  for  average  back  orders  and  NORS. 

(1)  X  represents  a  kit  consisting  of  individual  items  i  of  quantity  xj.. 

(2)  n  items  are  under  consideration  for  the  kit. 
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(3)  The  demand  of  an  item  is  a  Poisson  distributed  random  variable. 

(4)  The  failure  of  any  item  is  independent  of  the  failure  of  any  other 
item. 

(5)  When  a  demand  for  an  item  can  not  be  satisfied  from  the  stock  on 
hand,  the  required  part  is  consolidated  by  cannibilization. 

Let  F^(-)  denote  the  cumulative  Poisson  sums  for  the  item  i.  The  opera¬ 
tional  rate  given  no  more  than  k  planes  for  cannibilization  is  expressed  by 
the  statement 

bk  -  P[NORS  _<  k/X]  *  ir  Fi(xi  +  kaf) ,  k  =  0,1,2 . N 

i-1 

where  af  is  the  number  of  applications  of  the  part  i  on  an  airplane,  and  N  is 
the  total  number  of  aircraft  in  the  flying  unit.  Consequently 


Ab^ 

-  P[# 

of  NORS 

=  k] 

for  k  -  0,1,2,. 

,..,N 

then 

Ab0 

=  P[# 

of  NORS 

-  0] 

“  bo 
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Finally 

E(N0RS/X) 

N 

-  E 

(bk- 

bk-l) 

k-1 

An  equivalent  development  of  the  expectation  of  the  NORS  function  for  a 
given  kit  X  goes  as  follows. 

E(N0RS/X)  -  E  k  P[#  of  NORS  =  k] 
k-0 

E(N0RS/X)  *  I  (1  -  P[#  of  NORS  <  k]) 
k-0 

E (NORS /X)  -  E  (l  -  “  F±(xt  +  kai)) 
k-0  i-1 

An  approximation  to  the  above  function  for  computational  purposes  was 
proposed  byB.  L.  Miller  in  (6)  as  follows. 
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E(HORS/X)  «  E  (1  -  w  Fi(Xi  +  ka1))+  E  El-  F^Xi  +  kai) 
k-0  1=1  k=nri-l  1=1 

n 

where  m  Is  such  that  ir  Fi(xi  +  kai)  is  very  close  to  1  for  all  k>m.  See 

1=1 

Figure  3.1. 

The  expected  back  orders,  SDO,  for  a  given  kit  X  is  obtained  by 

E(SDO/X)  =  E  E  (l  -  Fi(X1  +  j)) 
i-1  j=o 

Reference  may  be  made  to  [3]  and  [8]  for  further  details. 

.  n 

If  item  i  costs  $ci,  the  cost  of  the  kit  is  then  E  ci  xi  or  in 

1=1 

vector  form,  CX. 

Mathematical  Properties  of  Performance  Measures 

A  few  properties  of  the  measures,  E(NORS/X)  and  E(SDO/X),  are  presented 
below.  They  will  be  used  later  in  the  development  of  solution  methods  and 
algorithms . 

Property  1 


Both  E(SDO/X)  and  E(NORS/X)  are  monotonically  decreasing  because  the 
first  differences  are  nonpositive. 


Define  Ak  E(SDO/X)  *  E(SDO/X  +  ek)  -  E(SDO/X)  where  ek  represents 
the  kth  unit  vector.  Then 

Ak  E(SDO/X)“  E  E  [1  +  FiCxi  +  j)]+  E  [l  +  Fk(xfc  +  1  +  j)]-  E  E  [l-Fi(xi+j)] 
i=l  j=0  j=0  i=l  j=0 

if»k 

00  oo 

=  E  (1  -  FfcCxk  +j)]  -  E  [1  -  Fi  (xfc  +j)]  which  yields 
j-1  J-o 


Ak  E(SD0/X)  •  -[1  -  Ffc  (x^)},  the  marginal  improvement  in  average  SDO 


when  a  one  unit  of  item  k  is  added  to  the  kit. 


Likewise,  the  marginal  improvement  in  average  NORS  is  calculated  as 
follows. 

Define  A*  E(NORS/X)  -  E(NORS/X  +  ek)  -  E(NORS/X).  Then 


Ak  E(NORS/X)  *  Ml  -  ?  Fi(xi  +  Jai)Fk(dk  +  jak  +1)]-  ?  [1  -  5  Fi(x±  +  jaj.)  ] 
j-o  i-1  j-o  j-1 

i-k 


■  S  (l  *  Fi(xi  +  3ai)l  tFk(xk  +  jakJ-FkCxk  +  ja^  +  1)]) 
j-o  i-1 
il*k 


*  -  I  Pk  (xk  Jak 
j-o 


“  Fi  (x±  +  ja±) 
i-1 
ij*k 


where  Pk(8) 


-^k  8 
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Property  2 

E(SDO/X)  >  E(NORS/X)  for  any  X 


Proof: 

E(SDO/D)  -  ?  S  [1  -  Fi(xi  +  j)]  +  Z  Z  [l-Fi(xi+j>] 

j-1  i-1  j-nrfl  i-1 

E(NORS/X)  -  Z  [1  -  x  Fi(xi  +  jai)  ]  +  Z  Z  [1  -  Fi(xi  +jai)] 
j-o  i-1  j-m+1  i-1 
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Define  DIF  -  F.<SDQ/X)  -  F.(NORS/X) 

-  ?  j?  [l-F1(i1  +  i)]-l  +  !  Ft  (x±  +  jai)l  +  E  E  [FiCxi  +  ja^-FiCxi  +j)J 
j«o  li-l  i-1  ’  j-nrfl  i-1 


-  E  in  -  1  -  E  FjCxi  +j)  +  J  FitXi  +  ja±)}  +  E  E  [F^  +jai)  -  F1(xi  +j)] 
j=*o  j-1  i-1  j— m+1  i— 1 


Since  a^  >_  1  for  all  i,  then  1  ■>  Fi(x^  +  jaj)  >_  F^(x^  +  j)  and 


n  n 

v  Fi(xi  +  Jai)  Z.  x  Fi(xi  +  j)t  further 
i=l  i-1 


$  Fi(xi  +  J)  -  5  Fi(xi  +  jai)  1  E  Fj(Xl  +  j)  -  5  FiCxi+j). 
i-1  i-1  i-1  i-1 


n  n 

Let  f(Y)  *  I  -  ir  y^.  The  maximum  of  this  function  is  (n-1).  The  proof 
i-1  i-1 


goes  as  follows: 

n  n  n-1 

Vf  -  (1-  x  ylf  . . 1  -  ir  yA, . . . ,  1  -  ir  y*) 
i-2  ii«k  i-1 


v2y 


o  a  j  2*  « 
32 1  O  •  • 


aln 

a2n 


ani  ®n2 


®nn 


where  j  =  -  tt  y± 
i?*k 

Vi 


and 


|akj | 1 


Clearly  rfc  V2Vr  <  0  reRn+.  Hence  f  is  concave  in  the  positive  unit 

hypercube.  Furthermore,  7f(y)  -  0  yields  the  solution  that  y±*  -  1  for 

all  i-1, . .  •  »n.  Consequently  T<1)  -  n-1  is  the  maximum  of  the  function. 

Finally  DIF  >  0,  since  (n-1)  -  [  E  Fi(xi  +  j)  “  /  Fi<xi  +  3ai>>  >  0 

i-1  i"l 

and 

E  [Fitxi  +  jai)  -  F1(x1  +  j)]  1  0  for  all  j  1  0, 
i-1 
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Property  3 

The  E(NOR S/X)  function  can  be  approximated  by  a  separable  function. 

This  is  accomplished  as  follows  (also  see  [3]). 

E(N0RS/X)2l  (nrt-1)  +  E  E  bv  log  F<  (x<  +  kat)  -  E  [1  -  FWxi  +  kat)  ] 

i-1  k-o  k-nrfl 

where  b^  *  n  +  ka^)  for  some  x  *  (Ifj  ,X2, . .  .If  ). 

i-1 

n 

where  n  is  such  that  it  f±(*i  +  ka^)  is  very  close  to  1  for  all  k>m. 

i-1 

m  n  m  n 

The  term  E  [1  -  it  FA(x^  +  ka*)  ]  can  be  rewritten  as  (nrt-1)  -  I  it  F^xj  +  kaj) 

k*o  i-1  k-o  i-1 

The  product  x  F^x^  +  kaj)  can  be  replaced  by  some  multiple  of  its  logarithm 
i-1 

in  the  following  manner. 

n  n 

-b  E  log  F^x*  +  kaj)  where  is  a  parameter  defined  by  bk  -  "  F^x*  +  ka*) 

i-1  i“l 

for  some  K  «  (Tfj ,  T£z  ,  .  • .  >  . 

Finally  then,  separability  is  accomplished  as  E(NORS/X)  -  (nrt-1)  + 

E  E  bk  log  Fi(xj_  +  ka^)  +  E  [1  -  Fi(xi  +  kaj)]. 
i=l  k-o  k-m+1 

This  results  plays  an  important  role  in  the  development  of  an  algorithm 
which  is  discussed  in  Section  IV. 
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Expected  Value  Calculations 


[  Calculation  of  the  expected  NORS  function  for  a  given  kit  X;  E(NORS/X)  = 

•  n 

X  (1  -  ir  F^(x^  +  ja^)  creates  computational  difficulties.  This  expectation 
J-l  i-1 

is  an  infinite  sum.  Since  it  is  convergent,  the  terms  become  negligibly  small 

^  for  sufficiently  large  j's.  There  are  two  ways  to  handle  the  computation. 

One  way  is  to  truncate  this  infinite  sum  for  a  sufficiently  large  value  M. 

The  second  way  is  to  develop  an  approximation  as  suggested  by  B.  L.  Miller  in 

^  [61.  Physically  the  number  MORS  airplanes  can  not  exceed  the  number  of  airplanes 

In  the  squadron.  This  then  establishes  a  natural  upper  bound  for  truncation 

Both  truncated  and  approximate  expected  MORS  formulas  are  provided  below,  and 

the  difference  between  them  is  evaluated. 

M  n 

a)  Truncation:  f~(X)  “X  (1  -  ir  Fi(x^  +  ja^)) 

j-0  i-1 

b)  Approximation:  f^(X)  •  X  (1  -  ir  F^Cx*  +jaj))  +  E  X  (1  ~  F^Cx^  +  jaj)) 

j-0  i-1  j-m+1  i-1 

where  M  is  the  number  of  planes  in  the  squadron,  a±  is  the  number  of  applications 

of  the  part  i  on  an  airplane,  and  m  is  such  that  9  Fi(xi  +  Jai)  is  very  close 

i-1 

to  1  for  all  j>m. 

The  difference  is  given  by  the  relation  A  -  f^(X)  -  f-j(X).  If  m  >  M,  then 

A  -  ?  (1  -  S  F^(xi  +  ja^)  +  X  X  (1  -  Fi(x^  +  ja^). 

j-M+1  i-1  j-nrfl  i-1 

If  m  _<  M,  then  A  -  X  X  (1  -  F^(x^  +  ja^)  -  X  (1  -  v  Ff  (x^  +  ja^)) 

j—m+1  i-1  j-m+1  i-1 

or  A  -  X  X  (1  -  Fi(x1  +  jaj)  -  X  X  (1  -  F^xi  +  ja*)) 
j-m+1  i-1  j-m+1  i-1 

because  the  approximation  is  valid  for  all  j>m.  Finally 
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A  -  E  E  (1  -  +  jai)) 

j-ttfl  i-1 

In  either  case,  A  is  a  positive  tern  and  A  may  not  be  negligibly 
small  especially  for  large  squadrons  with  some  items  having  heavy  demand 
rates . 

The  following  examples  illustrate  the  magnitude  of  the  difference 
between  the  truncated  and  approximated  formulas. 


Example  1 


Suppose  that  there  are  three  planes  in  a  squadron,  and  there  are  two 
essential  items  A,  B  with  demand  rates  y^  -  5  and  yfi  =■  5.5.  The  initial 
stock  levels  are  determined  to  be  x^  =  5  and  xB  ■  6.  For  simplicity, 
aA  -  aB  -  1. 

a)  Calculation  of  NORS  by  fT(X) 


_JL 

Fa  (5+1) 

Product 

1  -  Product 

0 

.616 

.686 

.42258 

.57742 

l 

.762 

.809 

.61646 

.38354 

2 

.867 

.894 

.77510 

. 22490 

fT(x)  =  1.18586 

b)  Calculation  of  NORS  by  f^(X) 
m  is  chosen  as  6 


JL 

lAiltll 

7jl(6+1) 

Product 

1  -  Prodt 

0 

.616 

.686 

.42258 

.57742 

l 

.762 

.809 

.61646 

.38354 

2 

.867 

.894 

.77510 

. 22490 

3 

.932 

.946 

. 88167 

.11833 

4 

.968 

.975 

. 94380 

.05620 

5 

.986 

.989 

.97515 

.02485 

6 

.995 

.996 

.99102 

.00898 

Sum  1.39422 


1  -  FB(6+j) 


It 


F 


7 

8 
9 

10 

11 


.00202 
.00070 
.00022 
.00007 
. 00002 


.00160 

.00060 

.00020 

.00006 

.00001 


Sum  .00303 


.00247  TSUM  =  .00550 


fA(X)  -  1.39972 


A  *  0.21385 


P 


Example  2 


Suppose  there  are  three  items  A,  B  and  C  with  pc  m  4.5  and  xc  -  5 
m  ■  6  again. 


a)  Calculation  by  truncated  formula 


FA<5+j) 

FB(6+J) 

FC(5+j) 

Product 

1  -  Product 

0 

.616 

.686 

.720 

.30426 

.69574 

1 

.762 

.809 

.884 

.52029 

.47971 

2 

.867 

.894 

.921 

.71387 

.28613 

fT(X)  -  1.46158 


b)  Calculation  by  Approximated  formula  fA(X) 


*A(5+j) 

F*<6+j) 

Fc(5+j) 

Product 

1  -  Product 

0 

.616 

.686 

.720 

.30426 

.69574 

1 

.762 

.809 

.844 

.52029 

.47971 

2 

.867 

.894 

.921 

.71387 

. 28613 

3 

.932 

.946 

.964 

.84993 

.15007 

4 

.968 

.975 

.985 

.92964 

.07036 

5 

.986 

.989 

.994 

.96930 

.03070 

6 

.995 

.996 

.998 

.98904 

.01096 

Sum  ®  1.72367 
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1  -  FA(5+j) 

1  ~  FB(6+j) 

1  "  Fc(5+j) 

7 

.00202 

.00169 

.00066 

8 

.00070 

.00060 

.00020 

9 

.00022 

.00020 

.00006 

10 

.00007 

.00006 

.00002 

11 

.00002 

.00001 

— 

.00094  TSUM  -  .00644 

fA(X)  -  1.73011 
A  -  .26853 


These  small  examples  show  that  there  is  an  error  which  Is  not  negligible. 
Even  If  the  number  of  airplanes  in  the  squadron  is  sufficiently  large,  one 
has  to  be  careful  to  reduce  the  magnitude  of  this  error.  In  addition,  there 
will  be  round-off  errors  for  large  size  problems  even  with  double-precision 
arithmetic.  It  is  essential  to  know  the  origin  of  the  errors  if  one  wishes 
to  reduce  the  magnitude  of  the  errors. 
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PART  IV 


< 


I 


MATHEMATICAL  MODELLING  AND  SOLUTION  METHODS  FOR  THE  WRSK  PROBLEM 
A  Multiobjective  Approach  to  the  WRSK  Problem 

In  simple  terms,  the  objective  of  the  WRSK  problem  is  to  determine 
the  quantity  of  each  part  in  the  kit  which  will  minimize  the  average  NORS, 
average  SDO,  and  the  cost  of  the  kit  simultaneously.  All  of  these  three 
objectives  conflict  with  each  other.  If  cost  is  minimized,  then  the  other 
two  measures  increase.  If  average  NORS  is  minimized  for  a  fixed  cost,  the 
average  SDO  turns  out  very  high.  Conversely,  the  average  NORS  comes  out 
high  when  the  average  SDO  is  minimized  for  a  fixed  budget  level.  The 
relationship  between  average  NORS,  SDO,  and  the  cost  was  discussed  in  [3] 
and  [8].  The  WRSK  problem  is  a  multiobjective  programming  problem,  and 
the  model  is  stated  mathematically  as 

Min  {E(NORS/X),  E(SDO/X),  CX} 
s.t.  X  £  0,  integer 

This  model  is  also  operationally  meaningful,  since  all  of  these 
measures  highly  relate  to  operations.  E(N0RS/X)  measures  the  expected 
number  of  cannibalized  airplanes.  E(SD0/X)  is  an  indicator  of  the  addi¬ 
tional  work  load  for  repairs,  and  CX  gives  the  cost  of  the  kit. 

The  following  solution  strategy  is  proposed  for  solving  this  multi¬ 
objective  programming  problem.  First  of  all,  the  problem  can  be  considered 
as  a  parametric  programming  in  terms  of  budget  levels  ’b'  as  shown  below. 

Min  (E(N0RS/X) ,  E(SD0/X)) 
s.t.  CX  £  b 

X  £  0,  integer. 
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This  reduced  problem  may  be  further  investigated  as  two  single¬ 
objective  problems  as  follows: 

(1)  Min  E(NORS/X) 

s.t.  CX  £  b 

X  _>  0,  integer 

and 

(2)  Min  E(SD0/X) 
s.t.  CX  £  b 

X  ^  0,  integer 

Obviously  the  optimal  solution  of  (1)  is  not  necessarily  optimal  to 
(2),  and  conversely.  If  X*  solves  (1)  then  E(SDO/X*)  turns  out  to  be 
very  high,  and  conversely  if  X*  solves  (2),  then  E(NORS/x*)  is  very  high. 

This  phenomenon  has  also  been  observed  by  WPAF  Logistic  Command  and  was 
stated  explicitly  in  [8],  the  Preliminary  Evaluation  of  D029  WRSK  Model. 

Since  E(SD0/X)  E(N0RS/X)  for  a  given  kit  X,  a  best  compromise  solution 
between  the  E(N0RS/X)  and  E(SD0/X)  is  sought  for  a  specified  budget  level. 
Figure  4.1  illustrates  the  relation  between  E(SD0/X)  and  E(N0RS/X)  for 
specified  budget  levels. 

Problems  (1)  and  (2)  can  be  merged  into  one  problem  (3)  by  the  use 
of  a  weight  w  as  follows. 

(3)  Min  (E(N0RS/X)  +  w  E(SD0/X)} 

CX  <  b 

X  ^  0,  integer 

The  weight  w  expresses  the  quantity  of  NORS  planes  equivalent  to  each 
SDO.  Prior  assessment  of  the  weight  is  required  to  obtain  the  best  compromise 
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Performance  Measures  Vs.  Various  Investment  Levels 
(Figure  taken  from  [3]) 
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Potential  Solution  Methods 

Many  potential  solution  methods  were  investigated.  Those  that  were 
found  significant  to  the  solutions  of  problem  (1),  (2),  and  (3)  are 
summarized  as  follows: 


A  Primal  Search  Algorithm 


.5 


A^.Ltt;  X)  -  L(X  +  e±;  X)  -  L(X;  X)  for  i  -  1,  ....  n  (3.7 


A±L(X;  X) 


-  Z  ((F.  (x  +  j  +  1) 
j-0  x  1 


F AX  +  j))  iX(\  +  j) 
1  1  k“l 
k/i 


Fi<*i  +  m  +  1)  +  XCA 


Since  0  ir  F  (x  +  j)  <_  1  for  0  <_  j  <  m,  then 
k*i  1  i 


:  V*!  +  J>  "  F1(k±  +  m  +  1)  +  Fi(xi)  +  1^0^  +  j)  -  1  + 


F^Cx^  +  m  +  1)  +  Xc^  <_  A^L(X;X  ) 


W  '  1  +  Xci  -  AiL<*i  *>  1  V*!  +  «  +  1)  -  1  +  XC± 


Therefore,  it  suffices  to  find  an  X°  *  (x°,  x“)  where  each 

in 


x°  satisfies  -1  +  F^  (x^  -  1)  +  XC^  <_  0  and 


- 1  +  F^(x^)  +  >0  i  *  lj  •••>  n* 


Step  1.  Search  for  a  better  feasible  solution;  set  k  =  k  +  1. 
Construct  a  favorable  direction 


dR  =  -Vf(X  )  -  p-pk-C 

where  7f(Xk)  =  (Ajf  (Xk) , . . .  ,Anf  (Xk)) ,  A^X11)  - 

-  I  f F. (xk  +  j  +  1)  -  F,(xk  +  j)l  H  F,(xk  +  j)  -  1  +  F.(xk  +  m  +  1) 
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and  pk  is  obtained  by  solving  (-VfCx1')  -  p^.C  -  0  and  p  is  left  as 
an  input  parameter  by  the  user.  dfc  is  normalized  such  the  max{jd  ^|; 


i  *  1,  . ..,  n)  *  1.  Calculate  an  a  by  a 


Th-c^I 
[Vn  C.dk  J 


where  [ • • J  denotes 


the  greatest  integer  function.  Then  a  point  ya  is  obtained  as  ya  =  Xk  + 
ad,  .  The  nearest  integer  point  V  to  y  is  obtained  by  v.  *  Ty  +  .5*1 

K  ait''  u{  J 

for  i  *  1,  ....  n. 

i)  If  V  is  feasible  and  f(V)  <  f (Xk) ,  set  X^+*  =  V,  to  to  Step  1. 
ii)  If  V  is  feasible,  but  f(V)  f (Xk) ,  then  check  the  points 

V  +  e£  for  i  =  1 ,  ...»  n.  Figure  4.3.  If  there  exists  an  i*  such 

k  j/i.1 

that  V  +  e.*  is  feasible  and  f(V  +  e..)  <  f(X  ),  set  X  *  V  + 

l*  l* 

e^*,  and  go  to  Step  1.  Otherwise,  go  to  Step  2. 


x2 


iii)  If  V  Is  not  feasible,  then  search  the  points  V  -  e^  for  i  -  1, 

. . . ,  n.  If  there  exists  an  i*  such  that  V  -  e^*  is  feasible 

i  k+1 

and  f(V  -  e^*)  <  f(X  ),  set  X  -  V  -  e^  and  go  to  Step  1. 
Otherwise  go  to  Step  2. 

1c 

Step  2.  Test  of  optimality.  Test  the  points  X  +  e^,  i  -  1,  ...»  n. 

If  there  exists  an  i*  such  that  X^  +  e^*  is  feasible  and  f(X^  +  e^*)  <  f(X^) 
k+1  k 

set  X  -  X  +  e^*  and  go  to  Step  1,  otherwise  terminate. 

The  Modified  Algorithm  of  Brooks  et  al. 

This  approach  requires  the  objective  function  f(X)  be  approximated 
by  a  separable  functiod  of  the  form 

n  m 

f (X)  -  Z  Eb.F.(x.jj) 
i-1  j-0  J 

where 

bj  -  C^x',  x’,  xp  >0 

for  some  X’  -  (x|,  x^,  x^) ,  and  for  some  function  3  (  .  ) . 

m 

Let  gi<x.)  -  Z  b  F  (x.;  j) 

j*0  J  1  1 

then 

f  (X)  -  Z  g  (x  ) 
i-1  1  1 


When  g1(X1>  is  strictly  decreasing,  then  the  implementation  of  marginal 
analysis  is  justified,  and  leads  to  the  following  series  of  steps. 


Step  0.  Initialization,  pick  b^,  b,,  ....  b  ,  set  k  ®  1. 

u  i  m 

Step  1.  Solve  the  problem. 


n 

min  Z  gi(x1> 
i*l 


Subject  to 


m 


Z  c,x 
i=l  1 


i±b 


_>  0  integer 

by  the  following  algorithm  which  is  known  as  marginal  analysis  [4]. 

Step  1  i)  Initialization:  X°  *  (0,  0,  ...»  0),  set  l  =  0. 
P+1  P 

Step  1  ii)  Y  =  Y  +  e^ ,  if  feasible,  where  j  is  any  index 
for  which  the  ratio 


gt(yi  +  1)  -  g^y^ 


is  maximum. 

P+1  k  P, 

Step  1  iii)  Stop  if  CY  >  b,  and  set  X  =  Y  ,  then  go  to  Step  2. 
Otherwise  go  to  Step  1  ii. 

Step  2  Calculate  b^  *  (X^)  for  j  =  0,  1 . m. 

k  k-1 

Step  3  If  bj  is  very  close  to  b^  for  all  j  ■  0,  1,  . . . ,  m, 

JL 

then  X  =  X  ,  and  stop.  Otherwise,  set  k  =  k  +  1,  and 
go  to  1. 

This  is  a  modified  version  of  the  algorithm  proposed  by  R.  B.  S. 
Brooks,  C.  A.  Gillen,  and  J.  Y.  Lu  on  page  24  of  [3].  The  Step  1  of  their 
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t 


approach  was  quite  different  from  the  one  described  above.  They  attempted 
to  find  the  saddle  point  solution  to  the  problem 


max  min  L(X:X)  ■  Eg^(x^)  +  X(£c^x^  -  b) . 

X>0  X>0 

Their  Step  1  was  as  follows: 

Step  1)  Find  x, ,  x»,  . . . ,  x  at  k  iteration  which  maximizes 
l  i.  n 

m  n  n 

Y.  b  Elog  F  (x  +  j)  -  X  I  c  x 
j*0  4=1  11  i=l 

Since  this  function  is  separable,  the  optimization  is  carried  out  for 
each  component  separately.  At  the  end  of  the  third  step,  the  Lagrange 
Multiplier  is  rearranged  as  suggested  in  [2],  and  the  algorithm  repeated 
again. 

A  Saddle  Point  Search  Algorithm 

This  algorithm  solves  the  problem 

max  min  L(X;X). 

X>0  Xenn 

The  algorithm  in  fact  searches  the  saddle  point  (X*,  X*)  of  the 
Lagrangian  function.  This  is,  however,  more  complicated  than  the 
previous  algorithms,  and  requires  subalgorithms  to  solve  the  primal 
and  dual  subproblems.  This  was  originally  suggested  by  H.  Everett  [2]. 
The  procedure  is  as  follows: 

Step  0.  Initialization.  Pick  X°  and  X  .  Set  k  =  0. 

o 

Step  1.  Solve  the  primal  subproblem. 

min  L(X;X.  ) 

xen" 

k+1 

to  obtain  X*  . 
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Step  2.  Solve  the  dual  subproblem 


(  max  LCX^Sx) 

A>0 

to  obtain  Aj^. 

Step  3.  Stop  if  is  very  close  to  A^.  Otherwise  set  k  -  k  +  1 

and  go  to  Step  1. 

B.  L.  Miller  developed  an  algorithm  to  solve  the  primal  subproblem 
for  discretely  convex  functions.  Since  a  detailed  description  of  this 
algorithm  is  presented  in  [6],  it  is  not  going  to  be  discussed  here. 

This  Is  a  very  complicated  algorithm.  The  author  claims  it  produces 
good  solutions. 

At  the  kth  iteration  (k  _>  1) ,  a  piecewise  convex  approximation  of 
L(X;A)  can  be  given  in  the  space  generated  by  the  points,  X®,  X*,  ...» 
X^+^  as  follows: 


L(X;A) 


l  y  L(X  ;A) 

t-0  c 


k+1 

where  I  y  -  1,  y  >  0  for  all  t  =  0,  1 . k+1. 

t-0 


Then  the  dual  subproblem  is  restated  as 


max 


Z  y.UX11;!) 
t-0  c 


Subject  to 


k+1 

i  yt 

t=o  1 


l 


A>0,  yfc  >.  0  for  t-0,  1,...,  k+1 
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Since  L(X;A)  ■  f(X)  +  A  CX  -  Ab,  then  the  problem 


max 


k+1 

E 

t»0 


k+1 
+  A  E 

t-0 


ytcx 


Ab 


k+1 

Subject  to  E  y  *  1 

t*0  c 


A>0,  Ut  2.  0  ^or  t  *  0,  1,...,  k+1 
can  be  considered  as  the  Lagrangian  of  the  linear  program 


k+1 

min  E  y 
t=0 


f(Xt) 


Subject  to 


k+1 

1 

t-0 


CX 


<  b 


k+1 


1 


Ut  >  0  for  t  «  0,  1,...,  k+1 

or  as  the  Lagrangian  of  the  dual  linear  program 

min  (Ab  +  v) 

Subject  to  ACXf  +  v  >  f(XC)  for  t  =  0,  1,...,  k+1 

A  >  0 

v  unrestricted 

Therefore  the  dual  linear  program  can  be  solved  Instead  of  the  sub¬ 


problem  in  Step  2.  This  can  be  done  as  suggested  by  Brooks  R.  B.  S.  and 
A.  M.  Geoffrion  in  [2]. 


Example  Problem  Solutions 
General 

A  war  readiness  spares  kit  problem  with  five  Items  Is  used  to  demon¬ 
strate  the  calculations  In  the  algorithms.  Computer  program  listings  are 
given  In  Appendix  A.  The  steady-state  Poisson  demand  rates  and  cost  of 
each  Item  is  given  in  Table  4.1. 

The  budget  is  assumed  to  be  M  ■  $25000.  In  addition  the  value  for 
m  of  Equation  is  chosen  to  be  5.  Also,  °  1  for  1=1,  2,  3,  4,  5. 

TABLE  4.1 

COSTS  AND  DEMAND  RATES  FOR  THE  SAMPLE  WRSK  PROBLEM 


Solution  by  Means  of  the  Primal  Search  Method 
Iteration  0.  Initialization:  Parameter  p  is  chosen  to  be  .85  and 
the  Lagrange  multiplier  is  specified  as  1  =  .0002.  (For  further  reference 
see  [9].)  This  iteration  seeks  an  X°  which  satisfies  the  condition  of  3.11 
page  4.6. 

Since  -1  +  F^l)  +  (.0002)  (2980)  =  -1  +  .37962  +  .596  =  -  .02438, 
and  -1  +  Fx(2)  +  (.0002)  (2980)  =  -1  +  .64963  +  .596  =  .24563,  xj  =  2. 
Repeating  the  same  argument  for  each  component,  yields  X°  =  (2,  2,  3,  6,  6) 
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Iteration  1.  Seeking  a  better  feasible  solution:  The  first 
differences,  Vf(X^),  at  the  point  X°  *  (2,  2,  3,  6,  6)  are  computed 
(by  Step  1,  page  4.6),  and  Af  (2,  2,  3,  6,  6)  -  (-.240,  -.106,  -.015,  -.163, 
-.035) 

P  *  ^12,-2  3,  6,  _6)_  c  =  >00008 

Since  p  *  .85,  then  after  normalization,  the  direction  is 
d2  =  (.610,  -.225,  -.268,  1.,  .190) 

and  a  is  evaluated  to  be  1.  Thus  a  point  U  =  (3,  2,  3,  7,  6)  is  obtained 
by  finding  the  nearest  integer  point  to  Y  =  (2.610,  1.775,  2.732,  7., 

6.190).  However  the  point  U  is  not  feasible,  cost  (U)  -  26398,  hence  a 

search  in  the  negative  direction  is  required.  For  that  purpose 

A1f(2,  2,  3,  7,  6)  =  -.262,  A2f(3,  1,  3,  7,  6)  =  -.333,  A3f(3,  2,  2,  7,  6) 

-  -.076,  A4f(3,  2,  3,  6,  6)  =  -.185,  and  A5f(3,  2,  3,  7,  5)  -  -.100  are 
calculated.  The  minimum  of  these  values  is  for  i  *  3,  but  (3,  2,  2,  7,  6) 
is  not  feasible.  The  second  minimum  value  is  for  i  »  5,  again  (3,  2,  3,  7, 
5)  is  not  feasible.  The  third  minimum  is  for  i  -  4,  fortunately  (3,  2,  3, 

6,  6)  is  feasible  and  f(3,  2,  3,  6,  6)  *=  .98571.  Therefore  X1  -  (3,  2,  3, 

6,  6)  with  CX  «  24898. 

Iteration  2.  As  in  Iteration  1,  a  search  direction  d2  ■  (-.356,  .362, 
-.037,  1.,  .216)  issuing  from  (3,  2,  3,  6,  6)  is  determined.  A  point 
U  -  (3,  2,  3,  7,  6)  is  obtained  similarly.  Since  U  is  infeasible,  and  a 
search  in  the  negative  direction  does  not  yield  a  distinct  and  better 
feasible  point  than  X^  *  (3,  2,  3,  6,  6),  a  search  in  the  positive  direction 
at  X^  is  employed.  All  the  points  of  the  form  X*  +  e^  for  i  ■  1,  2,...,  5 
are  infeasible.  Hence  ■  (3,  2,  3,  6,  6)  is  optimal.  This  small  prob¬ 


lem  takes  .24  seconds  of  execution  time. 


Solution  by  the  Modified  Solution  Method  of  Brooks  et  al. 


Iteration  1.  Choose  b.  -  b,  «■  ...  •  b.  “0  and  b_  =  1.  Set  x?  «  0 

0  1  4  5  i 

for  i  ■  1,  2,...,  5,  then  g^(x^)  *  b^  log  F^(x^  +  5),  and  calculating  the 
gi<x.  +  1)  -  g  (x  ) 

ratio  r.  - - r -  for  each  i  =  1 .  5  yields:  r.  =  0.00000496, 

i  Ci  1 

r,  *  0.00000202,  r,  *  0.00000270,  r,  -  0.00014200,  rc  *  0.00024952.  For 

2  3  4  5 

instance. 


r^  is  the  maximum,  hence  x2  =  1.  Updating  r^  yields:  r^  =  0.00000496, 

r2  =  0.00000202,  r3  =  0.00000270,  r 4  =  0.00014200,  r5  -  0.00011714.  The 

maximum  is  r^,  then  xj  ■  1.  Continuing  in  this  fashion  produces  X*  = 

(2,  2,  2,  8,  7)  with  costCX1)  ■  24801,  and  calculating  b  by  b^  »  if 

J  J  t-x 

F±(x*  +  J)  gives  bj  -  .419,  b*  =»  .726,  b*  =  .898,  bj  -  .967,  bj  -  .991, 
and  bj  =  .998. 

Iteration  2.  Repeating  the  same  steps  with  b^  j  "  0,  1 .  5  yields 

X2  =  (2,  2,  4,  7,  9)  with  cost(X2)  =  24915,  and  b2  =  .450,  b2  =  .728, 

b2  *  .891,  b2  =  .962,  b2  =  .988,  and  b2  *  .996.  Since  b^'s  are  not  within 
2 

.001  of  bj's»  8°  to  Iteration  3. 

2 

Iteration  3.  Repeating  the  same  steps  again  with  b^  for  j  =0,  1, 
2,..., 5  yields  X^  3  X2  and  b^  *  b2,  hence  X2  is  optimal  with  f(2,  2,  4, 

7,  9)  -  .98619. 

This  solution  however  is  a  near  optimal  solution.  The  primal  search 
procedure  has  produced  a  better  solution  than  this  one.  On  the  other 
hand  this  method  is  faster,  the  execution  time  to  solve  this  small  problem 
with  this  algorithm  is  .15  seconds. 
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In  this  algorithm  b^'s  are  treated  as  parameters.  They  have  practical 
implications  for  real  life  problems.  Bj  is  given  by  the  product  F±(X  +  j) 
for  some  X  and  j  *  0,  1,...,  m.  This  product  is  in  fact  equal  to  the  prob¬ 
ability  of  cannibalizing  at  most  j  airplanes  with  a  kit  X.  In  other  words 
bj  gives  the  operational  rate  given  j  airplanes  for  cannibalization  for  a 
kit  X.  Selecting  b^  “  b^  =  ...  *  ^®-l  *  ®  an(*  b^  *  1  implies  that  all  the 
planes  will  be  cannibalized  and  corresponds  to  the  most  pessimistic  view. 


Solution  by  the  Saddle  Point  Search  Algorithm 

Iteration  1.  A  starting  solution  of  X^  =  (10,  10,  10,  10,  10)  with 

f(X°)  ■  .02486  and  cost  (X^)  =  70380  and  a  Lagrange  multiplier  of 

An  -  .0001  is  chosen.  In  the  first  step  the  subproblem  min  L(X;An)  = 

XeIIn  0 

f (X)  +  Aq(CX  -  M)  is  solved  by  the  algorithm  proposed  by  B.  L.  Miller  in 
[6].  This  algorithm  produces  the  solution  X1  «  (2,  2,  2,  6,  6)  with 
f(XX)  -  1.28241  and  costCX1)  -  21456. 

Step  1.  Initialization:  An  initial  point  is  obtained  by  solving 
-1  +  T±(a  -  1)  +  A  C±  £0  and  -1  +  F^)  +  A  C±  >  0  and  setting  X±  = 
max  (0,  u>±  -  k)  for  all  i-1,  ....  n  where  k  is  an  input  parameter  and 
chosen  as  1  in  this  particular  example.  Since  -1  +  F^(2)  +  (.0001) 
(2980)  *  -1  +  .6496  +  .2980  =  -.0524  and  -1  +  Fj(3)  +  (.0001) (2980)  * 

-1  +  .8386  +  .2980  -  .1366,  then  Wj  -  3,  consequently  x°  =  2.  Repeating 

o  -  O  ft  O  / 

the  same  argument  for  each  component  yields  that  x^  <=  2,  x^  “Z,  “o» 

and  Xj  ■  6. 

Step  2.  Finding  A  Stationary  Vector:  Starting  from  X°  -  (2,  2, 

2,  6,  6),  a  point  X  is  sought  such  that  ^  g(X  -  e^  <_  0  and  g(X)  >  0 
for  all  i“l,  -1 •  For  instance  for  item  5,  A^  g(2,  2,  2,  6,  6)  is 
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calculated  as 

follows; 

A=  g(2,  2, 

2,  6,  6)  =  X  C5  +  F5(6  +5+1) 

-  1. 

5 

-  I 
j=*0 

(F5<6  +  j 

+  1)  -  F  (6  +  j)) 

3  i 

4 

ir  FjCXi  +  j) 

L=1 

> 

i 

F1(2+j) 

F2(2+j> 

F3<2+j) 

FA(6+j)  F5(7+j)  -  F5(6+j)  Product 

0 

.64963 

.80885 

.87949 

.76218  .  .03855 

.01358 

l 

.83864 

.93436 

.96623 

.86663  .01687 

.01107 

2 

.93787 

.98142 

.99225 

.93191  .00656 

.00558 

3 

.97955 

.99554 

.99850 

.96817  .  .00230 

.00217 

4 

.99414 

.99907 

.99975 

.98630  .00073 

.00072 

5 

.99851 

.99983 

.99996 

.99455  .00021 

.00020 

SUM  = 

.03332 

Then  a 

5  g(2,  2, 

2,  6,  6) 

=  .0345  + 

.99992  -  1.  -.03332  «  .00110. 

Since 

it  is  positive,  =  6.  This  argument  is  repeated  for  each  item  to 
obtain  a  stationary  vector  X  =  (2,  2,  2,  6,  6)  with  f(X)  =  1.28241, 
cost  (X)  -  21456,  and  g(X)  =  3.42801. 

Step  3.  A  heuristic  approach  to  obtain  improvement  in  "up-phase": 
n 

Let  S4  =  £  e^  -  e^,  then  the  ratios  ^  g(X  +  S^/a.^  g(x)  are  ranked  in 
k=i 

ascending  order,  and  the  vectors  of  the  form  X  +  will  be  tested  for 
improvement  where  is  the  vector  with  +1  for  variables  ranked  from 
1  to  i  and  0  for  the  ones  ranked  i  +  1  to  n.  Calculation  of  the  ratios 
gives  a^  g(X  +  gOO  =  —5 .61.  a2  g(X  +  S2)/  ^2  = 

A3  g(X  +  S3)/a3  g(X)  =  1.47,  A^  g(X  +  SA)/A4  g(X)  =•  1.44  and  A$  g(X  + 

S^) /a 5  gOO  =  1.64.  The  ordering  according  to  the  magnitudes  yields 
1,  2,  4,  3,  5  respectively.  For  instance,  for  i  =  4  =  (1,  1,  0,  1,  0) 

and  X  +  =  (3,  3,  2,  7,  6),  but  this  point  does  not  yield  a  better 

solution,  since  g(X  +  U^)  =  3.48107  is  greater  than  g(X)  =  3.42801. 
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Repeating  the  same  argument  for  each  component  results  that  the  current 
solution  is  stationary. 

The  same  argument  is  repeated  for  the  "down-phase"  too.  This  does 
not  yield  a  better  solution  either. 

Step  4.  A  Second  Initialization  and  Reducing  the  Set  of  Potential 
Improving  Points. 

First  the  up-phase  is  implemented.  For  instance  for  i  *  1 ,  the 

vector  (0,  1,  1,  1,  1)  will  be  checked  fcr  improvement  by  looking  at 

g(2 ,  3,  3,  7,  7).  Since  ^  g(2,  3,  3,  7,  7)  -  -.028  is  negative, 

the  test  fails.  For  i  =  2,  the  vector  (1,  0,  1,  1,  1)  is  tested  and 

* 

&2  8(3*  2,  3,  7,  7)  =  .005,  then  the  conclusion  is  that  can  not 
3.  Next  the  vector  (1,  O',  0,  1,  1)  is  used  for  the  third  component  and 
A3  g(3,  2,  2,  7,  7)  =  -.056.  Similarly  g(3,  2,  3,  6,  7)  =  -.069  and 
A5  g(3,  2,  3,  7,  6)  =  -.022.  The  test  fails  for  i  *  1,  3,  4,  and  5. 

Consequently  Part  1  is  used.  Calculation  of  the  ratios  Aj  g(2,  3,  3, 

7,  7)/Aj  g(2 ,  2,  2,  6,  6)  -  -5.61,A3  g(3,  2,  2,  7,  7)/a3  g(2,  2,  2,  6, 

6)  =  1.37,  Aa  g(3,  2,  3,  6,  7)/a4  g(2,  2,  2,  6,  6)  *  1.33  and  finally 

A^  g(3,  2,  3,  7,  6)/A3  g(2,  2,  2,  6,  6)  =  1.52  gives  that  the  maximum 

ratio  is  for  i  «  5.  By  contradiction  it  will  be  proved  that  *  1 
implies  s3  ■  1  where  X  +  S  is  the  point  to  be  tested.  Since  b^  =  -XC3 

-&3  g(2,  2,  3,  6,  7)  =  .042  >  0,  S3  can  not  be  1.  Repeating  this 

argument  for  i  »  1  and  i  =  4  does  not  yield  a  better  solution.  Conse¬ 
quently  the  down  phase  is  implemented.  Down-phase  does  not  generate  a 
better  feasible  solutior\.  Finally  the  conclusion  is  that  X  =  (2,  2,  2, 

6,  6)  is  a  stationary  point.  ' 


Then  the  linear  program 


min  .02486  u  +  1.28241  y2 
Subject  to  70380  ^  +  21456  ^  <  25000 

ux  +  w2  =  1 

Ml»  V2  >  0 

is  solved  to  obtain  X^  =  .0000257.  Since  X^  is  not  within  .000005  of 

XQ  another  iteration  is  performed. 

Iteration  2.  The  subproblem  min  L(X;X1)  =  f(X)  +  X. (CX  -  M)  is 

XelP  l 

2 

solved  again  in  the  first  step  by  the  same  method  to  obtain  X  **  (4,  4, 

4,  9,  9)  with  f (X2)  =  .17727  and  cost(X^)  =  37377.  In  the  second  step 
the  linear  program  min  .02486  y^  +  1.28261  y^  +  ^ 17727  y^. 

Subject  to  70330  y  +  21456  y2  +  37377  y3  <  25000 

yL  +  y2  +  y3  =  1 
>  0  i  =  1,  2,  3 

is  solved  to  obtain  X 2  =  .0000694.  Since  X2  is  not  within  .000005  of  X^, 
another  iteration  is  required.  Continuing  in  this  fashion  produces  the 
following  results: 

Iteration  3.  X3  =  (3,  3,  3,  8,  7)  with  f(X3)  =  .49955  and  cost 
(X3)  =  29994.  Further  X3  =  .0000917. 

Iteration  4.  X4  =  (2,  2,  3,  7,  7)  with  f(X4)  =  1.02459  and  cost 
(X4)  =  23763.  Also  X^  =  .0000848. 

Iteration  5.  X3  =  (3,  2,  3,  7,  7)  with  f(X3)  =  .75556  and  cost 
(X5)  =  26743.  Also  X$  =  .0000903. 


Iteration  6.  X  =  (2,  2,  3,  7,  6)  with  f(X  )  =  1.06282  and  cost 


f(X(Y)) 


(X6)  =  23418  and  Xfi  =  .0000903. 


This  algorithm  produced  a  solution  with  a  better  objective  value 
than  the  previous  ones,  but  it  is  an  infeasible  solution.  In  fact, 
this  algorithm  suffers  from  the  duality  gap.  In  the  first  step  of 


this  algorithm,  the  subproblem  min  L(X;X)  is  solved.  Since  L(X;X)  is 

Xen” 

discrete  in  X  e  11°,  there  is  a  range  for  X  values  which  produce  the 


same  answer.  From  the  results  of  previous  algorithms,  it  was  dis- 


covered  that  the  optimal  solution  was  X  =  (3,  2,  3,  6,  6).  Evidently 


its  corresponding  Lagrange  multiplier  should  lie  between  X^  =  .0000848 


and  X_.  =  .0000903.  Unfortunately  the  algori  thm  won'  t  generate  this 
value,  and  if  X  is  close  to  X^,  the  solution  to  subproblem  is  X  and 
if  X  is  close  to  X^,  then  the  solution  is  X6,  see  Figure  4.4. 


Figure  4.4.  Duality  Gap 


Selection  of  Solution  Method 


\ 

I 

The  objective  of  this  research  is  to  develop  and/or  modify  an  efficient 
solution  method  to  solve  the  WRSK  problem.  Therefore,  these  three  solution 
|  methods  have  been  under  investigation  to  observe  their  computational  require¬ 

ments.  The  result  of  this  investigation  shows  that  the  Modified  Algorithm 
of  Brooks  et  al.  is  more  efficient  in  terms  of  computational  time  than  the 
others.  This  algorithm  is  also  going  to  be  compared  with  the  current  D029 
Model  which  is  summarized  on  the  next  page.  The  result  of  this  comparison 
will  be  discussed  in  Section  V. 

The  Modified  Algorithm  of  Brooks  et  al.  is  a  greedy  algorithm  in  nature, 
and  will  be  refined  to  in  the  remainder  of  this  report  as  the  "Greedy  Algorithm". 
It  uses  marginal  analysis  as  a  subalgorithm.  It  is  easier  to  implement  in 
solving  the  WRSK  problem  and  is  faster  than  the  others.  The  selection  of  its 
initial  parameters  does  not  require  any  guesswork  from  the  user.  The  parameter 
values  obtained  for  the  requirement  kit  57-18  can  serve  that  purpose.  On  the 
other  hand,  marginal  analysis  adds  one  unit  of  the  item  which  will  yield  the 
highest  rate  of  return  in  terms  of  performance  measures  per  dollar  invested 
in  the  kit.  Approximation  of  the  E(NORS/X)  function  by  a  separable  function 
[as  discussed  in  Part  III]  eases  the  computations  and  in  fact  eliminates  the 
need  for  updating  all  of  the  marginal  returns  whenever  a  unit  of  some  item 
is  added  to  the  kit  as  is  necessary  in  the  current  D029  computer  program.  The 
speed  of  convergence  of  an  algorithm  is  very  important,  especially  when  solving 
large  size  problems,  and  this  algorithm  supplies  this  advantage. 

The  others  have  some  particular  advantages,  but  they  are  slower  than 
this  algorithm.  The  Saddle  Point  Search  Algorithm  is  theoretically  very 
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appealing  but  suffers  from  the  duality  gap  as  indicated  by  the  small  example 
of  the  last  section.  The  Primal  Search  Algorithm  uses  the  first  differences 
in  determining  the  search  direction  to  explore  better  feasible  solutions. 

All  three  algorithms  contain  some  heuristic  rules.  Thus  the  solutions  are 
not  necessarily  the  true  optimum  which,  for  the  same  reason,  is  true  for  the 
current  D029  Model.  The  comparison  of  their  solutions  either  with  each  other 
or  with  the  D029  solution  indicates  that  they  can  generate  good  near  optimal 
solutions. 


The  Current  D029  WRSK  Model 


This  model  and  its  associated  computer  program  is  currently  in  use 


by  WPAFLC.  This  solution  method  also  makes  use  of  marginal  analysis.  A 


detailed  description  of  this  algorithm  is  provided  in  [8],  A  brief 
description  of  the  algorithm  is  provided  here.  The  algorithm  solves  the 


problems  of  the  form 


I 

s.  t. 


Min  CX 

E(NORS/X)  £  a1 
E(SDO/X)  <  a2 
X  £  0  integer 


in  the  following  manner. 

Step  0.  Set  U  -  0,  X°  *  0 

k  k 

Step  1.  Choose  a  W  and  1^ 


4.22 


Ir  Ir 

Step  2.  Add  1^  units  to  the  kit  by  selecting  the  highest  I^t  ratios 


from 


E(NORS/Xk)  +  E(SDO/XU) 


k+1  k 

X  =  X  +  X  ei 


ieH, 


Nt 


k  1c 

where  H  is  the  index  set  of  the  highest  I  ratios. 

*Nt  Nt 

k+1  k+1  k+1 

Step  3.  Calculate  E(NORS/X  )  E(SDO/X  )  and  CX 

k+1  x 

Step  4.  Phase  1:  Terminate  if  CX  >  b 

where  b*  is  90%  of  the  cost  of  Kit  57-18 

tf+i 

Phase  2:  Terminate  if  E(N0RS/X  )  <  a.  and  E(SD0/X  )  <  a 


Otherwise  Set  k  =>  k  +  1  and  go  to  Step  1. 


PART  V 


A  SUMMARY  OF  COMPUTATIONAL  EXPERIENCE 

Introduction 

Computational  experimentation  was  devoted  to  the  investigation  of 
computational  requirements,  initialization  of  parameters  for  the  Creedy 
Algorithm  (Modified  Algorithm  of  Brooks  et  al)  and  the  IBM  version  of 
D029  Program.  F4D  Data  was  used  in  the  computational  analysis.  Some  of 
the  subroutines  remain  the  same  in  both  programs  in  order  to  make  mean¬ 
ingful  comparisons.  (The  FORTRAN  listings  of  the  programs  are  in  Appen¬ 
dix  B.2  &  C.2.)  The  Greedy  Algorithm  takes  16.53  seconds  of  CPU  time, 
and  the  IBM  version  of  D029  Program  takes  24.78  seconds  of  CPU  time  on 
IBM  370/168  at  O.S.U.  The  Greedy  Algorithm  is  obviously  faster  on  the 
F4D  problem  than  the  algorithm  currently  in  use. 

The  Parameters  bj  and  m  of  the  E(N0RS/X)  Approximation 

The  Greedy  Algorithm  requires  that  bj  and  m  be  Initialized  in  ad¬ 
vance.  A  thorough  inspection  of  the  WRSK  problem  suggested  that  initial 
b°j  values  can  be  calculated  from  the  requirement  kit  57-18  by  the  formula 

b°j  =  tt  F^Cxi  +  ja^)  j  =  0,1,2,...  K  =  Number  of  planes  and 

i=l  n  =  Total  number  of  items  in  consideration 

in  the  squadron,  where  x^  is  the  kit  57-18  quantity  for  item  i.  The  ini¬ 
tial  b°j  values  for  Phase  1  and  Phase  2  of  the  WRSK  Problem  are  displayed 
in  Table  5.1.  Each  bj  is  the  operational  rate  given  that  there  are  j  air¬ 
craft  for  cannibalization.  Hence,  b°j's  represent  operational  rates  under 
the  kit  57-18.  These  are  good  choices  for  bj's  since  the  b°j's  may  be  con- 


sidered  as  operational  rate  goals  as  a  consequence  of  the  AFLCR.  Since  the 
WRSK  problem  attempts  to  find  a  kit  which  will  perform  better  than  the  kit 
57-18,  these  b°j's,  which  can  be  observed  in  Figures  5.7a  and  5.7b  on  pages 
23  and  24  respectively,  are  the  target  operational  rates  for  the  optimal  kit. 
On  the  other  hand,  a  natural  selection  for  m  is  the  number  of  planes  in  the 
squadron.  An  inspection  of  the  Table  5.1  indicates  that  the  bj 's  approach 
1  for  some  lower  value  than  m.  Consequently,  using  Table  5.1,  one  can  choose 
t.i  =  12  for  Phase  1  and  still  get  the  same  answers. 

The  b*j's  displayed  in  Table  5.1  correspond  to  the  calculated  values 
from  the  optimal  kit  and  can  also  be  observed  in  Figures  5.9a  and  5.9b  respec¬ 
tively.  For  Phase  1  the  b*j's  are  very  close  to  the  b°j's.  However,  a  simi¬ 
larly  close  relationship  between  b*j  and  b°j  does  not  occur  in  Phase  2.  Note 
that  each  b*j  b°j  for  all  j.  Since  the  b°j  represent  the  operational  rate 
goals,  this  shows  that  the  optimal  kit  provides  higher  operational  rates  than 
the  goals. 

Instead  of  hueristically  choosing  the  initial  b°j's,  the  Greedy  program 
of  Appendix  B.2  uses  the  b°j’s  of  the  5718  kit.  It  also  uses  m=24,  the  number 
of  planes  in  the  squadron,  since  24  is  not  a  large  number. 


5.2 


TABLE  5.1 


COMPARISON  OF  bj  VALUES 


Phase  1 


5718  Kit 
b°j 

.00000010 
.00286799 
.06779460 
.25850234 
.48975927 
.68018747 
.81107932 
.89285344 
.94108090 
.96842135 
.98345607 
.99151443 
.99573534 
.99789912 
.99898552 
.99951989 
.99977739 
.99989891 
.99995505 
.99998044 
.99999167 
.99999653 
.99999859 
.99999944 
1  . 


Optimal  Kit 
b*j 

0.0 

.00172540 

.06441003 

.26115723 

.49893585 

.69276720 

.82435903 

.90472555 

.95038868 

.97499286 

.98772442 

.99410652 

.99722461 

.99871559 

.99941528 

.99973806 

.99988454 

.99994995 

.99997868 

.99999108 

.99999634 

.99999852 

.99999942 

.99999972 

.99999991 


Phase  2 


5718  Kit 
b°j 

0. 

I'. 

.  OOOOOOf.'  1 
.00011449 
.00861967 
.07484124 
.23535679 
.44260779 
.63162230 
.77280201 
.86636330 
.92392780 
.95769561 
.97688397 
.98754667 
.99337298 
.99651327 
.99818555 
.99906608 
.99952462 
.99976077 
.99988101 
.99994153 
.99997163 
1 . 


Optimal  Kit 
b*j 

0.0 

0.0 

.00000935 
.00376217 
.05835153 
.22116812 
.43810253 
.63279580 
.77581119 
.86953349 
.9266795* 
.95987318 
.97850019 
.98868473 
.99413854 
.99700800 
.99849417 
.99925267 
.99963432 
.99982364 
.99991621 
.99996080 
.99998195 
.99999182 
•  '■9J9h36 


Articulation  of  the  Weight  w 


The  weighting  factor  w  represents  the  relative  importance  of  E(SDO/X) 
function  compared  to  the  E(NORS/X)  function.  The  weighting  factor  is,  in 
fact,  a  measure  of  cannibalization  of  any  plane  for  each  back-ordered  part. 
Conversely,  1/w  measures  the  number  of  consolidated  parts  from  a  cannibilized 
airplane.  During  a  mission,  if  a  plane  is  grounded  for  lack  of  spare  parts, 
it  may  then  be  cannibalized  to  repair  others  by  removing  its  good  parts  and 
installing  these  parts  in  other  downed  planes.  Theoretically  100%  cannibalized 
is  possible.  This  process  increases  the  work  load  during  the  mission.  Obvi¬ 
ously,  a  reduction  in  the  work  load  is  required  by  the  mission  command.  In 
practice,  however,  there  is  empirical  knowledge  as  to  when  to  cannibalize  a 
plane,  which  plane  to  cannibalize,  and  how  many  parts  to  be  removed  from  a 
plane.  In  the  D029  program,  manipulation  of  the  parameter  RATIO  indicates 
that  this  empirical  knowledge  is  regularly  applied  (RATIO  in  the  Creedy  pro¬ 
gram  has  a  different  meaning).  It  has  also  been  observed  that  RATIO  =  50  pro¬ 
duces  better  results  than  any  other  value,  for  the  purpose  of  reaching  the 
E(NORS)  and  E(SDO)  goals  simultaneously.  The  operational  meaning  of  RATIO  = 

50  is  that,  on  the  average,  50  parts  are  removed  from  a  cannibalized  plane; 
thus  RATIO  corresponds  to  1/w.  This  would  also  mean  that  2%  of  an  airplane 
is  cannibalized  for  each  back  order.  It  is  suggested  by  this  research  that 
w  be  assessed  in  advance  by  evaluating  this  empirical  knowledge.  Tables  5.2, 
5.3,  5.4  and  Figures  5.1,  5.2,  5.3  show  the  relationship  between  E(N0RS/X) 
and  E(SD0/X)  for  different  weights.  These  relationships  based  on  limited 
experimentation,  appear  to  be  linear  effects. 


The  problem  is  to  determine  the  best  effective  weight.  Computational 
experimentation  with  F4H  data  indicates  that  w  =  2.25%  is  the  effective 
weight  for  the  initial  six-day  period,  and  that  1.5%  ±  v  <_  2.5 %  gives  a  range 
for  the  effective  weight  for  a  thirty-day  period.  This  means  that  45  back 
ordered  parts  would  be  cannibalized  during  the  initial  six-day  and  40  to  65 
back  ordered  parts  would  be  cannibalized  during  the  30  day  mission  period. 

This  observation  should  be  combined  with  the  empirical  knowledge  to  come  up 
with  an  effective  weight  in  determining  a  best  compromise  kit. 

Only  a  limited  effort  was  undertaken  to  determine  the  underlying  re¬ 
lationships  between  the  weighting  factor  w,  budgetary  levels,  objective  func¬ 
tion  values,  and  kit  configuration  sensitivity  to  these  parameters  and  measures. 
Tables  5.2,  5.3,  5.4  and  Figures  5.1,  5.2,  5.3  indicate  various  ways  that  re¬ 
lationships  can  be  exhibited  and  analyzed.  Analysis  and  study  of  the  relation¬ 
ships  exhibited  lead  to  the  selection  of  parameters  and  measures  that  are  used 
in  the  next  section  which  compares  various  solutions  which  have  been  computed. 


TABLE  5.2 


OBJECTIVE  FUNCTION  fCw)  VS.  WEIGHTING  FACTOR  w: 
PHASE  1  SOLUTIONS  FOR  A  BUDGET  OF  $3,635,906 


w(10“2) 

E(NORS/X) 

E(SDO/X) 

Y(v) 

f(w) 

2 

4.80602 

31.02199 

5.42646 

5.4346 

2.25 

4.84981 

29.36157 

5.51045 

5.5061 

2.5 

4.85617 

29.07728 

5.58315 

5.579 

3. 

4.87135 

28.57339 

5.72855 

5.724 

4. 

4.94188 

26.69839 

6.00982 

6.013 

V  (w)  = 

E(N0R/X)  +  v  • 

E(SDO) 

i  <w)  = 

28.94220  w  +  4. 

85578 

0  0 


w(10"2) 


Figure  5.1.  Plot  of  f(w)  Exhibited  in  Table  5.2 


5.6 


TABLE  5.3 


OBJECTIVE  FUNCTION  'i'(w)  VS.  WEIGHTING  FACTOR  w: 
PHASE  2  SOLUTIONS  WITH  A  BUDGET  OF  $7.5  x  106 


w(10~2) 

E(NORS/X) 

E(SDO/X) 

nw> 

.3 

6.56880 

261.25884 

7.35258 

7.36159 

1. 

6.91434 

113.02191 

8.04456 

8.05923 

1.25 

7.01317 

104.48303 

8.31921 

8.30838 

1.5 

7.12255 

96.53860 

8.57063 

8.55754 

2. 

7.15833 

94.41359 

9.04660 

9.05585 

2.5 

7.22784 

90.83792 

9.49879 

9.55417 

'V(w) 

=  E(NOR/X)  +  w 

•  E(SDO) 

'J’(w) 

=  7.06260  +  99. 

66263w 

TABLE  5.4 


OBJECTIVE  FUNCTION  H^w)  VS.  WEIGHTING  FACTOR  w: 
PHASE  2  SOLUTIONS  WITH  A  BUDGET  OF  $9  x  106 


wqo Cil 

1 

E(NORS/X) 

4.93755 

E(SDO/X) 

T(w) 

n») 

41.56367 

5.35319 

5.37412 

1.5 

5.01991 

38.17230 

5.59249 

5.57405 

1.75 

5.05138 

36.31621 

5.68691 

5.67401 

2. 

5.08688 

34.68301 

5.78054 

5.77398 

2.5 

5.12059 

V(v) 

nw) 

33.45378 

=  E(NORS/X) 

=  4.97427  + 

5.95693 

+  w  •  E(SDO) 

39. 98540 (w) 

5.97390 

5.7 


Figure  5.2.  Plot  of  E(NORS/X)  vs.  E(SDO/X)  from  Table 


Table  5.3 
Relationships 


Table  5.4 
Relationships 


; 


;\  Budget  =  $9,000,000 


Comparison  of  the  Solutions 


I 

Comparison  of  solutions  was  undertaken  relative  to  the  two  phases  of 
the  WRSK  Problem.  F4D  Data  was  used  as  sample  data,  and  is  provided  in 
|  Appendix  C.  Three  sets  of  solutions  are  shown  in  Figures  5.4  a,  b,  5.5  a, 

b,  and  5.6  a,  b.  They  are  the  Phase  I  and  Phase  II  output  results  of  the 
Greedy  Algorithm,  the  IBM  version  of  the  D029  Program  and  the  D029-CDC  Pro- 
|  gram  of  WPAFLC.  The  FORTRAN  listings  of  the  first  two  algorithms  are  pro¬ 

vided  in  Appendix  B.2  and  C.2,  and  are  the  final  products  of  this  current 
research.  It  should  be  understood  that  only  limited  experimentation  has 
taken  place  with  both  of  these  programs. 

The  frequency  of  the  differences  in  item  quantities  can  be  of  interest 
in  discussing  the  kits  obtained  by  these  algorithms.  In  Phase  I,  there  are 
|  180  items  under  consideration.  Table  5.5  indicates  that  73%  of  the  kit  items 

(for  the  three  comparisons)  have  the  same  quantity  values  with  a  standard  de¬ 
viation  of  1.4%.  Also  the  maximum  difference  between  the  individual  kit  items 
|  does  not  exceed  2.  One  can  also  note,  for  example,  that  the  Greedy  Algorithm 

kit  has  40  and  43  less  parts  in  the  kit  compared  to  the  D029  solutions  respectively. 

In  Phase  II,  there  are  254  items  under  consideration.  Table  5.6  shows 
the  frequency  of  the  discrepancy  of  the  kits.  75%  of  the  individual  items 
in  the  kits  have  the  same  quantity  values  with  a  standard  deviation  of  3%. 

This  indicates  that  the  kits  are  in  close  agreement.  In  this  case,  the 
Greedy  Algorithm  kit  has  19  and  18  less  parts  in  the  kit  compared  to  the  D029 
solutions  respectively.  However,  note  the  difference  of  nine  that  occurs  in 
column  1  and  column  3. 

This  difference  occurs  on  item  number  249  in  Phase  II.  WPAFLC's  D029-CDC 

Program  (Figure  5.6b)  yields  x*249=l,  on  the  other  hand  the  Greedy  Algorithm 

5.10 

/ 

/ 
t 

/ 


TABLE  5 .  5 


FREQUENCY  OF  DIFFERENCES  OF  KITS  IN  PHASE  I 


- - 

Greedy-D029  CDC 

1 - 

Greedy- D029  IBM 

D029  IBM-D029  CDC 

8 

1 

1 

6 

31 

46 

14 

134 

;  131 

129 

7 

3 

31 

-40 

-43 

+5 

i 

TABLE  5.6 

FREQUENCY  OF  DIFFERENCES  OF  KITS  IN  PHASE  II 


Greedy-D029  CDC 

1 

j  Greedy-D029  IBM 

2 

I 

3 

3 

34 

37 

196 

195 

18 

19 

— 

i  3 

1 

IL_  .  . 

-19 

i 

-18 

D029  IBM-D029  CDC 


(Figure  5.4b)  and  the  IBM  Version  of  the  D029  Program  (Figure  5.5b)  both 
produce  x*249=10.  This  item  has  a  moderate  demand  rate,  ji249=^*29728,  is 
not  expensive  at  $653,  but  has,  however,  a  high  number,  15,  of  applications 
on  an  aircraft  (Appendix  D.2.2).  Since  U249  =  ^ *  29728,  probability  of  demand 
not  exceeding  1  is  P(d  <_  1)  =  0.0056,  but  probability  of  demand  not  exceeding 
10  is  P(d  10)  =  .878.  Addition  of  9  more  items  increases  the  probability 
considerably,  thus  providing  more  protection  from  back  ordering  item  249. 

The  same  reasoning  applies  to  other  items  which  have  different  kit  values. 

If  the  difference  in  item  quantities  is  small,  the  probabilities  are  close 
to  each  other.  Investigation  has  not  been  undertaken  to  determine  why  the 
CDC  version  yields  X*24g=l.  Different  approaches  to  calculating  the  Poisson 
probabilities  may  be  a  contributing  factor. 

The  Greedy  Algorithm  solution  was  obtained  in  16.53  seconds  of  CPU  time 
while  the  D029-IBM  solution  was  obtained  in  24.78  seconds  of  CPU  time.  This 
indicates  the  potential  of  a  significant  improvement  in  computational  time  if 
the  Greedy  Algorithm  is  coded  for  a  CDC  FORTRAN  compiler  and  used  rather  than 
the  D029-CDC  program  currently  in  use. 

On  Figure  5.6b  it  can  be  noted  that  the  D029-CDC  program  has  obtained  a 
kit  with  E(NORS)  =  8.127,  E(SDO)  =  149.65,  and  cost  =  $7,108,451.15.  This 
cost  was  used  as  a  budget  parameter  for  the  Greedy  Algorithm  solution  which 
is  shown  on  Figure  5.4b.  The  solution  obtained  was  E(N0RS)  =  8.037,  E(SDO)  = 
144.96,  and  cost  =  $7,108,477.  For  $4  less  expenditure  of  budget  funds,  the 
Greedy  Algorithm  obtained  a  kit  yielding  better  performance  measures  than  the 
D029-CDC  solution. 

Figure  5.5b  illustrates  a  D029-IBM  solution  for  the  same  budget  for 
$7,562,983  used  by  the  D029-CDC  solution  (Figure  5.6b).  In  this  solution 


E(NORS)  =  8.130,  E(SDO)  =  145.30  and  cost  =  $7,007,734.  The  budget  Is  not 
utilized  to  its  full  potential.  The  objective  of  the  D029  approach  (Part  IV, 
page  4.22)  is  to  meet  the  NORS  and  SDO  goals  in  a  way  that  minimizes  cost, 
and  under  utilization  of  a  budget  can  result  since  the  actual  budget  is  the 
last  item  of  priority  that  is  checked  in  the  algorithmic  scheme  of  the  D029 
WRSK  solution  method. 

The  Greedy  solution  exhibited  in  Figure  5.7b  for  the  same  budget  as  in 
Figures  5.5b  and  5.6b  ($7,562,983)  yields  E(NORS)  =  7.158,  E(SDO)  =  96.81, 
and  cost  =  $7,562,944.  All  but  $39  of  the  available  budget  is  utilized. 

These  greatly  improved  performance  measures  are  as  expected  since  more  money 
has  gone  into  the  kit. 

The  WRSK  Model  solved  by  the  Greedy  Algorithm  (Part  IV,  page  4.1)  is  not 
the  WRSK  Model  solved  by  the  D029  Algorithm.  For  this  reason  it  is  difficult 
at  times  to  make  direct  comparisons  between  D029  solutions  and  Greedy  solutions. 
It  should  be  recognized  again  that,  in  essence,  the  goals  (constraints)  to  be 
satisfied  by  a  D029  solution,  are  the  elements  of  the  objective  function  (which 
is  to  be  minimized)  of  the  Greedy  model.  In  like  manner,  the  Greedy  model  budget 
constraint,  which  is  not  to  be  exceeded,  is  the  objective  function  to  be  mini¬ 
mized  in  the  D029  model. 

Figure  5.8b  illustrates  an  attempt  to  develop  a  Greedy  solution  that  is 
relatively  comparable  to  the  D029-IBM  solution  of  Figure  5. 5b.  Observe  that  a 
budget  of  $7,075,000  was  used  leading  to  a  solution  in  which  E(N0RS)  =  8.173, 
E(SDO)  =  145.09,  and  cost  =  $7,074,907  which  is  within  1%  of  the  D029-IBM  solu¬ 
tion  of  Figure  5.5b.  Since  two  different  models  are  being  used,  answers  exactly 
equal  to  each  other  are  not  be  expected.  However,  it  might  be  hypothesized  that 
a  closer  stopping  criteria  on  the  b*j  values  might  continue  to  improve  the 


Greedy  solution. 


\ 


All  of  the  algorithms  contain  some  heuristic  rules;  therefore  the  answers 
exhibited  in  this  section  are  best  categorized  as  near  optimal  solutions.  Both 


1 


algorithms  yield  satisfactory  answers,  but  with  differences  discussed  in  this 
Part  V  and  in  Part  VI. 
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Figure  5.6b  D029  CDC  Program-A  Phase  2  Solution  for  Kit  5718  Budget 
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SUMMARY  OF  RESULTS  ] 

•: 

i 

■S 

; 

This  research  provides  a  new  and  more  general  model,  the  "multiobjective  J 

model,"  for  the  WRSK  problem.  This  multiobjective  approach  resulted  in  an 
improved  measure  of  supply  performance  which  relates  E(NORS/X)  and  E(SDO/X) 

{ 

.'I 

for  a  WRSK  kit  in  a  meaningful  way.  J 

Three  potential  solution  algorithms  for  this  new  "multiobjective  model"  j 

were  Investigated.  The  Saddle  Point  Algorithm  was  the  most  general,  but  suf-  .j 

'•j 

fered  from  the  duality  gap.  Hie  Primal  Search  Algorithm  made  special  use  of  J 

the  properties  of  the  objective  function,  but  lacked  flexibility  to  adjust  j 

the  search  direct iori~in  an  effective  way.  j 

A  third  method  proposed  by  R.  B.  S.  Brooks  in  [3]  and  subsequently  mod- 
ified  during  this  research  lead  to  the  development  of  the  "Greedy  Algorithm"  j 

as  the  solution  procedure.  This  algorithm  requires  that  the  objective  function 
be  approximated  by  a  separable  function.  This  research  has  accomplished  the  t 

required  separability  of  the  multiobjective  model.  This  allows  the  use  of 
marginal  analysis  in  a  manner  that  achieves  the  highest  rate  of  return  per 
unit  of  performance  measure  per  dollar  Invested  in  the  kit  as  items  are  added 
to  the  kit. 

Computational  experimentation  has  shown  that  the  Greedy  Algorithm  is 
faster  than  the  other  two  potential  algorithms.  An  important  result  of  this 
research  has  been  that  the  Greedy  Algorithm  program  is  faster  and  provides  a 
better  solution  than  the  current  D029  using  F4D  data. 

The  Greedy  Algorithm  has  several  advantages.  It  requires  only  one 
ratio  of  marginal  return  be  updated  whenever  one  unit  of  an  item  is  added 
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to  a  kit.  This  is  an  immediate  consequence  of  having  achieved  a  separable 
approximation  of  the  objective  function  of  the  new  multiobjective  model. 

This  contributes  to  the  ability  of  the  Greedy  Algorithm  to  yield  an  improved 
kit  compared  to  the  D029-CDC  kit,  and  in  shorter  time  when  properly  compared 
to  the  D029-IBM  program. 

A  disadvantage  of  the  Greedy  Algorithm  is  that  it  is  constructed  on  the 
convergence  of  the  bj  parameters.  The  generation  of  a  new  set  of  parameter 
values  from  a  previous  set  takes  one  full  iteration.  The  parameter  values  of 
this  new  set  are  compared  to  the  previous  set.  Such  iterations  and  comparisons 
are  repeated  until  the  parameters  are  sufficiently  close,  at  which  time  the 
algorithm  terminates.  Therefore,  the  last  iteration  is  used  only  to  conclude 
that  the  optimal  solution  is  as  good  as  the  next  to  the  last  solution.  Thus 
there  is  always  an  extra  iteration  that  does  not  improve  the  kit  configuration. 
Another  disadvantage  might  be  that  kits  for  intermediate  budgets  might  be  de¬ 
sired,  but  must  await  determination  of  the  optimal  bj  parameters.  Also  the 
weighting  factor  w  must  be  determined  in  advance,  and  additional  experimentation 
may  be  required  to  determine  the  appropriate  value. 

Another  product  of  this  research  has  been  the  modification  of  the  current 
D029  program  for  compatibility  with  IBM  FORTRAN  compilers.  This  effort  has 
Provided  some  insight  into  the  problems  with  the  current  D029  program.  The 
word  packing  and  unpacking  required  by  the  CDC  compilers,  in  order  to  save 
memory,  consumes  a  great  deal  of  time.  The  storing  of  some  probability  values, 
and  the  interpolation  of  others,  when  needed,  does  not  appear  to  be  a  good 
policy.  Interpolation  underestimates  the  probabilities;  hence,  the  marginal 
rates  of  return  are  less  accurate  and  lead  to  kit  solutions  not  as  good  as 
those  obtained  by  the  Greedy  Algorithm. 
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The  storing  of  probabilities,  and  interpolation  of  others,  does  appear 

to  possibly  contribute  to  errors  in  a  kit  configuration.  The  Greedy  Algorithm 

and  the  D029-IBM  program  produces  X*4g  ■  10,  but  the  D029-CDC  program  produces 

X*  =  1.  The  source  of  this  error  should  be  understood  and  eliminated  from 
249 

the  current  D029-CDC  program. 

Comparison  of  kit  results,  indicates  that  the  kits  obtained  by  the  Greedy 
Algorithm  or  the  D029  program  are  similar.  Approximately  75%  of  the  kit  items 
are  present  in  the  same  amounts  and  the  remainings  deviate  only  by  1  or  2. 

This  is  confirmation  of  the  recognition  that  the  current  D029  model  is  accept¬ 
able.  The  Greedy  Algorithm  appears  to  offer  the  potential  of  substantial  time 
savings,  and  a  slightly  better  kit  configuration. 
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APPENDIX  A 


A.l  Computer  Code  for  the  Primal  Search  Algorithm 

A. 2  Computer  Code  for  the  Modified  Algorithm  of  Brooks,  et  al. 

A. 3  Computer  Code  for  the  Saddle  Point  Search  Algorithm 


A  List  of  the  Variables  in  the  Computer  Program 


ALFA 

basla 

BETA 

BUTCE 

C(100) 

CBETA 

CDOTD 

CDOTX 

CKARE 

COST 

DENOM 

DIF 

DIREC 

ENORS 

EPS 

IOX 

ITETA 

IX( 100) 
IY(100) 
KODE 

N 

NEWP 


Stepsize 

The  Subroutine  that  finds  an  initial  point 
Product  of  the  parameters  p  (user  determined)  and 
Budget 

Cost  of  item  t 
The  parameter  p^ 

The  dot  product  of  the  cost  and  direction  vectors 

The  dot  product  of  the  cost  vector  and  the  current 
solution,  also  the  cost  of  the  kit 

The  dot  product  the  cost  vector  with  itself 

Cost  of  the  current  solution 

The  dot  product  of  the  Cost  vector  the  vector  of  first 
differences  at  a  solution 

The  first  difference 

The  direction  vector 

E(NORS/X) 

Parameter  for  testing  the  sufficiently  closeness 
Integer  Optimal  Solution 

Indicator  of  the  search  direction.  If  ITETA  *  1,  the 
search  is  in  the  positive  direction,  if  ITETA  =  -1, 
it  is  in  the  negative  direction. 

Current  integer  solution 

Initial  solution 

Indicator  of  the  result  of  the  search,  "l"  means  success 
"0"  means  failure 

Total  number  of  items  in  the  problem 
New  integer  point 

Number  of  iterations  already  performed 


NITER 


NNORS  Number  of  NORS  function  evaluations 

NSTEP  (NITER)  Number  of  steps  within  iteration  NITER 
M  Upper  bound 

PAR  Parameter  p  (user  determined) 

PDFOP  Function  routine  which  calculates  the  individual 

poisson  terms 

POISON  (40,100)  The  matrix  which  contains  the  cumulative  Poisson  sums 

SIRALA  (X,IR)  The  subroutine  which  ranks  the  array  X  in  descending 
order.  IR(k)  gives  the  rank  of  X(k). 

TABLE  (POISON)  The  subroutine  which  fills  in  the  matrix  POISON 


ucos 

Cost  of  NEWP 

UNORS 

E(NORS/NEWP) 

VNORS 

Function  Routine  to 

calculate  E(N0RS/X)  funtion 

XLAG 

Lagrange  parameter 

(user  determined) 

XLAM(I) 

Poisson  demand  rate 

of  item  I. 

A. 1.3 


FORTRAN  Listing  of  the  Computer  Program 


1 

2 

3 


4 

5 

6 

7 

8 
9 


10 


11 

12 

13 


14 


15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 


1 JCB  TIME-05 

c  *»*«**»**»»»«**«»****«» ************ 

c  •  * 

c  *  * 

C  *  A  PRIMAL  SEARCH  ALGORITHM  * 


Cm  * 

C  *  * 

c 

C  MAIN  PROGRAM:  TO  REAO  IN  THE  DATA 
C  TO  WRITE  OUT  THE  SOLUTIONS 

C 

DIHENS  ION  IXI 1001,10X1 1001 . 1 Y( 1001 .OE LTAI 100) 

CtJHKON  N,  M/ BLOCK  1/ Cl  100  l/BLOCK  2/XLAHI  ICO)  /BL0CK3/BUTCE  ,  PAR/BLGCK4. 
•  NI  TER »NSTEPI  5001/ BL 0CK5 /POISON (40  ,103  )/COUNT /  NNORS 
EPS- .00005 

C  INPUT  THE  NECESSARY  INFORMATION 
C 

NNORS-0 

REACI5.90)  NOPR 
RE AOI 5* 1 00)  N.M'BUTCE 
REACI5.  110)  (C<  I),  I*  If  N  ) 

REACI5.110)  (XLAMII  ),I«t.N) 

REAOI  5,120)  PAR, XLAG 
C 

C  CALCULATE  TEE  PROBABILITIES 
C 

CALL  TABLE  IPO  I  SON) 

C 

C  F I  NO  AN  INITIAL  SOLUTION 
CALL  8ASLAIXLAC.  IY) 

00  10  I-l.N 
10  IXII)-IY(I) 


C 

c 

c 

c 

c 

c 


START  THE  0PTIMI2ATI0N 

CALL  OPT IMIIX,  IOX, ENORStCOST) 


PRINT  THE  RESULTS 


WRITE  1 6 ,400) 

WRITEI6,  405) 

WRIT  El 6t205) 

WR  ITEC6, 200) 

DO  20  I-l.N 
WR  ITS  16, 210) 

WR  I TE(  6, 201) 

WRITE(6,205) 

WR  I TE ( 6 ,220) 

WR  IT  El  6,  250) 

FORMAT  I  14) 

FORM  AT  1214, F10.0) 

FO RM AT 110F6.0) 

FORMAT (2  FI  0.0) 

FORMAT!  ///,38X,* INI TIAL* ,5X,*0PTI MAL*  ,/,5X,* ITEM* ,3X, 'CCST/UN IT* ,6 
IX,  •  DEMAND*, 4X, ‘SOLUTION*,  4X,  'SOLUTION ',//) 

201  F0RMAT(///,5X,1  OHEf  NORS/X  I -,  FI  0 .5  ,5X,9KC0ST  (X  )-i,  F10.2) 

205  FORMAT!//,  5X,  ••  OF  NOR  S  EVALUATIONS**  ,14) 

209  FORMAT!///, SXfllH#  OF  ITEMS-, I4.3X, 8H BUDGET -t,F10.2) 

210  FORMAT! 5 X,*A* ,  13,2 I2X.F10.2 ) ,2 I8X, 14  I ) 


20 


90 

100 

110 

120 

200 


NOPR 

N.BUTCE 


1  ,C! I) ,XLAH! 11,17 (I) ,IOX (I) 

E  NOR  S, COST 

NNORS 

NITER,  INSTEP!!), 1-1, NITER) 
PAR,  XLAG 
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34  229  FORMAT  1//.5X,  •*  OF  ITERATIONS-  M  A.  /,  5X,  •  •  OF  STEPS  I  N  EACH  I  TERA 

•IONS' , !5X,2014)I 

35  259  FORMAT!  //• 5X,  'LAMBDA" ' «F 10. 2*  5X» 1  LAGRANGE  MULTI  PLI ER-' , FlO.B./t 

36  A 09  FORMAT tIHl,//,10X.» A  PRIMAL  SEARCH  HE7H90') 

37  405  FORMAT!///, 5X, 'PROBLEM  NUMBER' *  14  ) 

38  STOP 

39  ENO 


AO 


41 

A2 


A3 

AA 

A5 

A6 


AT 

48 

49 

50 


51 


52 

53 

54 

55 

56 

57 

58 


59 

60 


61 


I 


62 

63 

64 

65 

66 

67 

68 


i 


69 

70 

71 

72 

73 

74 


SUBROUTINE  OPT 1H1 IX, IOX ,EN0RS , COST! 

C  •  * 

C  •  OPTIMIZATION  ROUTINE  * 

C  •  • 

DIMENS  ION  IX(IOO),  I  OX!  lOJI.OIRECf  100)  'DELTA!  100)  .NEW  11  00)  , 
1XNEWI100) , NEIGH! 190 ) , I Y ( 100 ) < I  RANK (100), XOEL (1)9 ) 

COMMON  N ,M/B LDCK1/C ( 100 ) /BL0CK2/XLAM! 1091 /BL0CK3/BUTCE,  PAR/ BLOCK* 
•NITER, NST EP(500) 

C 

niter-i 

ENORS-VNORS(IX) 

COST-PROD! IX) 

CKARE-OOT  PR( C ) 

C 

C  EVALUATE  THE  FIRST  DIFFERENCES 
200  CALL  TJREV(  IX,  DELTA) 

NS TEP (NITER  1-0 
OENCM-OOTPR(OELTA) 

CSET  A-DENCM/CK  ARE 
C 

C  EVALUATE  THE  DIRECTION  PARAMETER  PK 
BET  A-PAR*C6£T.A 
C 

C  DETERMINE  THE  SEARCH  DIRECTION 
DO  10  I- 1,  N 

01 REC(I ) -OELTA ( I) -B  ETA*CI I ) 

MS  A-  ABS  (  CIR  EC!  1)1 
1FII.EQ.1 )  GO  TO  15 
IFIWSA.LE.XS!)  GO  TO  10 
15  XSI-WSA 

10  CONTINUE 
C 

C  NORMALIZE  THE  DIRECTION 
DC  11  I-i.N 

11  DJREC(  I  )-DIREC!  D/XSI 

C  CALCULATE  TFE  STEP  SIZE  ALPHA 
CDOTD'DOTPRI DIREC) 

CDOTX-PROO! IX) 

XN-SORT ( FLOAT ( N) ) 

AMAX-I  8UTCE-C0CTX)/ !COOTO*XN) 

IFIAMAX.LT. 1.0)  AMAX-1.9 
KMAX«IFIX(AMAX*.5) 

AMAX-FLOAT(MMAX) 

AL  F  A*  AM  AX 

c 

C  OBTAIN  A  NEW  POINT 
109  DO  20  1-1, N 

HN-FLOAT ( I X(  I ) 1+ALF A*D 1REC ! 1)4.5 

IFIWN.GE.l.)  GO  TO  20 

WN-1.9 

20  NEVPII) -IFIXIWN) 

NSTEPIN ITER )-NSTEPI  NITER) ♦! 


A.l.S 


C  CALCULATE  ITS  COST  $  EtNORS/X) 

75  UCOS-PRODINE  WP I 

.76  UNORS- VNORS IN  EWP I 

C  If  INFEASIBLE  .  PERFORM  A  SEARCH  IN  THE  NEGATIVE  DIRECTION 

77  1FILICOS.GT  .BUTCE  I  GO  TO  30 

C  If  A  BETTER  FEASIBLE  POINT.  REPEAT  THE  STEPS 

78  IF IUNORS.LT. ENORSI  GO  TO  500 

79  ITETA-1 

BO  GO  TO  40 

61  30  IT  ETA— 1 

62  40  CALL  SEARCH  ITETA.UCOS, UNORS, ENORS, NEKP.KOOE) 

63  IFIKOOE.EO.l  )  GO  TO  500 

84  IFIALFA .EO.l.l  GO  TO  900 

85  ALFA-ALfA-l. 

86  GO  TO  100 

67  900  UCOf.'COST 

89  U.NORS-ENORS 

89  CALL  SEARCH! 1 .UCOS . UNORS. ENORS • IX .KOOEt 

90  IFIXOOE)  600,1000,600 

91  500  00  45  1-1. N 

92  45  lXID-NEWPlII 

93  609  COST-UCOS 

94  ENORS-UNORS 

95  NITER-N ITER*1 

95  GO  TO  290 

57  1000  00  50  1-1, h 

98  50  IOXIII-IXII) 

99  RETURN 

100  END 

101  SUBROUTINE  SEARCHI I TETA ,UCO S. UNORS, ENORS, NEWP.KOOE » 

C 

C  THIS  SUBROUTINE  PERFORMS  SEARCHES  EITHER  IN  POSITIVE  OR  NEGATIVE  01 
C 

102  DIMENSION  NEWP 1100 1, NEIGH!  190) , IRANKI 100> , DELTA  1 1 00) , IN  SET! 100) ,D 
•  LI  ICO) 

103  COMMON  N,M/BL0CX1/C(  190  l/BLOC  K2/XL  AMI  100)  /BLOCK  3/BUTCE  ,  PAR/BLOC  K4 
•NITER, NSTEPI500I 

C 

104  NST  EPIN IT  ER  I-NSTEPI NITER) 41 

105  IF  II  TETA.  EO.l )  GO  TO  9 

C 

C  CALCULATE  THE  FIRST  DIFFERENCES  FOR  NEGATIVE  DIRECTION 

106  CALL  TUR NO I NEWP, DELTA) 

107  GO  TO  11 

C  CALCULATE  TFE  FIRST  DIFFERENCES 

108  9  CALL  TUREVINEWP, DELTA) 

109  11  DO  10  I-l.N 

110  CELT  All) -FLOAT  I  IT ET A  )*DEL T A 1 1) 

111  10  DEL  1 1 l-DELTAl I) 

C 

C  RANK  THE  FIRST  DIFFERENCES 

112  CALL  SIRALAIOEL, IRANKI 

113  DO  20  I-l.N 

114  00  15  J-l  ,  N 

115  JV-J 

116  IFI IRANKIJ ).EQ  .1)  GO  TO  17 

117  15  CONTINUE 

lie  17  IVSET1II-JV 

119  20  CONTINUE 
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c 

C  CHECK  WHET  HER  THERE  IS  A  BETTER  FEASIBLE  POINT 
C 

120  03  25  l«I ,N 

121  JV  AL- IMS  ETUI 

122  COS J-UCOS  ♦  FLOAT! ITETA1-CI JVAL » 

123  IFICOSJ.GT.BUTCEI  SO  TO  25 

124  CO  35  JJ-l.N 

125  35  HE ICHIJJ) -NEWPIJJI 

126  1FINEWPI  JVAL) .EQ.0.AN3.ITETA.EQ.-1I  GO  TO  25 

127  NEICHIJVALI-NEWPIJVALIUTETA 

128  XMORS-UNORS-DELTAIJVALI 

129  IF (XNORS -ENORS  )  100,120.110 

130  25  CONTINJE 

131  GO  TO  120 

132  110  00  43  I-l.N 

133  40  NEWPin-NElCHI  I) 

134  UVORS-XNORS 

135  UCOS-COSJ 

C  IF  NO  BETTER  FEASI8LE  POINT *  FAILURE  <  K03E-0 

136  120  KODE-O 

137  RETURN 
C 

C  IF  A  BETTER  FEASIBLE  POINT,  SUCCESS  S  KODE-1 
135  103  KODE-1 

139  03  70  I-l.N 

140  70  NEKP  HI  -NEIGH!  1) 

141  UNOR$=XNORS 

142  UCOS-COSJ 

143  RETURN 

144  BIO 

145  FUNCTION  PROOIIXI 
C 

C  TO  CALCULATE  THE  DDT  PRODUCT  C.IX 
C 


146 

DIMENSION  IXIIOOI 

147 

l 

COMMON  N, M/BL0CK1/C (100 1 

148 

SUM- 0.0 

149 

00  10  I-l.N 

150 

IXI-lXIt  ) 

151 

10 

SUM-  SUM  ♦  FL  OA  T l  IXI  !*CII» 

152 

PROD-SU  M 

153 

RETURN 

154 

ENO 

155 

FUNCTION  DOTPRIX) 

C 

c 

TO 

1  CALCULATE  THE  DOT  PRODUCT 

c 

156 

DIMENSION  XI100) 

157 

COMMON  N .M/BL0CK1/: I  1001 

158 

SUM-0.0 

159 

00  10  J-l.N 

160 

SUM-SUM *CC J)*X(JI 

161 

10 

CO NT INJE 

162 

dotpr-sun 

163 

RETURN 

164 

ENO 

1 

i 


165  FUNCTION  VNORSIIX) 

C 

C  TO  EVALUATE  THE  EINORS/Xl 
C 


166 

167 

168 

C 

169 

173 

171 

172 

173 

174 

175 

176 

177 

178 

179 

183 

181 

24 

182 

183 

25 

184 

185 

186 

26 

187 

20 

188 

189 

193 

151 

43 

192 

193 

194 

195 

195 

187 

198 

33 

199 

203 

201 

202 

DIMENSION  IX  1 100 1, P Ml  100.  20}. PRO! 201 
COMMON  N.M/BL0CK5 /POISON (40 .1 00 )/ COUNT/NNORS 
DOUBLE  PRECISION  SUM.PRO, TOPL.PM.XW 

EPS-. 00001 
6PS0-. 00000001 
NNORS-NNORS+1 
Ml-M+1 
SUM-0.0+00 
00  20  JJ-1.M1 
PRO! JJI-1 .0+00 
t»  26  II-l.N 
IXJ-IXI  in+JJ-1 
IF  I  IXJ.GT.O)  GO  TO  24 
XX-POFOPI 1XJ, 111 
GO  TO  25 

1FI  IXJ.GT.40)  IXJ-40 

xx -poison i  ixj,  in 

IFIXX.LT. EPSO)  XX-0.0 
IF (PRO! JJ  I  .LT.EPSQ}  GO  TO  20 
PRO(JJ) -PRCC JJ J*XX 
CONTINUE 

SUM- SUM +  1 .-PRO  I JJ I- 
TOPL-O.O+OO 
DO  30  I-l.N 
KO 

MXX- 1  XI I  )  +  M+2  +  K 
WX-POFOPIMXX. I | 

WW- FLOAT (K  +  l }*WX 

IFtHX.LT. EPS)  GO  TO  30 

TOPL-TOPL+WH 

K-K+l 

GO  TO  40 

CONTINUE 

VNORS-SUM+TOPL 

RETURN 

BIO 

SUBROUTINE  TUREVt IX. DELTA) 


203 

204 

205 

206 

207 

208 
209 


C 

C  TO  EVALUATE  THE  FIRST  DIFFERENCES 
C 

DIMENSION  IXI 100) .DELTA(IOO) 
COMMON  N.M 
C 

DO  10  K-l.N 
DELTAK  l-DIFI  IX, K  ) 

13  CONTINJE 
RETURN 
END 


210  SUBROUTINE  TURNOI NEHP, DELTA ) 

C 

C  TO  EVALUATE  THE  FIRST  DIFFERENCES  IN  NEGATIVE  DIRECTION 
C 

211  DIMENSION  NE WP 1 100) , NE I GH I  100) .DEL TAI 100 ) 
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COMMON 'N  ,N 


00  10  I-l.N 
DO  20  J-l.N 
1FIJ.EQ.I)  CO  TO  25 
NE  ICM( J  >*KEWPIJ) 

GO  TO  20 

25  NEIGH  <J  >*N  EWP  ( J  J—  X 

lEINElGHIJI.LT  .0)  NEIGHIJJ-0 
20  CONTINUE 

DELTA! I )•  D 1F!N  E IGH,  I ) 

10  CONTINUE 
RETURN 
END 

FUNCTION  DIF! IX, K) 

TO  COMPUTE  TEE  FIRST  DIFFERENCE  OF  ITEM  K 

DIMENSION  IX! 1001 

COMMON  N, M/BL0CK5 /POISON 140,103) 

OOUBLE  PRECISION  PRO,SUM,EPSQ 

EPS0*1 .D-08 
M1*M*1 
S*JM-C.0»07 
DO  20  Jl*l .Ml 
PROM. 0*00 
J-Jl-1 
00  30  1*1  ,N 
IXJ-  IX!  1  >*J 
IFII.EO.KI  GO  TO  25 
IFIIXJ.EQ.O)  GO  TO  23 
IF!  IXJ.GT.39i  IXJ-39 
PRO*PRO*POISON!IXJ,  II 
GO  TO  30 
23  !XJ1=IXJ 
GO  TO  26 

25  IXJ1MXJ+1 

26  IFIPRO.LT. EPSO)  GO  TO  20 
PRO  =  PRO*  PCEOP I  IXJ1,  I) 

30  CONTINUE 
20  SUM* SUM ♦PRO 
IXK*IX(K)*K+1 
IF! IXH.GT.4G)  IXM-40 
D1 F*SUM  *1  .-POISON! I XM,  X  I 
RE  TURN 
END 

SUBROUTINE  T  ABLE! PO  I  SON  ) 

TO  FILL  IN  TEE  ARRAY  FOR  CUMULATIVE  POISSON  SUMS 

DIMENSION  PO I SCN! 40  * 130  I 
COMMON  N »M/B LOC K2/X LAM | 100) 

OOUBLE  PRECISION  EPS,  TERM,  SUM  X 

EP  S"  1 .0-06 
03  19  J*  1,  N 
Tc  RM* 1 . 0*00 


I 


761 

5UNX- 1.0*00 

262 

DO  20  I *1 *40 

263 

TERM-TERM-XLAMI  JI/FLOATII  1 

266 

K-I 

265 

IFITERH.LT.EPS)  GO  TO  25 

266 

SUMX-SUHX*TERM 

267 

POISON  I I*J)=EXPt-XLAM(J  H-SUNX 

266 

IF  (POISONU.J).GT.l  .0)  POISON!  I, J )*1 .9 

269 

20 

CONT  INUE 

270 

25 

IF! IX.E0.40)  GO  TO  10 

271 

DO  30  ! *  1  X*40 

272 

30 

POISON! I, JJ-1.0 

273 

10 

CONT INUE 

274 

RETURN 

275 

END 

276  FUNCTION  PDFOP1K,!) 

C 

C  TO  CALCULATE  JNOIVIDJAL  TERNS  OF  THE  P3ISSON  DISTRIBUTION 


I 


C 

277 

COMMON  /BLOCK 2 /XL AMI  133) 

278 

DOUBLE  PRECISION  SUM.WW 

279 

t 

SUM-0.9  *00 

280 

IF (K.LT .2)  GO  TO  10 

281 

DO  15  J-2.K 

282 

SUM*  SUM* ALOG I  FLOAT!  J ) ) 

283 

15 

CONT  INJE 

284 

10 

WW--XLAMI I  )♦  K*  A LOG!  XLAMI I ) ) -SUM 

285 

IFIKW.LT.-15.01  GO  TO  20 

286 

PDFUP-OEXPIWW) 

287 

GO  TO  25 

288 

20 

PDFOP-O .0 

289 

25 

RETURN 

290 

END 

291 

SUBROUTINE  SIRALAIX*  IR  I 

c 

THIS  SUBROUTINE  RANKS  AN  ARRAY  OF  ELEMENTS 

c 

c 

IRIK)  IS  THE  RANK  OF  ITEM  K  IN  THIS  3R0ERINI 

292 

L 

COMMON  N,M 

293 

r 

DIMENSION  XI  100)*  IR  I  100),  1SEQI  100),  Y!  109) 

294 

00  9  I-l.N 

295 

9 

ISEOID-I 

256 

DO  10  I -1  ,  N 

297 

XM IN- XIII 

298 

IMINJ-1 

299 

IM  IN-  I  SEQ  1 1) 

300 

Nl-N-  I-l 

*01 

DO  20  J-l  . N1 

302 

IF ( XI J)  .GT.XHINI  GO  TO  21 

303 

GO  TO  23 

304 

21 

XHIN-XI J) 

335 

IM  IN- ISEOI Jl 

30S 

IMINJ-J 

307 

20 

CONTINUE 

308 

IR  I  IM  IN  1  ■  I 

309 

IFIN1.E0.1  )  GO  TO  10 
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311 

03  30  J«1,N2 

312 

IMJ.LT.  I MIN J)  GO  TO  33 

313 

Jl«J*l 

31* 

ISEQUI-lSEOtJl  1 

31$ 

YC  JI-XIJll 

316 

GO  TO  29 

317 

33 

Y( J)*X( J ) 

318 

29 

XIJJ-YI  J> 

319 

30 

CONTI  NJE 

320 

10 

CONTINUE 

321 

RETURN 

322 

END 

323 

SUBROUTINE  BASLAI XL AC' IX 1 

c 

C  THIS  SUBROUT  IN  FINOS  AN  INITIAL  POINT  WITH  A  PREDETERMINED 
C  LAGRANGE  PARAMETER 
C 


32* 

DIMENSION  IXI130) 

325 

COMMON  N .M/BLOCKl/C 1100) /BLOC KS/P0IS0NI43 ,1001 

326 

V 

-C 

■ 

o 

s 

327 

Iw-0 

328 

A=  XL AG*C (  IJ*PDFOP<  IW,I)-1. 

329 

IF  I A.GE.O.O)  GO  TO  10 

330 

15 

IW* IW*1 

331 

B>  AfPDFOPI  IW,  I) 

332 

IF  (A.LE.O.O.ANO.B.GT.0.0)  GO  TO  10 

333 

A-B 

33* 

GO  TO  15 

335 

10 

IXUIMM 

336 

RETURN 

337 

END 
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A  List  of  the  Variables  in  the  Computer  Program 


! 

1 

9 


BO 

C(I) 

COST 

ENORS 

EPS 

FC100.50) 

ICOUN 

IMAX 

IX(100) 

M 

N 

NOPR 

NSTEP  (ICOUN) 

OPTIM 

PDFOP 

RATIO  (100) 

RMAX 

TABLE 

XLAM  (100) 

XINVT 

VNORS 


Coefficients  bn>  b,  , b 
0  1  m 

Cost  of  item  I 

Cost  of  the  current  solution  X 
E(N0RS/X) 

Parameter  for  testing  sufficiently  closeness 
Matrix  for  the  cumulative  poisson  sums 
Counter  for  the  number  of  Iterations 
The  index  I  for  which  RATIO(I)  is  maximum 
Current  solution 
An  upper  bound 

The  total  number  of  items  under  consideration 
Problem  number 

Number  of  steps  in  iteration  ICOUN 
The  optimization  routine 

Function  routine  which  calculates  the  individual  poisson 
terms 

Ratio  of  return  per  dollar  investment 
Ratio  (IMAX) 

The  subroutine  to  fill  in  the  matrix  F(100,50) 

Poisson  demand  rates 
Budget 

Function  routine  to  calculate  E(N0RS/X) 
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The  FORTRAN  Listing  of  the  Computer  Program 


P 


< 

C  * 

c 

C  MODIFIED  SOLUTION  METHOD  OF  BROCKS  ET  At 
C 


DIMENSION  IX (103)  , 80(1  001 

COMMON  N.M/BLOCKl/C  I  100),  XI NVT/BL OC K2/XLAMI 100)  /B L0CK4/F C 1 00,50 ) / 
•L0CK5/1C0UN.NSTEP  (50) 

EPS*  .00  1 
EPSC-EPS**2 

INPUT  THE  CATA 

READ! 5,90)  NOPR 
RE  A0(5,  100  )  N,  M.XINVT 
RE  A0(  5,110)  ICII),iM,NI 
REACI5,  110)  (XLAMtl  I.I-l.N) 

HR  IT  El  6,  400) 

VR  I TEI  6, 405)  NOPR 
WRITEI6, LOS)  N.XMVT 
ML  -H*l 

INITIALIZE  TEE  PARAMETERS  B 
00  LO  J-l.M 

14  io  eou)-o.o 

15  B0IMD-1.0 
C 

C  CALCULATE  ThE  PROBABILITIES 
15  CALL  TA6LEIF) 

17  ICOUN-1 
C 

C  PERFORM  THE  OPTIMIZATION 
C 

18  1000  CALL  OPTIMIBO, IX, COST) 

C 

C  CHECK  IF  B-NEW  IS  SUFFICIENTLY  CLOSE  TO  B-OLO 
C 


19 

SUMSQ-0  .0 

20 

00  20  Jl-l.Ml 

21 

PRO-  l  .0 

22 

03  40  15-1  ,N 

23 

IXJ-IXI  I  5 )  +  J 1 

24 

IF  I  IXJ.GT  .50  )  GO  TO  40 

25 

PR0-PR0*F(15,IXJ) 

26 

40 

CONT  INUE 

27 

1FIPR0.LT .EPSO)  PRO* 0. 0 

28 

SUHSQ-SUMStHIIBOUl  1-PR0)**Z  ) 

29 

60  UD-PRO 

30 

20 

CONT INJ  E 

C 

C  IF  SO,  TERMINATE  AND  PRINT  THE  RESULTS 
C 

31  IFtSUMSC.LT.EPS)  GO  TO  500 
•.  C  OTHERWISE  REPEAT  THE  STEPS 

32  ICOUN-ICOUN*! 


A. 2. 3 


CO  70  1009 
SOO  WRITE! 6.  240) 

03  3C  1-1,1* 

30  WR  IT E!6,  250)  ],C(  D.XUMI  I  ),I  X(  1  > 

ENORS -7N0RSI  IX) 

WR I TE( 6, 201 )  E NCR S .COST 
WK  IT EI6,  330  )  (  BO!  J  I ,  J-  1,  M  l) 

MUTE  (6, 301)  ICOUNf  I  NS  T  EP  (  I  ) ,  1-1 ,  ICOJN I 
93  FORMAT!  14  I 
100  FORMAT  !  2  14, F  19.0) 

105  FORMAT!//, 5X,*  N-*  ,14, 5X,*  BUDGET-**  ,  H  0  .2  I 
110  FORMAT! 10FB.01 

231  FORMAT  !// ,  5X,  1  OHE  iNORS/X)-  ,  F10 .5,  2X  .9HC0ST!  X)  «S,F  1  0.2) 

240  F3RMATI ///,  38X,  *  OPTI  HAL*  ,/,5X,*  ITEM*,  3X,*  COST/ JNIT'.SX,  'DEMAND*,  4> 
1, ‘SOLUTION*,// ) 

250  FORMAT ! SX « '  f ' ,  I3,2!2X,F13.2),8X,I4) 

300  FORMAT!//,  5X,  *  B  :*  ,  1  OF  1  0.3) 

331  FORMAT  !///,  5X,  16HI  OF  ITERATIONS-  ,  14,  //»!  5X,  •  f  Or  STEPS  IN  EACH  I* 
•ERATIONJ*  ,2014)) 

400  FORMAT  (1  Hi,//,  SX,  •  THE  MODIFIED  SOLUTION  METHOD  OF  BROOKS  ET  AL'*/. 
♦  I 

405  FORMAT! ///,5X,‘ PROBLEM  NUMBER', 14  I 
STOP 
END 

SUBROUTINE  OPT  1MI BO, IX, COST) 

OPTIMIZATION  ROUTINE  PERFORMS  THE  MARGINAL  ANALYSIS 
DIMENS  ION  1X11001,  BO! 130),  RATIO!  1031, MR!  100) 

COMMON  N  ,M/BL0CK1/C  1 103  I,  X  INVT/BLOCKV  /F!  133 , 53  )/  BLOCK 5/  ICOUN.NSTEl 
*150) 

Ml  -M*l 
COST-O.O 
09  13  I-l.N 
I X  I  1 1  -3 


CALCULATE  THE  RATIO  OF  RETURN  PER  EACH  DOLLAR  INVESTED  FOR  EACH  I  TE 
SUM-0.0 
00  20  J-  1, Ml 
JJ-J+1 

WW-FU,  JJl/FI  I,  J) 

20  SUM«SUM*eO(J  )*AL3GI  WW) 

10  RAT  10! 1 1 -SUH/C ! 1 1 

NSTEP I ICOUNl-l 
200  NOR-3 

NSTEP1IC0UNI  -NSTEP!  ICOUNl  +  l 
00  21  I-l.N 

21  MRID— 1 
103  IMR-0 


FIND  THE  MAX  1MJM  RATIO  *  ADO  ONE  UNIT  3F  THIS  ITEM  TO  THE  KIT 


03  20  I-l.N 

IF(HRil).EO.l)  GO  TO  25 
IMR-I MR* 1 

IF! IMR.EO.l)  GO  TO  26 
IFIRAT  101  1  l.LE.RMAX)  G3  T9  25 
26  HMAX-RAT 10! I ) 


79  !m  AX  ■  I 

60  25  CONTINUE 

Si  MRIlMAXI-l 

c 

C  CALCULATE  THE  COST  OF  THE  NEW  KIT 

62  XCOST-COST+CI  I  HA  XI 
C 

C  IF  COST  EXCEEDS  THE  BUDGET.  FIND  THE  IT  EH  WHICH  W  ILL  VIELD  MAXIMUM 
C  WITHOUT  EXCEEDING  THE  BUDGET*  THEN  BE  TURN  THE  MAIN  PAOG. 

C 

S3  IF(XCOST.GT.XINVT)  GO  TO  30 

64  C9ST-XCOST 

85  IX  UMAX  J«  I XI  1MAXI +1 

66  SUM-0.0 

87  00  35  J-l.Mi 

63  IX  J-IXI  INAXI+J 

89  IXJ1-1XJ*1 

90  IF(  IXJ.GE.SO>  GO  TO  35 

91  WK-FI  IHAX,  IXJ1  I/FI  IMAX,  IXJ1 

92  SUM-SUM* BC  l  J>*  ALOGIWWI 

93  35  CONTINUE 

C  OTHERWISE  U  PCAT 6>  THE  RATIDIIMAX1,  REPEAT  THE  STEFS 

94  RATIOI I MAX I" SUM/C (  I  MAX) 

95  GO  TO  200 

96  30  N0R«N0R*1 

97  IF  (N0K.LTt.NI  GO  TO  100 

93  RETJRN 

99  END 


100 


101 

102 

103 

104 

105 
10S 

107 

108 
109 

no 

111 

112 

113 

114 

115 

116 

117 

118 

119 

120 
121 
122 


SUBROUTINE  TABLEIF ) 

C 

C  CALCULATING  CUMULATIVE  POISSON  SUNS 
C 

DIMENSION  FI  100.  50) 

CUMMON  N .M/BLOClU  /CI1001.XI NVT/ BLOCK? /X L AMI  133  ) 
DOUBLE  PRECISION  EP  S. TERM,  SUM 
C 

EPS-1. D-06 
DO  10  I-l.N 
FI  1,1  I*  POFCPIO, 1 1 
TERM-1. 0*00 
SUM-1 .0*00 
DO  30  J-2,50 

TERM»TERH*XLAM< ll/FlOATfJ-lI 
JX-J 

IF  I  TERM. LT . EPS  I  GO  TO  25 
SUM- SUM* TERM 

FI  1, Jl-EXPI-XLAMt I) ) -SUM 
1FIFII, Jl.GT.l.OI  FII.J1-1.0 
30  CONTINUE 

25  IFIJX.EQ.50i  GC  TO  10 
03  35  J-JX.50 
35  FII.JI-1.0 
10  CONTINUE 
RETURN 
END 


123  FUNCTION  POFOP I K, I ) 

C 

C  CALCULATING  THE  INDIVIDUAL  POISSON  TERMS 
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c 


12* 

COMMON  /BLCCK2/XLAM 11001 

125 

f* 

DOUBLE  PRECISION  SUM,  WM 

126 

L 

SUM-O.O^OO 

127 

IFIK.LT.2)  CO  TO  10 

128 

00  15  J-2 ,K 

129 

SUM- SUH+ALOGI FLOAT!  1)1 

139 

15 

CO  NT  1NJE 

131 

10 

WW--XLAMI  1>4K«AL0S(  XLAMII  M-SUM 

132 

IF  IWW.LT. -15.0)  GO  TO  20 

133 

POFOP-OEXPIWW) 

13* 

GO  TO  25 

135 

20 

POFOP-9 ,0 

136 

25 

RE  TURN 

137 

END 

138 

f 

FUNCTION  VNORS  1 IX  ) 

c 

f 

CALCULATING  THE  EINORS/XI  FUNCTION 

139 

€ 

DIMENSION  1X11001 .PR0I20I 

1*0 

COMMON  N,M/BLOCK*/F I  100,50) 

1*1 

OOUBLE  PRECIS  ICN  SUM, PRO,  TOPL 

1*2 

t 

EPS- .00001 

1*3 

Ml-M+1 

1** 

SUM-0.0400 

1*5 

DO  10  J- 1»N1 

1*6 

PROIJI-1.C400 

1*7 

00  20  I -1 , K 

1*8 

IXJ-  1X1  I  )4J 

1*3 

20 

PRO! J)*  PRGIJ )*F(I»IXJ) 

150 

10 

SUM-S  IM41  .-PROI  Jl 

151 

TOPL- 0. 0400 

152 

00  30  J -1 , N 

153 

K-0 

15* 

*9 

MXK-  1X1  I ) «M4j4X 

155 

WX»POrOP<FXK,I 1 

156 

KW-FLOATC K41 |*WX 

157 

IF  (WX.LT  .EPS  )  GO  TO  30 

158 

TQPL-T0PL4WW 

159 

K-K41 

160 

GO  TO  *0 

161 

30 

CONTINUE 

162 

VNORS- SUM4T0PL 

163 

RETURN 

16* 

END 

SENTRY 
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A  List  of  the  Variables  in  the  Computer  Program 


A(50) 

AXM 

AXP 

BASLA 

BUTCE 

C(100) 

COST 

CX(50) 

CXM 

CXP 

DIF 

DOTPR 

DUALPR 

ENORS 

INPUT 

ISTAR 

IX(100) 

1X0(100) 

1X1(100) 

1X2(100) 

1X0(100) 

IY(100) 

KPAR 

M 


The  cost  of  the  current  solution  for  subroutine  DUALPR 
Minimum  of  A(I)  which  exceed  BUTCE;  A(IXM) 

Maximum  of  A(l)  which  are  less  than  BUTCE;  A(IXP) 

Subroutine  which  finds  an  initial  point 

Budget 

Cost  of  the  item 

Total  cost  of  the  current  solution 
E(N0RS/X)  (used  only  in  DUALPR) 

CX(IXM) 

CX(IXP) 

Function  routine  for  calculating  first  differences 
The  dot  product  of  two  vectors 

Subroutine  which  solves  the  dual  problem  to  obtain  the 
Lagrange  multiplier 

E(N0RS/X) 

The  subroutine  which  reads  in  the  necessary  data  and 
initializes  the  arrays 

Starting  solution  X° 

Current  solution 

Solution  obtained  in  the  Subroutine  BASLA 
Solution  obtained  in  the  Subroutine  SEARCH 
Solution  obtained  in  the  Subroutine  TARAMA 
Final  and/or  optimal  solution 
Initial  solution 

Parameter  for  finding  an  initial  point  in  BASLA 
An  upper  bound 


« 


N 

Total  number  of  items  under  consideration 

NARA 

Number  of  times  Subroutine  ARAMA  called 

NBAS 

Number  of  time  Subroutine  BASLA  called 

NITER 

Number  of  iterations 

NOPG 

Number  of  points  generated  by  the  Saddle  Point  Search 
Algorithm 

NOPR 

Problem  number 

NPAl 

Number  of  times  Part  1  called 

NPA2 

Number  of  times  Part  2  called 

NPA3 

Number  of  times  Part  3  called 

NSEA 

Number  of  times  SEARCH  called 

NVNOR 

Number  of  E(NORS/X)  function  evaluations 

OPTIM 

Optimization  Routine 

PARTI 

The  subroutine  which  performs  a  test  by  calculating  an 
upper  bound  to  the  linear  program  in  determining 
whether  a  component  can  be  increased  by  one 

PART2 

The  subroutine  which  calculates  an  upper  bound  for  the 
sum  of  first  differences 

PART3 

The  subroutine  which  performs  the  same  test  as  in  PARTI 
but  sets  up  the  linear  program  explicitly 

PDFOP 

Function  routine  for  calculating  individual  poisson 
terms 

PLAG 

Current  Lagrange  multiplier 

POISON 

Matrix  for  the  cumulative  poisson  terms 

POLD 

Previous  Lagrange  multiplier 

SEARCH 

The  subroutine  that  makes  a  one  dimensional  search  for 
each  component  separately 

SIRALA  (X,IR) 

The  subroutine  that  ranks  the  array  X  in  descending 
order 

TABLE 

The  subroutine  which  fills  in  the  matrix  POISON 
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TARAMA 

VNORS 

VXO 

VX1 

VX2 

XLAM(IOO) 

XLPRO 


The  subroutine  which  makes  an  attempt  for  improvement 
heuristlcally 

The  function  routine  to  calculate  E(NORS/X)  +  X  CX 
E(N0RS/IX0) 

E(N0RS/IX1) 

E(N0RS/1X2) 

Poisson  demand  rates 

The  Linear  Programming  Routine  which  solves  problems 
of  the  type  max  CX  subject  to  AX  ^  b  X  >_  0. 
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FORTRAN  Listing  of  the  Computer  Program 


f JOB  TIME-OS 


c  •  * 

c  •  * 

C  *  A  SADDLE  POINT  SEARCH  ALGORITHM  * 
C  *  * 

c  •  • 

C  ft****.*.*****,**,. *•*„*********•**** 

c 


C  MAIN  PROGRAM 
C 

DIMENSION  IX  (1001.  I  Yl  1001 
INPUT  THE  DATA 

CALL  INPUT  (NOPR.KPAR (PLAG.  1YI 
PERFORM  THE  OPTIMIZATION 

CALL  OPTIMIPLAGiKPAR.IY.  E  NOR  St  IX, COST  I 
PRINT  THE  RESULTS 

CALL  QTPUT(PLAG,NOPR.I Y.ENORS. I X, CO ST  I 

STOP 
ENO 

SUBROUT INE  INPUT! NOPR, KPAR, PL  AG. IYI 
TO  READ  IN  TFE  DATA 
DIMENSION  1Y(1P0) 

COMMON  N, M /BLOCK  1/C C  103  I/BLOC K2/XLAP.I  1001 /9LOCK3/BUTCE 
COMMON  /COUNT/ NITER,  NBAS  (SO  I.  NS  EAI 50 1  .NT  AR I  50),  NARAISO)  .NPA1I50I 
*NP A2( SO  I , NPA3I SOI .NVNOR/TAB/POI SONI 50 .IDS) 

RE  AD IS, 90  )  NOPR 
READ!  9,1001  N..M.3UTC: 

RE  ADIS,  113)  ICC  1).  1*1,  N) 

RE  AOI 9,11  31  ( XLAM (I ) , 1* 1 , N) 

REAOI5, 1501  KPAR, ISTAR, PLAG 

INITIALIZE  THE  COUNTERS 
NITER-1 
NVN9R-3 
DO  19  1-1, SO 
NB  AS  (  11-0 
MS  EA(  I  )*0 
NT  AR( II >9 
NARAIII-9 
NPA1  m«o 
NPA2(II«9 
NPA3I I)*0 
19  CONTINUE 
00  10  I-l.N 
13  IY  (!)  *  1ST  Ak 

CALCULATE  THE  CUMULATIVE  POISSON  SUMS 


CALL  TABL E (PO 1  SON) 


90  FORMAT  (14) 

100  FORNATI2  14, FI  3. Ot 
110  FORMAT! 10FB.0) 

150  FORMAT  12  14,  F13.0) 

RE  TURN 
END 

SUBROUTINE  OT  PUT  I  FLAG,  NOPR,  I Y,  ENOR  S,  l  X,C3ST) 

WRITING  OUT  THE  RESULTS 
DIMENSION  1 Y II  JO) (IX (1001 

COMMON  NfM/BLOCKl/CI  100) /BLOCK 2/ XL AM( 103 ) /BL0CK3 / BUTC E 

COMMON  /  COUNT/ NIT  ER,NBAS(53I,NSEA  (50)  tNTARI  53),NARA(50> ,NPAl(50)  , 

•  NPA2I53) , NPA3 ( 50) , NVNOR/T  AB/POISONIS) tl33 I 

WRIT  EI6.3G3) 

Wft  I  TE (6,435)  NCPR 
WRITE  (6,  235)  N.BUTCE 
WP.ITE(6,22  J) 

03  20  I  *  1  ,N 

20  WR  IT E( 6,  225)  I ,C(  I )  ,  XL AHI  I)  ,1  Yl  I)  , I  X(  11 
WW  I  TE(  6  *21  3)  ENORStCOST 
WAITE!  6, 23  3)  P LAS 
WRITE  (6,243)  NITER,  NVN3R 
WRI TE (6 ,245) 

D3  30  I- 1, NITER 

30  WR  IT  E(6,  25  )  I  I,  NBAS  ( I),  NSEAI I  )  ,NTAA(  I  ),NAR4(  I  J.NPAld  )  ,NPA2(!)  ,NP/ 
13(1) 

205  F3RMATI5X,  'NUMBER  3F  I TEXS«‘  ,  1 4,2X, * B UDGE T*»  •  ,F  13. 2) 

220  FORMAT  (///  .3  8X,  •  IN  1TIAL  ',  5X,  ‘OPTIMAL  •  ,/,  5X,  ‘ITEM*  ,  3X,  ‘COST/INI  T*  ,6 
1  X,  •  DEMAND'  ,4X,  •  SOLUTION'  ,4X ,‘  SOLUTI  ON*  ,//  ) 

225  FORM  ATI  5 X,  •••,  13,21 2X.F10.2), 2( BX.I  4) ) 

210  FORMAT! //,5X,‘ El  NOR S/X ) «• , F 10. 5, 5X, ‘COST (X)  =  »• ,  FI  0.2 ) 

230  FORMAT! /,5X, ‘LAGRANGE  MULTI PLI ER* ' , FI 0.6) 

240  FO.*  A  AT  (// ,  5X  ,  1  »  OF  ITERAT I3N S- • , I  4, // , 5X,  • •  OF  NOR  $  E  VA  LUATI  ONS«‘  , 
1 14  ) 

245  FORMAT!//, 5X, • $  OF  STEPS  IN  EACH  I  TE RATI 0N‘ ,//,5X ,* 1 TE RATION  f‘,4 

♦  X  ,  •  ST  EP1  SIEP2  STEP3  STEP  4  PUT  1  PART2  PA  RT3  '  ,/) 

25C  F3RMAT(12X,I4,2X,7(3X,I41) 

300  FORMAT  1 1H1, //,  5X, ‘THE  SA03LE  POINT  SEARCH  ALGORITHM',//) 

405  FORMAT!/// ,5X, ‘PROBLEM  NJMBER * , 14,/ I 
RETURN 
ENO 

SUBROUTINE  OPT IM( PL AG, KPAR, IYC.VX A,  IXO.COST) 

OPT  1 M I  Z  AT  ION  RO:JT  INE 

DIMENSION  1X3(  1031,1X1(1331  ,1X0(133),  I Y01103 ) , CX (50 1, At 50  ),  1X2(100 
*1 

COMMON  N, M/BLOC K3/ BUTCE/BL0CK4/NQPG/BL0CK 1/C ( 130  I 
COMMON  /COUNT/NITER  ,  NB  A  S<  53 )  ,  NSE  A  ( 50)  ,NTARC5D,  NARA  (5  3 )  .NPA1I53  ), 
*NP A2 (5) ) ,NPA3t  50 ),  N VNJR 


c 

I 


73 

74 

75 

76 


77 

78 


79 


80 

81 

82 

83 


84 

85 

86 
87 
85 
89 


90 

91 

92 

93 

94 


95 


96 

97 

98 

99 
100 
101 


C 

C  STEPIt  FIND  AN  INITIAL  POINT 
C 

500  CALL  BASLA(PLAG,KPAR,1  XOI 
If  (NITER. GT.l)  GO  TO  20 
00  25  1-1,6 
25  ivom-ixom 
C 

C  STEP2 :  FIND  A  STATIONARY  POINT 
C 

20  CALL  SEARCMPLAS,  1X3,  1X1) 

VX1  -VNOP.SC  PLAG  •  1X1 1 
C 

C  ST  EPS :  A  HEURISTIC  1MPR3V&ENT  TEST 
C 

CALL  TARAHA(PLAG, IX1,VX1,I X2, VX2I 
C 

C  COMPARE  THE  NEW  POINT  WITH  THE  INCUMBENT  POINT 
If ( VX2- VX 1 1  14,12,12 

14  no  13  I-  1,  N 
13  1X0(11-1X2(1) 

GO  TO  20 
C 

C  STEP4J  KE0UCT10N  OF  POTENTIAL  CANDIDATES  BY  CONTRADICTION 
C 

12  CALL  ARAM/.IPLAG,  IX2.VX2,  IXO.VXOI 
IF  (VXO-VXll  10,15,15 

13  00  11  I-  i,N 

11  1X3111-1X0(1)  * 

GO  TO  20 

15  N3PC-N3PG+1 
C 

C  READJUST  THE  C0EFF1CE INT  OF  THE  LINEAR  PROGRAM 
C 

A ( NOPG) -  DOT PR (  1X0) 

COST-A(NOPG) 

VXA-VHORSIZZ, IXO) 

CX (NOPG) -VXA 
POLO-PLAG 
C 

C  SOLVE  THE  OUAL  LP  TO  OBTAIN  THE  LAGRANGE  MULTIPLIER 
C 

CALL  OJ  ALPRI CX, A,  PL  AG) 

C 

C  IF  THE  NEW  MULTIPLIER  IS  SUFFICIENTLT  CLOSE  TO  OLD  ONE,  TERMINATE 
C  OTHERWISE  REPEAT  THE  ITERATIONS 
C 

FARK«ABS(PCLO-PLAGI 
IF(FARK.LT  .EPS  )  GO  TO  1000 
Nl  TER-N1 TERfl 
GO  TO  500 
1000  RETURN 
END 


) 


i 


H 


i 


102  SUBROUT  1NE  OUALPR(CX,A,  X) 

C 

C  SOLVE  THE  DUAL  TO  OBTAIN  A  LAGRANGE  MULTIPLIER 
C 

103  01  MENS  ION  CX(59),A(50> 

104  COMMON  /BlCCK3/BUT:E/BL0CK4/N0PG 
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oo  io  i>i.r.0PC 
IF  (ACM -8U1CEI  20.29.33 
20  IY«IV*1 

iriIY.EQ.1)  CO  TO  25 
IF  (AI1I-AXPI  19.19.25 
25  AX  P«  A  ( 1 1 
CXP-CXl I  ) 

GO  TO  10 
30  lXsIX+1 

IF  (  IX.EO.l  I  CO  TO  35 
IFIAIII-AXMI  35.19.19 
35  AX  M«  All) 

CXH«CX(II 

10  CONTINUE 

X*  l CXP-CX  M  )/  I AJW-AXP  ) 

RE  TURN 
END 

SUBROUTINE  BASLAI  XL  AC.  KPAR.  IX 1 

FINOS  AN  INITIAL  POINT 

DIMENSION  IX!  1901 
COMMON  N.M/BL3CK1/CI  1 J 3  I 

COMMON  /COUNT /NITER. NS AS (53  I .NS EA 1 50) • NT ARCS 9 1.  NARAI53I .NPA1I 53) 
•  NPA2I53)  .NPA3I  S9I.NVN0R 
NS  AS  (  N  IT  ER  >■  NBAS  INI  TER  )♦! 

DO  10  1*1 «N 
lvt  ■  0 

A=XLAG*CtI)  +  PDFOPI  IN,  I  1-1. 

IFIA.CE.9. 31  GO  TO  11 
15  lU*Ih»l 

B*  A*POFOP(  IN,  I  1 

IMA.LE.O.O.  AN0.8.CT.3.3)  CO  TO  11 
A*  6 

GO  TO  15 

11  I  «*'*  I W-KP  AR 

IF  (  IN  .LT  .3)  IW-0 
10  IX(I)»IW 
RETURN 
ENO 

SUBROUTINE  SEARCHIXLAG, IX , 1X01 

PERFORMS  ONE  DIMENSIONAL  SEARCH  FDR  EACH  COMPONENT  SEPARATELY 

DIMENSION  I XI  1301,1  X0( 103) .1 V ( 109) 

COMMON  N.M 

COMMON  /COUNT/ NITER ,N3AS(53).NSEA[>3 ) .NT  AR (5) I,  NAR A (50 ) ,  NPA1I 53 ) 
*NP  A2( 53 I.NPA3C  50I.N  VNQR 

NSEAINI  TERI  *NS  EA(  NITER  1  *1 
D3  10  I-l.N 

CALCULATE  THE  FIRST  DIFFERENCE 
X1*0IF(XLAG,  IX, II 


C  IF  IT  .IS  POSITIVE  ,  THE  POINT  XIII  IS  STATIONARY.  OTHERWISE 

C  FIND  XIII  SUCH  THU  ITS  FIRST  DIFFERENCE  IS  POSITIVE 
C 

150  IF!  XII  2C.2S.25 

151  20  DO  11  J*  1.  N 

152  IF  1 1 . CO. Jl  CO  TO  12 

153  IY  IJ  l»  IXI J  I 

ISA  GO  TO  11 

1SS  12  I  Y(  J)  *  I  XI  J  )+l 

ISA  11  COKT  IN'JE 

1ST  YlaDIFIXLAC.IY.il 

158  IFIY1I  39.36.  36 

159  30  00  13  J-l.N 

160  IF  1 1- Jl  13.16.13 

161  16  IX  I J  I  *  I X I J  1+1 

162  13  CONTPUE 

163  Xl-Yl 

166  CO  TO  20 

165  25  00  15  J«1.N 

166  IF  II. EQ. Jl  GO  TO  16 

167  IY  IJ  |a  IXI  J  I 

168  CO  TO  1 5 

169  16  IVIJIalXI Jl-l 

170  IFIIYIJI  .LT.OI  G3  TO  36 

171  15  CONTINUE 

172  Y?=OI FI  XLAC. IY.II 

173  IFIY2I  36,  36.  69 

176  60  00  60  J-l.N 

175  60  IXMI-IYIJI 

17*.  XI  >  Y2 

177  G3  TO  25 

178  36  03  65  J-l.N 

179  65  IX9(J)aJX|J) 

180  10  C3NT  IN'JE 

18!  RETURN 

182  END 

183  SUBROUTINE  TAR AHA  I  XLAG.IX.VX.I X0>  VXOI 
C 

C  OBTAINS  IMPROVEMENT  HEUR  1ST  I CALLY 
C 

186  01  HENS  ION  1X1100).  I  Yl  100).  RATIO!  1001 .  IR4NKI  1001  ,1X0(100) 

1 85  C3MM0N  H.H 

186  COMMON  /COUNT/NITER , NBAS! 501. NSEAI50)  .NTARIS3I, KARAI50) ,NPA1 (531, 
•  NPA2I53  ),  NPA3I  50),  SVN3R 

C 

187  NTARINITER  l=NTAR( NI TER) ♦ 1 

188  IT  ET  A=1 

1»9  100  00  10  I  * 1  ,N 

190  00  20  J«  1.  N 

191  IFIJ.EQ.I)  GO  TO  25 

192  IFIITETA. EC.-l .ANO.IXI JI.EQ.3)  GO  70  25 

193  IT  (J  )»  IX  IJ  )♦  IT  ETA 

196  CO  TO  20 

195  25  1Y(J)aIX(J) 

196  20  CONTINUE 

197  ALFB«OI  FIXLAC.IY.I) 

199  AL rAaOIFI XLAC,  IX,  I  I 

199  RAT I0(I)-ALF8/ ALFA 

200  10  CONTINUE 


RtNK  THE  RATIO  OF  FIRST  DIFFERENCES  OF  EACH  COMPONENT 

CALL  SIR  AC  AIR  AT  IOt  IRAN*  I 
DO  SO  I  -1  »  R 
00  35  J-  1*N 

IF(  IRAN* IJI.GT  >1)  CO  TO  36 
IF  dTETA.  EO.-l  .ANO.  lX(J).EQ.D)  CO  TO  38 
IVIUI- IXIJdITETA 
CO  TO  35 
38  lYIJ)-IX(j) 

35  COST  1NUE 

CALCULATE  ITS  EIRORS/X) 

VY- VNORS  I XLACt  IV) 

COMPARE  WITH  THE  INCUMBENT  SOLUTION 
IFIVV.LT.VXI  CO  TO  59 
30  CO NT INJE 

IFI ITETA.EO.-ll  SO  TO  90 
IT  ET  A*-l 
CO  TO  100 
50  00  60  1-1, N 
60  1X01  II-  IV  (II 
VXO-VV 
CO  TO  119 
90  00  70  I »  l »  N 
70  IXOdl-IXII) 
vxo-vx 
110  RETJRN 
END 

SUBROUTINE  AR  AMAI  XL  AC,  1 X,  VX,I  XO,  VXO) 

REDUCES  THE  SET  OF  POTENTIAL  CANDIDATES  BV  CONTRADICTION 

DIMENSION  IX (1001,1X01 100 1, IVI  100  I, IZ ( 1001 , 1  SIR  1 1 00) ,  JARR I  100 1 ,M0 
*CI 100) 

COMMON  N ,  M/  BLO  CX  5/ 1 T ET  A/B  LOCK  8/NO  ZE  ,  J  ART  /BLOC  K6/WLAM 
COMMON  /COUNT/ NITER, NBAS  153), NS EAI 50). NT AR15 3),  NARAI53I ,NPA1(  53 ) , 
*NPA2<  57I.NPA3I  59),NVN0R 
NARA(NIT  ER)»NARAIN1TER)*1 

WL  AM  -  XL  AG 
IT  ETA-1 
500  I C CUN-0 
N9ZE-0 
IKODE-O 

A 00  00  1  0  I -1  ,N 

IF  INOZE.EQ.O)  GO  TO  19 
DO  15  I  I  -1  •  NOZ  E 
IF  I  JARR  (  m.EQ.I)  GO  TO  13 
15  CONTINUE 
19  00  20  J-l «N 

IFI  J  .E3.I )  GO  TO  25 

IFIITETA.E'J.-l.AND.IXIJ  J.EO.O)  G3  13  25 
IF  INCIZE.  EO.O)  GO  TO  29 
03  26  J  J  -  l.NOZE 
IF IJARRI JJ).EQ J  )  GO  TO  25 
28  CONTINUE 
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247 

248 

249 

250 

251 

252 

253 
754 

255 

256 

257 

258 

259 

260 
2A1 
262 

263 

264 

265 

266 
267 
768 

269 

270 
771 

272 

273 

274 

275 

276 

277 
27B 
279 


280 

781 

282 

283 

284 
235 
7R6 
287 
?83 

289 

290 
281 

292 

293 

294 

295 

296 

297 

298 
293 

300 

301 

302 

303 
3P4 


29  1Y  IJ  )■  IX  (J  1+  IT  ETA 
00  TO  20 

25  IY<J)«IXCJ) 

20  CONT  tMJE 

ALFA2-0IF1  XlAG.IY.I  ) 

IF  (ALFA2  I  11,12,12 

11  ICOUN-ICOUMl 
ALFA1>04F{  XLAC.IX.It 
RATIO'AL  FA?/  ALFA1  . 

ISIRClCQUNl-l 

IFt  ICOW.EQ.il  GO  TO  41 
1FIRATI0.LT.RMAX)  GO  TO  19 
41  JcNAX-RATIC 
JJ-I 
03  TO  10 

12  NOZF«N0ZE*l 
JARMNOZEI-I 

10  CONTINUE 

lFlIK33E.CC.il  GO  TO  110 
IF(  IC0Uh.EQ.31  GO  TO  110 
03  33  I«1,N 
IF I1.E3. JJ)  GO  TO  33 
CO  TO  35 

33  IF ( IT  ET  A. fcO.-l  .AND. IX ( 1 1.EQ.0 I  GOTO  35 
I  Y(  I  )  ■  IX(  IIMTETA 
GO  TO  33 
35  IYin-IXII) 

30  C3NT  INUE 

IF  I 1C0UN.E0.1)  GO  TO  63 
N-JC  *0 

00  50  J-l,  ICOUN 
IF  (  IS  18  C  J  1  .EQ.JJ  1  GO  T3  SO 
JL-ISIRUI 
C 

C  PROVE  BY  CONTRADICTION  THAT  JL-CONPOMEN T  CAN'T  SE 
CALL  PARTI (IY.JJ.JL , IZ, ICOOE) 

IH ICOOE. EO. 01  GO  TO  55 
DO  59  Jl-l.h 

59  IYUlIzIZlJll 
GO  TO  50 

55  NUC*N*JC»1 
MOUC  IN'IC)  *  J 
50  CONTINUE 

IFINUC.EC.JI  GO  TO  60 
IFlNUC.GT.il  GO  TO  40 
JUC*KOUC(Nl'C) 

IF1 1TETA  .EC.-l .ANO. IY1JUC1.  EQ.01  GO  TO  60 
IY1JUC1*  IYIJUC  1  I T  E  TA 
GO  TO  63 

40  CALL  PART2 1NUC  »MOUC , IV, XLB) 

IF1XL8.GE.VX1  GO  TO  100 

CALL  PART311Y,VX,JJf N'JC, MOUC,  IZ) 

00  61  J1«1,N 
61  !Y  (Jll- IZ<  Jll 

60  VY »VN03S 1XLAG,  IY) 

IF 1VY.GE.VX1  GO  TO  100 
VXO-VV 
DO  65  J«1,N 
65  IXOUI’IYUI 
GO  TO  300 


> 
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100  1C9UN*IC0'JN-1 

IF(  ICOUN.EQ.9)  63  T3  110 

N0ZE«NTZE*1 

JARR I N3Z  E)  03 

IK00E*1 

on  to  490 

110  IF  1 1 TETA. EQ.-l I  GO  TO  233 
IT  ETA—  l 
CO  TO  530 
?00  00  70  l«l,N 
70  ixom«ixm 

V  XO*V  X 
300  RE7URN 
END 

SUBROUT  1  NE  P AAT1 (IX , J J , JL ,  IY ,  I COOEI 

PERFORMS  A  TEST  BY  SIMPLY  EVALUATING  AN  UPPER  SOUN)  FOR  THE  LINEAR 
IF  ICOOE-1  ,  TEST  SUCCEE9ED 
IF  1CD0E*0  TEST  FA  I  LEO 

DIMENSION  1X(103),IY1130),IZ(133)  ,XCI  50),  XAI5  3) ,  J  ARRt  100) 

C3PM0.N  N  ,  M/BL0CK5/I  TETA /BLOC K6 /ML AM/B L0CK8/NQZE  ,  JARR 
COMMON  /COJNT/ NITER,  NBAS!  50),NSEA(  50)  ,NTAR(  53) ,  NARAI5  3)  ,NPA1I50)  , 
• NP  A2  C 53)  ,NPA3 (53) ,NVNOR/TAB/POISONC53 ,133  I 

NPA1  (NIT  ER  )*NPA1|  NITER  )♦! 

Ml «M+ 1 
XLAC-HLAM 

BO— OIFIXLAG,  IX.JLI 
IF (B0.CE.3.0)  CO  TO  59 
OC  13  I*  1,  N 
IF (l.EQ.JJ)  CD  TO  11 
IZIII-IXI II 
QC  TO  13 

11  121  I)  I X ( 1 1-1 TETA 

13  CONTINUE 
.  31  —DIF ( X L SC,  IZ.JJI 
DD  29  1  * l ,  Ml 
IXS-IX(JL)«I 
T1=P0F0P( IXS,  JL) 

T2  =  P0 1  SON ( IXS, JL) 

XC(I)— T1/T2 
IXP*1X( JJ)«1*1 
P1=P0F3P ( IXP.JJI 
P2*P0  ISONI  IXP, JJ) 

XA  ( I )  *— PI  / P2 
20  CONTINUE 
TER1 “l . 

Tc  R2  *1  • 

IFIBI.LT. 3.0)  CD  TO  22 
WIO  .3 
SO  TO  23 

22  Ml* I1.1)*BD/Bl 

23  TERH«K1 *B I 
OC  30  1*1,  Ml 
WI*XCI I )-Ml*XA( I ) 

IF(WI)  39,30,35 

35  03  49  X- 1,N 

DD  37  KK» l.NOZE 
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357 

IF(JARRIKKI.EO.K)  SO  TO  36 

358 

37 

CONT  IN  J E 

359 

IXK*1 XI K)+ !♦! 

363 

TER1-TER 1*P0I SONI I  XK.KI 

361 

CO  TO  38 

362 

36 

IXK-IXI K>*1 

363 

TER1-TER  1*P0IS0NI  IXK,K> 

36V 

39 

IFIK.ES.JJ)  CO  TO  39 

36$ 

1 Y  K*  I XI  K 1  ♦  1-1 

366 

GO  TO  41 

367 

39 

1YK*  IXI K  )  *  1*1 

363 

41 

TSR2*TER2*P0IS0N<  IYK.KI 

360 

40 

ONT  INUE 

373 

Ul*TFRl-TER2 

371 

TERN*  TERM  ♦,JI*WI 

372 

33 

CONT  1NU  E 

373 

WM  “TERM- 80 

374 

IC0DE*0 

375 

IFIViNl  50,55,  55 

376 

5D 

ICQDE'l 

377 

IF  I  IT  ETA  .EQ.-l  .AND.  IXIJLI.EQ.O)  CO  TO  55 

378 

IX  l  JL )  ■  IX  I  JL  )  +  ITETA 

379 

55 

DO  60  1*1,6 

383 

60 

iYin-ixm 

381 

RETURN 

382 

END 

383 

SUBROUTINE  PART2INJC.M0UC,  I  Y,  XLBI 

c 

C  IF  PARTI  FAILS,  THIS  FINDS  AN  UPPER- BOUND  FOR  THE  SUM  OF  FIRST 
C  DIFFERENCES  OF  THE  UNOETERM IMED  C3EFFI? IENTS 
C 


384 

DIMENSION  IYI  lDOl.MOlCI  109)  ,IZ(103) 

3B5 

C'J KMON  N ,  M/  B  L  0  CK  5/  I T  ET  A  /BL  3 CK  6  /WL  AM 

386 

* 

COMMON  /COUNT/MTER  ,N3ASI53),NSEA|53),NTARI5?  I,  NARAI5D  ) ,  NPA1 1  53 
♦  NPA2I50  I.NPA3I  50), N  VN3R 

387 

w 

NPA2  INI  TER)*NPA2I  NITER  )♦! 

388 

SJM* 0.0 

389 

DO  10  1*1, NUC 

393 

IXK*MOUC<  11 

39i 

IXKl*  IXK-1 

35? 

IF  I  IXK1  . EQ.O)  GO  TO  15 

393 

00  2C  J*  1,  IXK1 

39V 

IFIITETA.  EO.-l  .AND.  IYIJ  ).EQ.9)  G3  T3  25 

395 

IZIJ)=IYIJ)*ITETA 

396 

GO  TO  20 

397 

25 

IZ(J)*irlJ) 

358 

23 

CONT  I  hD  E 

399 

NXK* IX< 

40? 

GO  TO  35 

4-«l 

15 

NXK*  1 

4.2 

35 

00  30  J-NXK.N 

4?3 

33 

I  Z  (  J  )  *  l  Y 1  J  ) 

4C4 

TcRMOIF|WLAM,IZ,IXK|- 

435 

IF  ITERM.LT.D.O)  GD  TD  40 

4C6 

TERM'D. 0 

407 

40 

SUM*  SUM  ♦  TERM 

4?8 

13 

CONT  INUE 

409 

XI 8* SUM* VNORSIHLAM, IY) 

I. 
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410 

RE7URN 

411 

ENO 

412 

SUBROUTINE  PAR  T31 IX  ,VX  ,JJ,NUC ,  MOUC,  U  1 

c 

C 

PERFORMS  THE  SAME  TEST  AS  »A»T1,  BUT  SETS  UP  THE  LINEAR  PROGRAM  EX 

413 

C 

DIMENS  ION  1X1  100),  1 Y (  10 0),  I Z(  1031  ,JARRt  1301  ,MOUCt  1001  .K0UCI103)  , 
•10,10) 

414 

COMMON  N.M/BL0CK5/1 TETA/B LOCKB/NOZE , J ARR/BLCCK9/NR0W, NC0L/BL0CK6 
•'.AM 

41 5 

COMMON  /COUNT/NITER, NBASI50), NSEAf 50), NT AR{ 531,  NARAI50) ,NPA1 (50) 
*NP A2I  53 ) *NPA3I 50) .NVNOR/TAB/POISONC 53,100) 

416 

NP  A3  (NIT  £R  )»NP  A3( N  17ER  )4l 

417 

L 

NUCB-NUC 

413 

X  L4G-WLAM 

419 

NROM-M-2 

426 

NCOL-MO 

421 

NC0L1-NC0L-1 

422 

00  9  I* 1 »N 

423 

1F1 I.E3.JJ  I  GO  TO  8 

424 

IZ ( I)-IXt I 1 

425 

a  TO  9 

426 

8 

12  (I)- IX  (  1,1-  IT  ETA 

427 

9 

CONTINUE 

426 

00  1G  I SAY-1 ,  NUCB 

423 

JL-HOKI  ISAY  1 

430 

Yl  1 ,1)*0. 3 

431 

DO  11  J-2.NC0L1 

432 

TCRl-1.0 

433 

TER  2-1.0 

434 

03  40  K-l.N 

435 

00  37  KK-l.NOZE 

436 

IMJAJtai  KKI.EQ.K)  SO  TO  36 

437 

37 

CONT 1NJ E 

438 

I XK-I Z ( Kl  +  J 

439 

GO  TO  38 

440 

36 

IX  K«I2IK1*J-1 

441 

38 

TER1-TER1- POI SCNI JXK,K) 

442 

If  (K.EO.JJ  1  GO  TO  39 

443 

IYK»I21K1*J-1 

444 

GO  TO  41 

445 

39 

1YK-  12  (K  )*J 

446 

41 

TER2-TER2* POISON! IYK.K) 

447 

40 

CONTINUE 

448 

11 

Y(1,J  >*T  ER1-TER2 

449 

Yl  1  ,NCOU  —  DIF  (XL  AG,  12,  JJ) 

450 

00  15  I-2.NR0W 

451 

IXS-IXl  JLMI-1 

452 

Ul  ,1)  -  PDF  OP  (  I  XS,  JL) /POISON  (IX  S,  JLI 

453 

IXP-  IXUJ  MI-2 

454 

]  XR* I X I J J  M  1-1 

455 

15 

Yl  I ,  NCOL  )*r'DFQP(IXP,JJ)/PO)SON(IXP,JJ) 

456 

DO  25  I- 2, NROrf 

457 

DO  25  J-2  , NC0L1 

458 

iriJ.E9.ll  GO  TO  27 

459 

Yl I. 41-0  .0 

460 

GO  TO  25 

461 

27 

YC  1,  J  )*  -  1  .9 
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I 


AM 

25 

CONTINUE 

463 

BO—DIFI  XLAG,  IX.JLI 

A6A 

00  213  1*  If  K'ROM 

465 

210 

CONTINUE 

466 

CALL  XLPROIY.VOBJI 

467 

FARK-VOBJ-BO 

468 

IF IFARK.GT.O.O)  GO  TO  10 

469 

IFIIXIJL  J.EO.O.AVD.ITETA.EO.-l)  GO  TO 

473 

IX ( JL I* !X( JL I  * 1TETA 

471 

50 

KUC*N«IC-1 

472 

IFIKUC.EQ.jl  GO  TO  10 

473 

NKUC-0 

474 

03  5  5  J-l.NUC 

475 

IF  IMOJCI  J)  .EO.JL)  GO  TO  56 

476 

1FINKUC.EQ.1I  GO  TO  57 

477 

KOUCI JI-MOUCI J 1 

473 

GO  TO  55 

479 

56 

NK  UC  *  1 

480 

GO  TO  55 

48* 

57 

JI*J-1 

48? 

K3UCIJ1) *M0UCI J) 

483 

55 

CONTINUE 

484 

CALL  PAR  T2  IKUC,  OUC  .IX.XLB) 

465 

IFIXLB.GE.VXI  GO  TO  503 

486 

NJ  C*KUC 

487 

00  60  1*1*  MUC 

488 

60 

H3UCI II =KCUC ( 1 1 

469 

10 

CONT  INUE 

490 

500 

DO  70  1*1  .N 

491 

70 

IYI  II-IXCII 

492 

RETURN 

493 

ENO 

494 

FUNCTION  0  IF ( XLAGt I  X«K  1 

c 

c 

CALCULATES  TFE  FIRST  OIFERENC ES 

495 

V 

DIMENSION  IX  1  1001 

456 

COMMON  H,M/BL0CK1/C 1100 1/T AB/POIS 0NI5 0 

45  7 

DOUBLE  PRECISION  PRO, SUM 

493 

Ml=M4l 

459 

SUM»0.0*00 

500 

03  20  J  1*1, HI 

531 

PR0*1 .0+00 

502 

03  30  1*1, N 

503 

IX  J-  IX  (  1  I  ♦  J 1 

504 

IFI  I.EQ.K1  GO  TO  25 

505 

PR0*PR0*P0iS0Nl IXJ, I 1 

506 

GO  TO  33 

507 

25 

PPO  =  PRO*  POFOPI IXJ, 1 1 

508 

30 

CONT  INJE 

*09 

23 

sjm»sjm+pro 

;10 

IXM*IXIK1*MM 

511 

01  F*  -S’JM-1  .*P0  ISONI  IXM,KI*XCAG«C(KI 

512 

RETURN 

513 

ENO 

514 

FUNCTION  VNORSIXLAG,  I X  1 

c 

c  evaluates  the  function  e<nors/xi*l.c.x 
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■ 


515 

516 

sir 

516 

51* 

520 

521 

522 

523 

524 

525 

526 
522 

528 

529 

530 

531 

532 

533 

534 

535 

536 

537 
533 

539 

540 

541 

542 

543 
5  44 

545 

546 


DIMENSION  1XC100I»P Ml  139,29  I •  PRO! 23)  . 

C3MM0N  N,N 

COMMON  /COUNT  /  Nil  EX  (M8  AS  I  50  )<  NS  Ell  50)  ,NTA  RI53) ,  NARAI59)  ,NPA1I  531 
♦  V>A2<  53) .NPA3I  SO) ,NVN0R/TAB/P01S0NI53 (123 )  ' 

DOUBLE  PRECISION  SUM,  PROt  TOPL  •  PM ,  XM 

NVNOR  ».NVNCR*1 

EPS*  .03001 

Ml*M*l 

SUM*0.  3*  33 

DO  20  JJ«1.M1 

PRO  I JJI* 1 .003 

DO  25  11*1  ,N 

IXJ*  IXI  IIMJJ-1 

IF  I  IXJ.GT  .9)  GO  TO  24 

XX*PDF3PI  IXJ.I1) 

GO  TO  25 

24  XX “POISON I IXJ, II) 

25  PRO  I J  J  >  =PXO(  J  J  I* XX 
23  SJM=SJM*1  ,-PROIJJ  ) 

TOPL *3. D, 39 
DO  30  1*1, N 
K*  3 

40  MXK“IXd)*K*2«-K 
WX*PDFOP(MXK, I  ) 

WJ “FLOAT (K+1)*WX 
IF (WX.LT. EPS)  GO  TO  33 
TOPL»TOPL*WN 
K«K4l 
GO  TO  40 
33  CONTINUE 

VNORS*SUM.»TOPL  ■•XL  AG *03  TP1 1 IX) 

RETURN 

END 


547 


548 

549 

550 

551 
55? 

553 

554 

555 

556 

557 
556 

559 

560 

561 

562 

563 

564 

565 
564 

567 

568 


SUBROUTINE  TAB  LEI PO ISON ) 

C 

C  CALCULATES  TFE  CUMULATIVE  POISSON  SUNS 
C 

DIMENSION  POI SON!  5C, 133) 

COMMON  M , M/ BL0CK2/XL AM  1133) 

DOUBLE  PRECISION  EPS, TERM, SUKX 

EPS* 1 .0-06 

03  13  J*1,N 

TE  RM* 1 . D+0  3 

SJ“X“ 1. D*03 

DO  20  1*1,50 

TFRM“TERM*XLAMIJI/FLOATIl  ) 

I<*  I 

1MTERM.LT.EPS)  GO  TO  25 
SUMX*  S'JK  XdERM 

POISONII,  J)=EXPI-XLAM(J  ))*SUMX 
IF  (POI  SCNU,  J)  .GT.l  .0)  POISONI  I,  J)“l  .3 
23  CONTINUE 

25  I F  (  IX.EO.5'3)  GO  TO  13 
03  30  1  * ■ X,59 
33  POISONd,  J  1*1.0 
10  CONTINUE 
RETURN 
E*4D 


4 


'*1 
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569 


FUNCTION  FCFOPIK.il 


C 

C  CALCULATES  TFE  INDIVIDUAL  POISSON  TEAMS 
C 


5T0 

COMMON  /8LCC K2/XLAM 1 1 301 

571 

DOUBLE  PRECISION  SUM.HW 

572 

SUM-0. 0*00 

577 

IFIK.LT.2t  CO  TO  10 

575 

DO  15  J*2,K 

575 

SUM* SUM^ALOGI FLOAT!  Jll 

575 

15 

CO  NT  1NUE 

577 

10 

KW—  XLAMI  I  )-K*ALOGIXL  AM  (ID-SUM 

578 

IFIWW.LT. -15. 01  SO  TO  20 

579 

POrOP-DEXPIWW) 

580 

CO  TO  25 

581 

23 

P0F9P  0.0 

582 

25 

RETURN 

583 

END 

585 

FUNCTION  OOTPRI  IX) 

CALCULATES  THE  DOT  PRODUCT  C.X 


585 

DIMENSION  IX  (  100) 

585 

COMMON  N.M/8L0CK1/C  (100  ) 

587 

SJM-0.0  • 

588 

on  10  i-i.n 

589 

XX-FLOAT ( I X( 1 1  1 

590 

13  SUM* SUM *XX*C  (  I ) 

591 

03  TPR-SUM 

592 

RETURN 

593 

END 

595 

SUBROUTINE  S I RALA (X . I R) 

C 

C  THIS  ARRANGES  AN  ARRAY  OF  ELEMENTS  IN  ASCENDING 
C  IRIK)  IS  THE  RANK  OF  ITEM  K  INTHE  LIST 
C 


595 

DIMENSION  XI  1301,  1R  (  100),  ISEOI  100),  Yt 

596 

COMMON  N, M 

597 

03  9  1-1, N 

593 

3 

IS  EQI  1)*  I 

599 

DO  10  I -l ,N 

600 

XM  IN*  XI  1 1 

601 

IM  IMJ-l 

602 

IMIN-ISE0I1) 

603 

Nl-N-I*l 

605 

00  20  J  * l , Nl 

605 

tF(XIJ)  .LT.XMINI  SO  TO  21 

606 

GO  TO  20 

607 

21 

XMIN*XI J) 

(>08 

IM  IN*  ISEOI J) 

603 

IM  IN  J*  J 

610 

2  3 

CONTINUE 

611 

HI  ININ  l-I 

612 

I F IN  1  .60.11  GO  TO  13 

613 

N2  *N1-1 

615 

D?  30  J*  1,N2 

615 

IF(  J.LT  .  I M I NJ  )  GO  T  3  33 

616 

Jl-JM 

0R3ER 
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617 

ISEOUI-ISEOiJll 

618 

Vi  J)«Xi  Jl) 

619 

GO  TO  29 

62} 

33 

VI J)-X(J) 

621 

29 

XI  JJ -Yl  J) 

622 

33 

CONT INJE 

623 

ID 

CONTINUE 

624 

RETURN 

625 

END 

626 

SUBROUTINE  XLPRO(Y,VOajI 

e 

c 

c 

c 


SOLVES  A  LINEAR  PROGRAM  OF  THE  TYPE 
MAX  CX  ST  AX<B.X>0 


627 

DIMENSION  Yi  10*23  It  Xt  10i  20 1 «  ISNBI  lDI.ISBVi  lDl.HRDUilOl 

628 

COMMON  / BLQCK9/NR0W » NC QL 

629 

V 

is  av  (ii-o 

63D 

I SNB  ( 1  i  “NCC'L+1 

631 

03  10  I-2.NR0W 

632 

10 

1SB/I  II-  I-2+-NC0L 

633 

00  11  J-2.KC0L 

634 

11 

ISN8(JI«J-1 

635 

NR0M1-NR0U-1 

636 

NCOLl-NCOL-l 

637 

NC3L2-NC0L  *1 

638 

NCOL3-NC0L*NR0Hl 

639 

03  12  I-l.NROW 

64  • 

00  12  J-NCOL  2* NC3L3 

641 

IM=I-1 

642 

JM- J-NCOL 

643 

IF(JM.EQ.IM)  GO  TO  13 

644 

Y(I, JI-0.0 

645 

G3  '0  12 

645 

13 

Y(I,J)-1.3 

647 

12 

CONTINUE 

643 

100 

NNC-0 

649 

NPV-3 

65D 

NPC-0 

651 

NN  V-3 

652 

00  15  J-2.NC0L 

653 

IF CY(1, Jl .LT.O.O)  GO  TO  16 

654 

NNC*NNC*1 

655 

1 F ( NNC . F  0  •  NC  OL 1  1  GO  T3  533 

656 

GO  TO  15 

657 

16 

NPC-NPC*  1 

658 

IF(NRC.EO.l)  GC  TO  17 

6  59 

IF IY(  1  *  J I . GT.YM1NJ  SO  TO  15 

660 

17 

YM  IN- Y(l* J » 

661 

JM IN- J 

662 

15 

CONTINUE 

663 

IE V*  ISNSiJMIN  I 

664 

NM  V  A  •  1 

665 

03  2  0  I-2.NR0H 

565 

IFiY  (  I*  JMIN)  .LE.3.0  J  G3  T3  25 

667 

NYV-NNVM 

668 

7ETA* YC  I,  1  l/Yi  I.JNINI 

669 

IF |NNV. EO.l)  GO  TO  23 

670 

IF (TETA-TM1NI  23*20*20 
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671 

672 

673 
67V 
675 
67b 
67  7 
67ft 
679 
667 
661 
6B? 
663 
686 

695 

<  e6 
667 
68S 
#69 
690 
651 
692 
653 
666 

696 
696 
657 
699 
669 
70C 

701 

702 

703 
706 

705 

706 
7C7 
703 

709 

710 

711 

712 

713 
716 

715 

716 

717 

718 

719 
722 

721 

722 
7  23 
726 

725 

726 


23  TM1N-TETA 
1MIN-I 

M-tOuINMVA  I »I Ml N 
GO  TO  23 

25  NPV«NPV*1 

iriN»v.ro.:JROtfll  G3  T3  653 
23  CONTINUE 

03  26  I-2.NR0W 

IF  17(1,  JH  Ml. LE. 0.0  I  GO  TO  26 
TETA-VII  ,1  l/VI  l,  JMINI 
IF ( TETA .GT.TM1NI  GO  TO  26 
IF(  I.EO.  IMJNI  GO  T3  26 
NM VA«NNVA*1 
MROV  (NMVAt-1 

26  C3NTINJE 

IFINMVA.EO.il  GO  TO  230 
03  80  Is  l.NKVA 
IK*  MROW  (II 
JPOSO 

00  80  J*NCQL  2.  NC3L  3 
X(1K,JI-Y(  IK. J  l/Y (IK.JMINI 
IF(X(  IK,  Jl  .GT.O.OI  GO  TO  85 
G1  TO  30 
85  J?jS*JPOSM 

1FUP3S.GT.1I  GO  TO  83 
IF ( IK .GT .MROW (111  GO  TO  87 
Gj  TO  93 

87  1FIJ-JKMHI  83,93.93 

90  IF  (XIIM1N,  JKMINI.LE.XI1K,  Jll  GO  T)  S3 
93  1 9  I N ■ IK 
JKMIN»J 
80  CONTINUE 
200  IDV*  1  S3  V(  I  MIN  I 
1SNBIJM IN  1* 10V 
IS9VI IMI N ) *1 EV 
03  60  I  *  1 ,  NP-  0  W 
03  60  J* 1, NCOL  3 

IF I I .  EO. 1  MINI  GO  TO  55 

IFIJ. E3.JM  IN  I  GO  TO  63 

XI i , J1*Y  I  l.J  l-IYI l,  JMIN  »*rl ININ, JI/YC  1MIN, JMINI  1 
GO  TO  60 

55  IF  IJ.E3.JM1N}  CO  TO  57 

XI 1 ,J|*Y(  l.J )/Y( 1M IN, JMIN) 

GO  TO  60 

57  XI  l.J  1*1  ./  Y  |  I,  J} 

GO  TO  60 

63  XI  I,JI*-Y(  I.JI/YI  IM1N,  JMINI 
63  COUTINJE 

013  73  1  *1  ,  6P0W 
DO  70  J*1,KC0L3 
73  71  I,  JI=X  I  I.J  I 
GO  70  103 
650  -3ITEI6,  25:i 

250  F0FMATI5X, ‘ONBCUNOEO  SOLUTION*! 

500  RETURN 
PIO 


iEMRV 
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THE  SADDLE  POINT  SEARCH  ALGORITHM 


PROBLEM  NUMBER  1 


NUMBER 

OF  ITEMS- 

5  PUOGET-* 

25000.03 

I  NI  Tl  AL 

OPTIMAL 

ITEM 

COST/JNIT 

DEMAND 

SOLUTION 

SOLUTION 

1 

2980 .00 

2.10 

2 

2 

2 

1751 .00 

1  .50 

2 

2 

3 

462.00 

1.20 

2 

3 

4 

1500.00 

5.00 

6 

7 

5 

345.00 

3.50 

6 

6 

EIN9RS/XI-  1.06282  COST  IX)* *  2S418.00 

t 

LAGRANCE  MULT  IPIISR-  0.000099 

•  OF  ITERATIONS-  6 

•  OF  NORS  EVALUATIONS-  141 


t  OF  STEPS  IN  EACH  ITERATION 


ITERATION  f 

STEP  1 

STEP2 

STEP3 

STEP4 

PARTI 

P4RT2  PART3 

1 

1 

1 

1 

1 

5 

0 

0 

2 

1 

5 

5 

2 

0 

0 

0 

3 

1 

4 

4 

2 

3 

0 

0 

4 

1 

2 

2 

1 

1 

0 

0 

5 

1 

2 

2 

1 

0 

0 

0 

6 

1 

2 

2 

1 

2 

0 

0 

STATEMENTS  EXECUTED-  273 

42  7 

CORE  USAGE 

OBJECT  CODE- 

31  504  BYTES,  ARRAY  AREA 

-  40572 

BYTE 

DIAGNOSTICS  NUMBER  OF  ERRORS-  Ot  NUMBER  OF  WARNINGS-  0.  * 

CCPP11E  TIME-  0.52  SEC,  EX  ECUT  ION  TIME-  5.02  SEC,  12.27.33  TIE 

CtSTOP 
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APPENDIX  B.l 

DETAIL  DOCUMENTATION  OF 
GREEDY  ALGORITHM 
COMPUTER  PROGRAM 


APPENDIX  B. 1 


i 


DETAIL  DOCUMENTATION  OF  GREEDY  ALGORITHM 
COMPUTER  PROGRAM 


General 


In  addition  to  the  explanatory  comment  statements  that  are  provided 
in  the  program  listing  of  the  algorithm  (See  Appendix  B.2),  formulations 
of  computational  formulas,  and  detailed  statement  by  statement  interpreta¬ 
tions  are  provided  in  this  appendix  for  several  critical  parts  of  the  pro¬ 
gram.  Since  some  of  the  subroutines  are  either  essentially  the  same  as 
those  of  the  D029  program,  or  rather  simple  and  straight  forward,  the 
effort  below  is  devoted  mainly  to  three  subroutines:  TABLE,  OPTIM,  and 
MARGAN.  A  flow  chart  describing  the  relationship  between  OPTIM  and  MARGAN 
is  also  provided. 

Terminology  and  Notations  Used  in  Notes  on  Subprograms 


a*  -  NAP(I) 
bi  =  B0(I+1) 

q±  =  KIT(I) 

w 

Fi(x) 

ITOT 

MUP 

NAP(I)=a1 

NMAX± 

NOPS  »  UE 
NORS  =  NMCS 
POISON  (I,J) 


number  of  applications  of  item  i 

conditional  operational  rate,  given  i  planes 
available  for  canibilization 

quantity  of  item  i  in  the  kit 

weighting  factor  of  AE(SDO)  used  in  r±  calculation 
probability  of  demand  <_  x  for  item  i 
total  number  of  items  considered 
upper  bound  of  subscript  of  bi 

maximum  possible  number  of  demands  for  item  i 
number  of  planes  in  squadran 

number  of  Not  Operationally  Ready  Supply  aircraft 
cumulative  Poisson  demand  Fi(j-l)  for  item  i 
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r^“RATIO(I)  Weighted  average  improvement  of  E(NORS)  &  E(SDO) 
per  dollar  invested 

SDO  number  of  Stock  Due  Outs,  or  backorders 

S(N)  probability  of  NORS  <  N-l 

UE  -  NOPS 

Notes  on  Main  Program 

This  main  program  contains  the  various  subroutine  calling  statements 
shown  below. 


Notes  on  Subroutine  INPUT 


This  subroutine  reads  input  data,  calculates  average  demand  rates,  and 
obtains  KIT  5718  and  its  target  value.  The  code  is  essentially  the  same  as 
the  code  for  the  D029  algorithm.  LRU-RR,  LRU-RRR,  SRU-RR,  and  SRU-RRR  are 
differentiated  and  processed  in  different  manners  as  required. 


Notes  on  Subroutine  TABLE 


This  subroutine  calculates  the  individual  cumulative  Poisson  probabilities 
for  each  item. 

AX  “XI 

Poisson  probability  is  Pj(x)  *  -I— ^ -  ,  x  -  0,1,2,...  and  I  =  1,2,...1T0T 

hence,  Pj  (x+1)  =  PI  (x) 

e-g-  Pi  (0)  =  e‘AI,  Pl  (1)=  At  •  Pi(0),  Pt(2)  Pr  (l),...etc. 


Then  the  cumulative  Poisson  probability,  Fi(x),  can  be  expressed  as  follows: 
x  XI  e~AI 


Fx(x)  =  I  — - -  *  Fj  (x-1)  +  P!(x) 

j=0 


By  employing  the  two  dimensional  array  POISON  (500,60),  Fi(j)  of  item  i 


can 

be  stored 

in  POISON  (I,  j+1),  since  POISON  (I,  1)  =  Fj  (0)  for 

The  main 

segment 

of  code  for  this  subroutine  is  shown  below. 

ISN 

0007 

TE«t»5.  0-  12 

ISN 

0008 

00  10  I<l«I TOT 

ISN 

0009 

nn««0BL£<  XNU<  II  1 

ISN 

0010 

TEH»*0£XPt-MNN) 

ISN 

0011 

POISON*  I«  1  »=TERH 

ISN 

0012 

DO  20  0-2*60 

ISN 

0013 

I FI  TERM.LT.TERL.ANJ.PQISONII.J-  11  *GT • 

ISN 

00  IS 

TERP>TERH«NMM/DFL04  Tl  J  -  l » 

ISN 

0016 

POISON!  I*  J»«P0IS0NI  1*0-1  1-tlERM 

ISN 

0017 

20 

C  ON  T I  NO  E 

ISN 

0018 

GO  TO  10 

ISN 

0019 

30 

00  40  JX« J*  60 

ISN 

0020 

40 

POISON!  I*  JX  »«  1*D  +  0J 

ISN 

0021 

10 

CONTINUE 

ISN 

0022 

RETURN 

ISN 

0023 

end 

Line  7  gives  the  criteria  for  stopping  further  computations  of  the  upper- 
tail  of  the  Poisson  distribution.  At  lines  10  and  11,  POISON  (I,  1)  =  Ft  (0)  = 


e-AI.  Line  13  checks  if  the  stopping  criteria  is  met.  Lines  15  and  16  calculate 


all  Pj(x)  and  Fj(x)  for  each  item.  Lines  19  and  20  assign  to  all  the  re¬ 
maining  cumulative  probability  terms  the  value  of  1. 

Notes  on  Subroutine  INKIT 


This  subroutine  is  essentially  the  same  as  that  of  D029.  It  appropriately 
increases  KIT  quantities  to  ensure  that  the  probability  of  any  item  being  the 
cause  of  half  a  squadron  being  down  is  less  than  0.50. 

Notes  on  Subroutine  OPTIM 


This  routine  calls  marginal  analysis  subroutine,  MARGAN,  and  updates  the 
parameters  b^'s.  It  terminates  if  the  new  set  of  bj's  are  sufficiently  close 
to  the  values  of  the  previous  set. 


B0(J)  «  bj-i 


ITOT 

X 


ITOT 

=  v  POISON  [KIT(£)  +  NAP(£)  *  (J-l)  +1] 
f  =  l 


The  measure  of  closeness  between  two  sets  of  bj_’s  and  b^'s  is  the  sum  of 

squares  =  E  (bj_  -  b^)2. 
i 

Referring  to  the  code  of  the  main  body  of  this  subroutine  on  the  next 
page  and  the  combined  flow  chart  that  follows,  the  algorithm  can  be  clearly 


explained. 


XSN  0017 


c 


100  CALL  MARGANt  KIT  (PRICE  I 


XSN  0010 


ISN  0033 
XSN  0035 


C 

c 

c 

c 


NXTER*NITER+1 
STEP  2 

UPDATE  THE  PARAMETERS  B 


XSN 

0019 

SUHSG*0.0*00 

ISN 

0020 

DO  10  J*1.MUP 

ISN 

0021 

PR0*1.D»00 

ISN 

0022 

00  20  I *1.1  TOT 

ISN 

0023 

KlTPJ*KIT«X  j- 1  )«NAP(  IM-J 

ISN 

0024 

IFIA1TPJ.GT.60)  GO  TO  20 

ISN 

0026 

PR0*PM0«P0IS0M  X.KITPJ1 

XSN 

0027 

IF<  PRO.  LT.EPSttl  GO  TO  15 

XSN 

0029 

20 

CONTINUE 

ISN 

0030 

15 

5UMSW*SUMSU*<  (B0<  J!-PRO»**2l 

ISN 

0031 

B  04  JI*PR0 

ISN 

0032 

10 

CONTINUE 

C 

C 

c 

c 

c 


STEP  J 

COMPARE  THE  NEN  VALUES  CP  B  MllH  PREVIOUS  VALUES 
TERMINATE  XF  THE  V  ARE  SUFFICIENTLY  CLOSE 

IFiSUMSU.LT.EPS!  GO  TC  30 
GO  TO  100 


ISN.  0036 
XSN  0037 
ISN  0038 


30  CALL  CAlCI  K  IT  tENOHS  .SCO  ) 
RETURN 
END 


Line  17  calls  the  marginal  analysis  subroutine  in  order  to  find  the 
optimum  kit  under  the  current  bi  parameters.  Once  this  is  done,  it  comes 
back  and  starts  to  update  the  bj's. 

Line  24  provides  a  safeguard  against  running  out  of  the  limit  of  the 
Poisson  array.  At  line  26,  b^'s  are  calculated  by  multiplying  all  pertinent 
cumulative  probability  terms  together.  Line  27  gives  a  lower  bound  check  on 
bi  to  save  further  calculations,  since  when  bi  £  EPSQ  =  5x10  11 ,  we  can  treat 
it  as  zero  without  practically  affecting  actual  results.  Line  30  gives  a 
measure  of  the  difference  between  the  latest  two  consecutive  sets  of  bi's. 
Line  33  checks  to  see  if  the  difference  is  small  enough  to  stop.  If  not,  the 
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algorithm  continues  iterations  until  the  above  stopping  criteria  is  satisfied. 
Finally,  before  returning  to  the  main  program,  the  performance  of  the  optimal 
kit  is  evaluated. 


Start  With  Initial 
Kit,  Calculate  All 
RATIOS 


26-33  j _ 

Find  the  Item  With 
Max.  RATIO,  IMAX 


Initialize 


Obtain  the  Price  of 
the  Kit  With  (1jmax+^ 


19-32 


Update  b^'s 


Price  Budget 

s.  ?  ^ 


E(b  -  b'  )2  <  EPS 
i 


37-50  l 


Update  the  Kit,  Price, 
and  RATIO  of  Item, 
IMAX 


1 


Set 

RATIO  (IMAX)=0 


Call  CALC  to 
Evaluate  the 
Kit  Performance 


Return 


MAX  RATIO  >  0 


Subroutine  MARGAN 


Subroutine  OPTIM 


Figure  B.l.  Flowcharts  for  Subroutines  MARGAN  and  OPTIM 
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i 


C 

i 


Notes  on  Subroutine  MARGAN 


This  subroutine  performs  the  marginal  analysis.  It  starts  with  the 
initial  kit  and  calculates  the  RATIO  of  return,  (r^),  per  dollar  invested 
for  each  item  i.  This  ratio  is  the  weighted  average  of  the  improvements 


of  E(NORS)  and  E(SDO). 


=  A|E(NORS)  +  wAjE(SDO) 
ri  Ci 

flTOT 

where  AjECNORS)  =  [  E  ^MUP/>  +  x> 

i=l  *- 

flTOT 

-L  E  (MUPn  +  1) 
1=1 


MUP  I  TOT  -| 

-  E  bj  E  In  Fp(q  „  +  ja„) 

j=o  e=i  *  1  1  J 

MUP  ITOT  -J 

-  E  b-s  E  In  F„(q^  +  ja„) 

j=0  f=l  K  1  1  J 


where  q^  =  +  1  when  £  =  i 


qt  "  q£ 


when  £  t  i 


MUP  r  1 

AjECNORS)  =  E  bi  Un  Fi  (qA  +  1  +  jai)  -  £n  Fi  (qi  +  ja^ 
j=0  L  J 


MUP 


=  E  b ,  £n 


Fi(qi  +  +  i) 


J“0 


j  Fi(<Ii  +  jai) 


MUP  POISON  (i,  qt  +  jai  +  2) 

=  Z  b4  £n 
j=0 


POISON  (i,  qt  +  jai  +  1) 

Since  POISON  (I,  1)  -  Fi  (0) 


E(SDO)  =IT2T  IT  JT[l~F£(qe  +j)] 

£=  1  x=<i£  £=1  j=0 


where  JUP^=  (NMAX  -  1)  - 

«  [KIT  (£)  +  ( NOPS- 1 ) *NAP (£)  +  l]  -  1  -  KIT(£) 
=  ( NOPS- 1 ) *NAP ( f )  +  1 


1 
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"i 

\n 


ITOT  JUPf  r  ,  IT0T  JUP  „  _ 

Then  A^SDO)  =  L  i  C  [l-F„  (q„+j)J  -  E  E  [l-F^(q|  +j)] 

*-i  j=o  *  u  1  j=a 

where  q^  =  q^+  1  when  i  *  i 
q£  =  q^  when  ^  i 


=  JEPl  jl-Fi(qi  +  j)  -Cl-FtCqi  +  1  +  j)]| 
j=0  ' 


JUPi  r  n 

i.e.  A1E(SDO)  =  s  1  [Fi(qi  +  1  +  j)  -  Fi(qi  +  j)] 

j-o 


-  Fi(qi  +  1  +  JUP^  -  F i (q  i ) 

or  =  POISON  fi,  q±  +  1  +  NMAXt  -  1  -  q±  +  l)-POISON  (i,  q±  +  1) 
=  POISON  (i,  NMAXi  +  1)  -  POISON  (i,  KIT(i)  +  1) 


After  comparing  each  ratio  ri  and  obtaining  the  maximum  one,  we  increase 
its  kit  quantity  by  one.  Termination  occurs  if  the  kit's  price  is  greater 


than  budget,  or  the  procedure  is  repeated  to  update  the  kit's  value  and  the 


ratio  until  the  budget  is  exceeded. 


Ratio  Calculation  Segment: 


C  START  WITH  THE  INITIAL  *IT  ANO  CALCULATE  THE  %Ay i-0 

C  PER  OOLLAR  INVESTED  FOR  EACH  I  TEN 

C 


ISN 

0009 

PRICE -PRC  IN 

Tf 

ISN 

0010 

00  10  I  *  1  •  I  TOT 

ISN 

0011 

KITH  IiKITIM  II 

ISN 

0012 

NHAXl-ftlTlI  !♦  (NOP  S«  1 1  IN  API  1 1  42 

ISN 

0013 

1  Ft  NNAX  1*  GE*6  01  NMAXl=60 

ISN 

0015 

I  XS*KITt  I  *♦  1 

ISN 

0016 

SUN  «DBlE!  PAR  1*1  POISON  1 1*  NNAX  ]  l-PO  ISON  I  I  *1  XS  II 

— 

ISN 

0017 

00  20  j«l*MOP 

ISN 

0016 

I  XJ«KITf  I  !♦<  J-1I*N4P(  1141 

ISN 

0019 

1FI  I  X  J»  GE  •  6  01  CO  TJ  2  1 

ISN 

0021 

N«a POISON!  I»IXJ4  1  1/ POISON!  1#  1  X J  1 

ISN 

0022 

GO  TO  20 

ISN 

0023 

21 

1*  0400 

ISN 

0024 

20 

SUrt*SUM4B0t  J!  «0LDCi  *«  1 

ISN 

002S 

10 

RAT  1011 I  *  SUM/ DOLE 1C  OST!  Ill 

Line  9  and  11  start  with  the  initial  kit  and  its  price.  Line  15  increases 


KIT  quantity  by  1,  once  for  every  item.  Line  16  gives  the  vAE(SDO),  which  is 
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the  weighted  improvement  of  E(SDO).  Line  19  provides  a  safeguard  against 
running  out  of  the  POISON  array  limits.  Line  21  and  24  calculate  AE(NORS), 
improvements  of  E(NORS),  and  add  it  to  w-aE(SDO).  Line  25  computes  the  per¬ 
formance  improvement  measurement,  RATIO. 


The  Remaining  Segments: 


C 

C  FIND  THE  MAXIMUM  RATIO 

C 


ISN 

0026 

too 

00  25  I  *  1. 1  TO  T 

ISN 

0027 

I FI  I.EU.ll  00  TO  3  J 

ISN 

0029 

IFI  RATIO!  II.lT.RMAX  *  GO  TO  25 

ISN 

0031 

30 

IMA  X*I 

ISN 

0032 

RMAX*RAT10l  II 

ISN 

0033 

25 

CON  TINUE 

ISN 

0034 

r 

XPR  1C E«f>*  ICE* COST  UMAX) 

c 

TERMINATE  IF  the  T<RGET<  XRE  ACHIEVEC 

ISN 

0035 

C 

c 

IFI  xPRiCE «GT. BUDGET )  GO  TO  5C 

c 

c 

IF  NOT*  UPOATE  THE  KIT  kALOE  A  AO  THE  WAT  10 

ISN 

0037 

KIT  II  MAXI  *KIT  II  MAXI  *1 

ISN 

0038 

P  RI C£  *XPR ICE 

ISN 

0039 

NMA  X1**ITI  JMAXI*I  NJ  PS  -  1  >  *N  API  IM  AX  1*2 

ISN 

00  40 

1FCNMAX1.GT.60I  N  MA  XI  ■  60 

ISN 

0042 

IXS  *K  IT  1  1  MAX  1*1 

ISN 

0043 

SUM*OtfC.Et  PAR  1*1  POISON  1  I  Max  «NA  AX  1  1  -POISO  M  IM  AX  •  JXS  1  1 

ISN 

0044 

00  40  J ■  1 »  HUP 

ISN 

0045 

I  XJ  *K  IT!  IMAxl  *1  J-  II  4N  API  IM  |Xl*l 

ISN 

0046 

iFt  1XJ.  GE.60I  GO  Td  43 

ISN 

0046 

M  *3  PO IS  ON C I MA  X. I X J* 1 1/POlSCNIlMAX .1 X  J I 

ISN 

0049 

40 

SUM*SUM*flOC  Ji  *l)L0GC  MN  1 

ISN 

0030 

45 

RAT  IOC  IMAXlaSUM/OBuEl  CUSTI  IMAxl  I 

ISN 

0051 

GO  TO  100 

ISN 

0052 

c 

50 

RAT 10llMAXt*0.0*00 

ISN- 

0053 

V 

IFI  RMAX.GT.  0.  0*  001  GO  TO  1(0 

ISN 

0053 

RET  URN 

ISN 

0056 

ENU 

Lines  26  to  33  find  the  maximum  RATIO  value.  Line  34  then  calculates 
the  corresponding  price.  At  line  35,  the  budget  is  compared  with  kit  price. 
If  the  kit  price  is  exceeded,  continue  to  check  the  feasibility  of 
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other  possible  kits.  If  not  exceeded,  the  kit  value  and  the  RATIO  are 

updated.  Lines  37  to  50  are  similar  to  previous  segments,  lines  9  to  25, 

except  this  time  calculation  of  RATIO  is  for  the  most  promising  item  only. 

With  the  updated  RATIO  value,  this  item  enters  competition  with  all 

other  items  again.  If  the  kit  with  9imax  +  1  is  infeasible,  then  line  52 

sets  its  ratio  to  zero  in  order  to  eliminate  q.  +1  from  further 

imax 

investigation.  Line  53  then  checks  if  any  other  possibilities  for  consider¬ 
ation  are  left  out.  If  all  possibilities  are  exhausted,  then  imax  =  0; 
hence  the  program  exits  from  MARGAN. 


Notes  on  Subroutine  CALC 


This  subroutine  calculates  the  expected  NORS  and  expected  STOCK  DUE  OUT 


values. 


UE  UE  _  _ 

E(NORS)  =  Z  YP(N0RS  *  Y)  *  Z  Y  [P(N0RS  £  Y)  -  P(N0RS  £  Y-l)] 

Y=1  Y=1 


=  UE  •  P(N0RS  £  UE)  -  Z  P(N0RS  _<  j) 

j=0 

ITOT  UE*1  ITOT 

*  UE  •  tt  Fj^qi  +  UE  -a*)  Z  it  F-^qi  +  ja-t) 
i-1  j«0  i=l 


UE  *  S(UE  +  1)  -  Z  S(j) 


where  S(N)  =  P(N0RS  <  N  -  1) 


ITOT 

w  F1(q1  +  (N-l)*ai) 

i-1 


or  -  *  POISON  (i,  9i  +  (N-l)  *  ai+1) 

i-1 
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i.e.  E(NORS)  =  NOPS  *  S(NOPS  +  1)  -  Z  S(j) 

j=l 

where  NOPS*UE»number  of  planes  in  the  squadron. 

ITOT  NMAXj-1  ITOT  JUP»  _ 

E(SDO)  =  Z  Z  [l-F,,(x)]  =  Z  Z  *[l-F.(q*+  j)] 
£-1  tml  j=0  1  « 

where  JUP^.  =  (NMAXf  -  1)  -q£ 

-  [(q^  +(NOPS-l)*a£  +  l)-l]  -q^ 

=  (NOPS  -  l)*a£  +  1 

e(sdo)  =  I?0T  Ji!P^  [i-poison  (  £,  q«  +  j)] 

*-l  j=0  z 


E(NORS)  Comput 


an  Segment: 


ISN  0009 


001G 

0011 

0012 

OU13 

0014 

0015 

0016 

0010 

0019 

0021 

0022 

0023 

0024 

0023 

0026 

0027 


EPS*5*O-l0 

ElNORS/Xl  COMPUTATION 

SUM*  0*0  + 00 
SINOPS* 

00  10  X*l  (NOPS 
SfK  t*  liDt  ou 
00  20  1  *  1 • ITOT 
1  XJ«K  1T<  1  »♦<  K-l  JoNA  P<  I  »♦  1 
IP!  IXJ.GT.60I  00  TJ  20 
SIX««S<  Kl  *PQlSOM  U  IKJt 
I  Fi  Si  K»  •  LT.  EPS*  GO  TO  15 
20  CONTINUE 

SUM«SUM*S<  XI 
GO  TO  10 
IS  SIX  »«  0.  D+00 
10  CONTINUE 

UNO  RS*DFLOAT(  NOPS  Si  HOPS*  11-SUM 
ENORS-VNOftS 


Line  18  gives  the  probability  of  NORS  j<  N-l,  S(N).  Line  19  provides  a 
check  for  very  small  S(N)  values  to  save  further  computational  effort.  Lin 
24  just  set3  any  small  value  (<  5  x  10  *®)to  zero  for  practical  purposes.  A 


line  26,  E(NORS)  is  calculated  by  using  the  formula  shown  previously. 
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c 

c 


e<  suo/m  coMpur  at  ion 


ISN 

0026 

S00  P*  0«  0*  00 

ISN 

0029 

00  30  I*l*IT0T 

ISN 

0030 

jUP*i  NO  PS  -  It  *NAPi  II  *1 

XSN 

0031 

00  60  J»l*dUP 

ISN 

0032 

1  XJ*K  IT  <  I  t+J 

XSN 

0033 

Z FI  IXJ«GT»60*  GO  TO  SO 

ISN 

0035 

60 

SDO  P*S00P  ♦  i«D  *00~PJ  IS  DN  <  I »  IX  J 1 

ISN 

0036 

50 

CONTINUE 

IS* 

0037 

SOQ  *SDOP 

ISN 

003d 

RETURN 

ISN 

0039 

ENO 

Line  35  gives  the  calculation  of  E(SDO) ,  which  is  also  formulated  in 
the  above  section.  At  lines  28  and  35,  SDOP  is  used  instead  of  SDO  to  avoid 
alternating  the  current  value  stored  in  SDO.  Line  37  then  updates  the  SDO 


value. 


Notes  on  Subroutine  OTPUT 

This  subroutine  prints  out  the  final  results  of  the  optimization.  It 
includes:  the  optimal  kit,  its  price  and  performance  indices  E(NORS)  and 
E(SDO);  parameters  set  b^'s;  weighting  factor  of  RATIO  computation;  kit  57-18 
and  its  relevant  values.  Number  of  iterations  is  also  printed. 


APPENDIX  B.2 

FORTRAN  LISTING  OF 
GREEDY  ALGORITHM 
COMPUTER  PROGRAM 
AND  SOLUTION  OUTPUT 
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GHEE DV  ALGORITHM 


ISN  0002 


DIMENSION  KIT  1500* 


ISN  0005 
ISN  0004 
ISN  00C5 
ISN  OOCO 
ISN  0007 


CALL  In  POT  It  K  IT  .PHI  CE  I 
CALL  TABLE 

CALL  INKITI  KIT. PRICE# 

CALL  OPTIMI  K1T.PR1CE.EN0RS.SC0* 
CALL  0TPUT1  I  KIT.M  IC E » EmO AS. SO  0* 

CALL  1NPUT2I  K  IT  .PRICE  I 
CALL  TABLE 

CALL  INKITI  KIT.PRICEl 

CALL  OPTINI  KIT.PRICE.cNORS.SCO* 

CALL  OT  POT  2  I  KIT  •  P4  IC  E»  E  NOFS.  S00» 


SUBROUTINE  INPUT  K  K  IT  (PRICE  I 
DIMENSION  KITISOOI.  OFT!  301 
COMMON  PUISON I SOO  *601 »BOt SC  I 

COMMON  /BLOCK  l/ITOT  (NITER  .MOPS.  MU  P.BU06 ET » SKF .  PAR 
COMMON  /BLOCK  2/  NAPIS«0*»CtSTI500».KMUI  3001 
COMMON  /BLOCK 4/  KTSTlSi SOO  I. .GOAL  .SGOAl  .BUDKT 
COMMON  /PHASE/  IPHASE 
COMMON  /INPUT /NDAYS <i*FP .  NSUO 
DOUBLE  PRECISION  PJISON.BO 

ORTA  KRU/*LNU«/tRR4/<  AMR* /  .OFT/ 7* 1. IS • • 53. 22*. *2/ 


ISN  0012 
ISN  Ot)  15 
ISN  OU 14 
ISN  00  IS 


IPHASE* 1 
N  IT  ER*0 
PRICE-0.0 
ITEM-0 


ISN  0010 


ISN  0017 


RE40IS.109I  PAR 


REA0I3.10I)  SFPD»  NJ  AT  S»  OF  F  *  MOP  S»  13U0.  S  KF 


ISN  001A 
ISN  0019 

ism  ooac 


HOP  ■MOPS*  1 
10  ITE  M» ITEM.l 

20  NEAOI  1* 102<END*S0I  CAP.NRC .NAP1  IT  EM* . TOR. CO ST  I  ITEM !• T »P 
ONL V  LRU  items  bill  BE  PNOCESSCO 


ISN  0021 
ISN  0023 
ISN  0024 
ISN  002S 


IF!  TTP.EB.XRU1  GO  T.0  SO 
GO  TO  20 
50  KMUlITEMIaO.O 
L  EN  T  ■  NS  UO 


ISN  0020 


IFICAP.EU.RRRI  LENT  »N  SUO  .NRC 
OEM  AN  0  AMO  KIT  VALUES 


ISN  0020 
ISN  0029 
ISN  005E 
ISN  OOSl 
ISN  0052 


00  40  IOAT»l»L*NT 

40  A  *AM  ITS  N*  ■  AMU!  I  IE  Ml  *f  LO  ATI  RAN  ITE  M  M  .TOR.  OF  Tl  IOATI 
lITUTIMMO 
GO  TO  10 
90  I T0T*ITEN«1 


6.2.2 


c 

C  COMPUTATION  OP  THE  REUU  IRENE NT  KIT  3718  4  TARGET  VALUES 

C 


ISN 

0033 

a  UO  OE  T>  0.  0 

ISN 

0034 

00  60  I-laITOT 

ISN 

0035 

KT57101  IlalFllllNW  II  •SIF*«SI 

ISN 

0036 

IP!  KTS7161  Il.LT.ll  XT  37  161  Hal 

ISN 

0036 

6U0  GE  T- BUOGET aCOSTt 1 1 «FLO* ll *T37l (1 1 1 1 

ISN 

OOSS 

c 

60 

CONTINUE 

ISN 

0040 

BUOKT  -BUDGET  - 

ISN 

0041 

6UO«CT-.9*BUOGET 

ISN 

0042 

c 

101 

P  ORNATl PSa2» IK* I2a IKa  PSa  2a  IX • 12 »T 2Sa 12* T26tF6a4» 

ISN 

0043 

102 

PON MAT! Tl9aA3*I2f 14 a!2»T32 jF S.4.T 46*P6a  2* A3  1 

ISN 

0044 

10S 

FORMAT!  141 

ISN 

004S 

109 

FORMAT!  P  10*  01 

ISN 

0046 

110 

FORMAT!  1006a  01 

ISN 

0047 

RET  CRN 

ISN 

00  46 

END 

ISN 

0002 

subroutine  otputk* it.pricc *enors <sooi 

ISN 

0003 

DIMENSION  KIT  ! S001*  SN!  4  1 

ISN 

0004 

COMMON  POISON!  S  00(6  01  »BOl  SO 

ISN 

0005 

COMMON  /BLOCK 1/ITOT aNI TEN# AOPSa MU F aBUDbET a  5 KF a  FAR 

ISN 

0006 

COMMON  /MLUCK2/  NAP! SOO I (CCS1! 500  ItXMU!  9001 

ISN 

0007 

COMMON  /BLOCK  4/  KT5  7161  S 00 1(  V60  AL  (SGO  AL  (BUD  KT 

ISN 

0006 

COMMON  /PHASE/  IP  HA  SE 

ISN 

0009 

V 

c 

OOUBLe  PRECISION  PJISJN(BO 

ISN 

0010 

V 

r 

ITEM-0 

ISN 

0011 

V 

MAI TE!6(200! 

ISN 

0012 

■Ml TE!6(2 10 1  ITOT 

C 

c 

PRINT  OUT  THE  TARGET  VALUES 

ISN 

0013 

BRI  T£( 6(205  i  !KT5T16!  1 1.1-1,  ITOT# 

ISN 

0014 

■  RS TE 16  (220  i  BUOKT  (VGOAL. 3GCAL 

ISN 

0015 

•  RITE!  6  (2301  BUDGET. 

ISN 

0016 

C 

■RITE!6( 100 1  PRICE. EnCRS. 310 

ISN 

0017 

•  RITE! 6. 1101  NITER. I  BO! J  1  a «-  I.NUP  1 

ISN 

0016 

•RI TEl 6(1091  PAR 

ISN 

0019 

NRI TEl 6( 1051 

ISN 

0U2C 

c 

•  RI  TEl  6(  120 1  !KIT111(  I-ldlOII 

ISN 

0021 

V 

RE-INO  1 

ISN 

0022 

t 

100 

FORMAT!  SX( 'PRICE!  XI  •» •. P 15.2.5X. • El N0RS/X1- • . P  lOaS (5X .'El SOO/Xl-' . 

»F  10  a  5  (/  /I 

ISN 

0023 

IOS 

FORMAT! //(5X(*1NIT1Al  KlT'.//» 

ISN 

0024 

109 

FORMAT!  //(SX(  'HEIGHT  «*.PlC.3.//l 

ISN 

0«|2S 

110 

FORMAT! /(  15 X( 'NUMBER  JF  IT  EK  ATION  S' a  14.  //  (  5  X(  •  «  J1  •  .  1C!  2 X(  F 10.  81  1 

ISN 

0026 

120 

FORMAT!  SX.25I4I 

ISN 

0027 

200 

FORMAT! // a 15X (*  GREE OY  ALGOFI TUN' . //. 1SX  PHASE  l'.///» 

ISN 

0026 

209 

FORMAT!  Sx  (2  SI  41 

ISN 

0029 

210 

FORMAT! //dSX  ('INITIAL  KIT  COMPUT  ATION*  .//.  SX(  14,  •  ITEMS  IN  COMPUTATION 

IT  AT  ION*  *//*  20Xa'KlT  57l8aW/l 

ISN 

0036 

220 

F0RNATI//.15X ('TARGET  S'.SX.'FERFO  FMANCE  OF  STIC  XlT' ( // . SX. 'COST 

1  OF  571S  XIT-4*(F15#2  (5X  .  •  FORSbOA  t«*  (FlOaS.SX.'SOOUOAL**  »F10.  5.//1 

ISN 

003 1 

230 

FORMAT!  SX.' BUDGET  UP  1NITIILKIT  IS  *  '.FlS.  2.//1 

ISN 

0032 

SO 

IPHASE-2 

ISN 

0033 

MET  URN 

ISN 

0034 

EMO 
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ISN 

0002 

r 

SUBROUTINE  INPUT24K IT. PRICE  1 

C 

c 

NC*0  IN  THE  PEN  T  INC  NT  OlTS  FFOF  THE  T4PES  AND  C  4*0  S 

ISM 

0003 

DIMENSION  KtTISOOl*  SNI  41. TM  41.  OP  SHI  301  .OHOSI  3(1  .RlIMT  (  30  l.RPHSI  30* 

ISN 

00(4 

CONPON  PUtSONOOO.Oei  .00150 

ISN 

OOOS 

CONNOR  /BLOCK  1/ ITOT .NITER. 10FS. NUt.BUOGCT • SKF • FAR 

ISN 

0006 

CONNOR  /BLOCK  2/  N  AP  <S  Oilt.C  OSTl  9  00  l.KNUI  9001 

ISN 

0007 

CON  NON  /BLOCK*/  KT5  710I  500I.V60AL  .SGOAL.BUOKT 

ISN 

OOOS 

CON  NON  /PHASE/  IPHASE 

ISN 

0009 

CONNON  /INPUT /NBA TS .ORP.NSLO 

ISN 

0010 

DOUBLE  PRECISION  POISON.BO 

ISN 

0011 

DATA  XRU/ •LRU'/.SRJ/* SRU* / .RRR/ *R RR*  / 

ISN 

0012 

OATA  OPGN/7-1.1S..S3.  22*. 5 2/ 

ISN 

0013 

V 

NtT  ER«0 

ISN 

00 1* 

PRICE«0.0 

ISN 

0015 

c 

ITENaO 

ISN 

0016 

V 

r 

REA 04  9.  1091  PAR 

ISN 

0017 

V 

r 

10 

ITEH>ITEN41 

ISN 

0010 

20 

REA 04  1.  111. END-2001  4  SN4I1.I-1.4I .CAP.  NRC > NA Pt ITEM  I . TOR, OCR. COST 

•« ITEN i»  Tip 

ISN 

0019 

c 

BRR— TOR— OOR 

V 

c 

r 

ONLV  LRU  t  SRU  ITENS  MILL  BE  PROCESSED 

ISN 

0020 

V 

IF4TVP.Eu.SRUI  GO  TO  30 

ISN 

0022 

IHTVP.NE.KRUI  GO  TO  20 

ISN 

002* 

INU-K1TI  ITEN1 

ISN 

0029 

GO  TO  17 

ISN 

0026 

30 

INU  — 0 

ISN 

0027 

17 

T  OT  flU-O. 0 

ISN 

0026 

DO  HO  10  AVa  l.NOA  VS 

ISN 

0029 

c 

110 

0N0S4  10AvlaFL0AT4  NAP4  1TENI  la  TOR  aOFGN4  ID  IV  1 

C 

e 

RRR  PROCESSING 

ISN 

0030 

IF4  CAP.  NE.RRR  1  GO  TiO  IS 

ISN 

0032 

XF4  TVP.EU.XflU.ANO.lNU.LT.il  iNUaj 

ISN 

003* 

K TO TaNSUU.NflC 

ISN 

0039 

1F4 KTOT.GT. 30  1  KTOT -3 0 

ISN 

0037 

l>0  120  lOAVa  1  .K  TOT 

ISN 

0030 

120 

RPRS4  I0AVI-0.  0 

ISN 

0039 

RPR HO L> 0. 0 

ISN 

00*0 

1F4  xTOT.uE.30l  GO  TO  |32 

ISN 

00*2 

K  TO  T*NS UU.l 

ISN 

0043 

00  12S  la  1.  KTOT 

ISN 

0044 

129 

RPRHOLaRPRHOL«FLOAT(NAP| IT EN M • BR k.DPGH 41 1 

ISN 

00*9 

RPR 34  NSUO+NROIMAPAmOL 

ISN 

0046 

jTOT-NSUD.NRC  *2 

ISN 

00*7 

00  130  I0AV«jT0T»N0AV3 

ISN 

0040 

130 

RPR  34  10  A V  »«Fl0AT4  NAP4  ITEN1  laCRRaO  FGN4  ID  AV-NAC  1 

ISN 

00*9 

132 

RON  T4  llaON0S4  11-RPRS4  |l 

ISN 

OOSO 

00  139  10 AV  a2 .NOA VS 

ISN 

0051 

139 

N0NT4 10 A V laRONT  4 IOA  V- 1 I.0NIS4I0AV I-RPRS4I0A VI 

ISN 

0092 

T  OT  Ru-ROHTI  t* 

ISN 

0083 

00  140  IUAV«2  .NOATS 

ISN 

0094 

XH  RONTIIOAVI.GT.  TjTRfll  TO  IROaRQN  Tt  TO  AV  1 

ISN 

0096 

160 

CONTINUE 

ISN 

0097 

GO  TO  TO 

B.  2.4 


I 


C 

c 

c 

86  PROCESS  INC 

ISN 

0099 

V 

19 

00  160  I0AT-1*N0AVS 

ISN 

0099 

160 

TOTRO*TOTRU«OMOS!  IOAV  I 

c 

0EM8N0  AND  KIT  VALUES 

ISN 

ooeo 

V 

70 

XNUllTCNI-TOTRS 

ISN 

0061 

KIT1ITENI-INU 

ISN 

0062 

PRICE-PRICE  ♦COST*  IT  EM  t  -FLOAT  »  INOI 

ISN 

0065 

r 

60  TO  10 

ISN 

00<« 

V 

200 

IT0T-1TEN-1 

V. 

c 

#» 

COMPUTATION  OF  THE  REQUIREMENT  KIT  S7l6  6  TARGET  VALUES 

ISN 

0069 

V 

BUOGET-O.  0 

ISN 

00(6 

00  90  I*1*IT0T 

ISN 

0067 

KTS7181  Il-IFIX!  AMU  II«SKF**5I 

ISN 

9068 

IF«  KTS7  1811  l.LT*  l»  KT  97  181  11-1 

ISN 

0070 

BUJGET-UUOGET  ♦COST!  II  4FL0A1I  kTSTiHIII 

ISN 

0071 

SO 

CONTINUE 

ISN 

0072 

f* 

BUOKT-ttUOGET 

ISN 

0073 

102 

F0HMAT1  SA *4A4*SA*  1 4  1 

ISN 

0074 

109 

FORMAT!  14  1 

ISN 

0079 

108 

FORMAT!  1006*01 

ISN 

0076 

109 

FORMAT!  F 10*  01 

ISN 

0077 

111 

FORMAT!  4A4*Tl9.A3tlX*  I1*1X*I2*T32  *FS*4*  IX  *FS*  4  *T46  *F8  •  2*  A  31 

ISN 

0078 

RETURN 

ISN 

0079 

ENO 

ISN 

0002 

SUBROUTINE  OT  PUT21X IT  *PRICf*ENORS *S00 1 

V 

c 

THIS  SUBROUTINE  PRINT  THE  FIFAL  RCSlLTS  OF  THE  OPTIMIZATION 

ISN 

0003 

OIMENSION  K IT !  9001 

ISN 

0004 

COMMON  POISONlSOOtoOl *8019(1 

ISN 

0009 

COMMON  /BLOCK t/ITOT *NI TEN* (OPS*  MU  F* BUDGET  *  SKF • FAR 

ISN 

0006 

COMMON  /BLOCK  2/  N  AP  !  9  OtflaCOST  (9  00  I*KNU1  9001 

ISN 

0007 

c 

COMMON  /BLOCK  4/  K  T9  71 8!  900 1  •  *60  AL  *SGO  AL  *BUD  KT 

ISN 

0008 

t 

double  precision  PJISONtBO 

ISN 

0009 

V 

MRITE!6*99I 

V 

c 

PRINT  OUT  THE  TARGET  VALUES 

ISN 

0010 

C 

mRI T£!6  *9  50  1  VGOAL* Sb QAL • BLOAT 

ISN 

0011 

MRI  TE16*990I  1XT5718!  Il.I-l.ITOT* 

ISN 

0012 

L 

•  RITE  16*90#  I TOT* NITER 

ISN 

0013 

nRI TE! 6 • 10U 1  ENORS*  SO 0* PR ICC 

ISN 

0014 

•  RITE!6*20U#  1 K  IT!  1 1*  I-  1  *  I  TOT  # 

ISN 

0019 

•RI TE16*  3001  (BO!  Jl  »  J-l.MUFI 

ISN 

0016 

/• 

•R1TE!6*S99I  PAR 

ISN 

0017 

90 

FORMAT! //*9X**T0TAc  NUMBER  OF  ITEEMS  UNDER  CONSIDERATIONS*  *  I* . // . 

♦SK**NUMBER  OF  ITERATIONS**  14* ///l 

ISN 

0016 

99 

FORMAT! //•  15X**CRCE OT  ALG08 ITHM  •* /*  19 A*  *PHA  SE  <•*//! 

ISN 

0019 

100 

FOR NAT! SA* • E!  NORS/A 1- *  *  F 10 *9 »SA •* El  900/ XI «• *F 1<* 5* 5X» ‘PRICE! A ••»• • 

•F  19*9*/// • 19X  »*OPTI RA  L  KIT»*//I 

ISN 

0020 

200 

FORMAT!  5A*25I4I 

ISN 

0021 

300 

FORMA  T!//»5A.*6*./.!lO!  2X»FlO*OII  I 

ISN 

0022 

999 

FORMAT! //•  9 X* ’HEIGHT  «**Flf*3*//l 

ISN 

0023 

990 

FORMAT! //*5X*  *TARG6TS**//5>»*  NORSGOAL  «* »F  10* !» 9X * *SO0GOAL » • *Fl 0* 

•9*9X«*PRICEm»**FiS*  2.//I 

ISN 

0024 

990 

FONNAT!  Sft*  291  41 

I*H_£029 

RETURN 

ISN 

0026 

ENU 

B.2.5 


ISM  0002 


SUBROUTINE  TABLE 


c  computation  of  the  cjmulativi  sums  cf  poisson  distributor 
c 


ISN 

0003 

CONPON  POISONlSOO.oQI .ttOISC! 

ISN 

0004 

COMMON  /BLOCK  1/ ITOT (NI TER » 10FS, MU FaBUOG ET. SKF . PA A 

ISN 

0003 

c 

COMMON  /BLOCKS/  NAP!3C0!(C0ST<  5  00*.  AMU  300! 

ISN 

0006 

V 

c 

DOUBLE  PRECISION  PJ  IS  ONaTE  PM.  —  -,  SUM.  TE  RL 

ISN 

0007 

V. 

TENL-Sa  D-  12 

ISN 

0006 

00  10  I- la ITOT 

ISN 

0009 

PNN-OBLEl  AMU!  I!  ! 

ISN 

0010 

TERP-CE  API-MNNl 

ISN 

0011 

P01 SO Nila  11-TERM 

ISN 

0012 

DO  20  0-2(60 

ISN 

0013 

1F<  TERM.LT.  TERL.ANP  .PiUSONI  I  (J- 11  .GT.9.  0-01 1  GC  TO 

ISN 

0013 

TERP-TERM-NNN/0FLO4  T<  j-1! 

ISN 

0016 

POISON!  I.  JI-POISON!  I.  J-llalERM 

ISN 

0017 

20 

CONTINUE 

ISN 

0010 

GO  TO  10 

ISN 

0019 

30 

DO  40  UK-U(60 

ISN 

0020 

40 

POISON!  ta  JA  1-  1*0 a OJ 

ISN 

0021 

10 

CONTINUE 

ISN 

0022 

RETURN 

ISN 

0023 

ENO 

ISM  0002 


ISN  0003 
ISM  00C4 
ISM  0003 
ISN  0006 
ISM  0007 

ISN  0000 

ISN  0009 
ISN  0010 
ISN  0011 
ISN  0012 
ISN  0014 
ISN  0016 
ISN  0017 
ISN  0011 
ISN  0019 
(SN  0020 

ISN  0021 
ISN  0022 


SUBROUTINE  INF  IT! KIT  a  PR  ICE  I 
C 

c  increase  ait  values  to  ensure  that  the  probability  of  any  item 

C  -Ml  CM  MILL  HAVE  MAlF  of  the  fleet  DO-N  is  less  than  so* 

c 

DIMENSION  KIT 49001 

COMMON  POISON!  900(601  »B0€  SO 

COMMON  /BLOCK  1/ITOT .NITER. FOPS. MU F ( BUDGET » S KF . PA R 
COMMON  /BL0CK2/  NAP  1300  1.CC5T1500  I.KMUt  300 » 

COMMON  /BLOCK  3/  K  IT  IN  1  3  00 1  (PRC  I N 
C 

DOUBLE  PRECISION  PlISONaBO 
C 

NCPS-lNOPSail/2 
00  10  ITEM*  i(  ITOT 
20  NUMaNEPSaNAPt  ITEM  la  K1  Tl  ITEM  I 
IF!  MUMtGT (601  GO  T4  1C 
IF! POISON! ITEM, NUMl,uT..S»  GO  TO  10 
KlTl  I  TEN! -KIT!  ITEM!  n 
PRICE-PRICE  aCOSTI  ITEM! 

GO  TO  20 

10  kit  INIITEMI-KITI  ITEM! 

PRC  IN-PRICE 

c 

RETURN 

ENO 


B.  2.6 


ISM  0002 


ISM  0003 
ISM  0004 
ISM  OOOS 
ISM  00C6 

ism  ooor 

ISM  0000 

ISM  0009 
ISM  0010 

ISM  0011 


ISN  0012 
ISN  0013 
ISN  0014 
ISN  0019 
ISN  0016 


ISN  0017 
ISN  0016 


ISN  0019 
ISN  0020 
ISN  002 1 
ISN  0022 
ISN  0023 
ISN  0024 
ISN  0026 
ISN  0027 
ISN  0029 
ISN  0030 
ISN  0031 
ISN  0032 


ISN  0033 
ISN  003S 


SUOROUTINE  OPTINl  Ki  T. PRICE  .EkORS.  SOOt 
C 

C  NS  IN  OPT  IN  II AT IJ  L  ROUTINE 

C 

OtMENSION  KIT4S00*.  SI  SOI 
COMMON  P01S0N4S00*6  01  •  AM  5(1 

COMMON  /BLOCK  t/ITOf .NITER. fcOPSt MUF. BUDGET • SKF * PAR 
COMMON  /BLOCK 2/  MAPI  9  C0».  CCSTl  900  t.XMUl  800* 

COMMON  /BLOCK  4/  KTSTlCl 500  1.460 AL  .SGOAL .BUOKT 
COMMON  /PHASE/  IPH1SE 
C 

DOUBLE  PRECISION  PJ ISON. BO .SLMSS. EPSO.P MO.E PS . S 
COMMON  /Sl/S 
C 

OAT  A  EPS. EPSU/1.D-J7.  S.Q-i  1/ 

C 

C  COMPUTE  THE  PERFORMANCE  CF  Mil  57-1 «  ANO  INITIATE  PARAMETERS 

C  0 

C 

CALL  CALC4KT571B. VuOAu.SGOALI 
00  1  J*  1.  MUP 
1  B  01 u»*S4 J I 

MRI  TE (6*300  I  480«  J«  .JM.NUM 

300  FORMAT*  •  i*/T9.* INIT IAl  VALUES  OF  PARAMETERS  B  •  /  I  1QI  2A.F  10.  81  *  I 

C 

C  START  OPTIMIZATION 

C 

C  STEP  1 

C 

ioo  call  marganikit. price i 
c 

NITER*NITER*1 

C 

C  STEP  2 

C  UPOATE  THE  PARAMETERS  B 

c 

SUM 99*0.0*00 
00  10  J»1.MUP 
PRO  *  1.0  *00 
00  20  I *1. ITOT 
KITPJ«K1TII  i«IU-lt*NAPI  II*  1 
IFI AITPu.GT.60l  GO  TO  20 
PRO*PRO*POISaN(  I.KITPJ1 
IFI  PRO.LT.EPSO*  GO  TO  15 
20  CONTINUE 

15  5 UM  Su*SUMSU*4 4  BOt 01  — P RO  t •• 2  I 
B 04  U»*PRO 
10  CONTINUE 
C 

C  STEP  3 

C  COMPARE  THE  NEM  VALUES  CF  0  MIlH  PREVIOUS  VALUES 

C  TERMINATE  IF  THET  ARE  SUFFICIENTLT  CLOSE 

C 

IFI  SUMSO.LT. EPS*  GO  TC  30 
GO  TO  100 
C 

30  CALL  CAlC4MIT.EN0RS.SC0* 

RETURN 

EMO 


ISN  0036 
ISN  0037 
ISN  0030 


ZSN  0002 


c 

c 

c 


SUBROUTINE  MAR6AN(KIT<PRICEI 


THIS  SUBROUTINE  PERFORMS  THE  MARGINAL  ANALYSIS 


I 


ISN  0003 
ISN  0004 
ISN  0009 
ISN  0000 
ISN  0007 

ISN  0000 


DIMENSION  RATIOI  SOO  I*  KlTlSIO  I 
COMMON  POISON! 500*0  01 *BQ4 SCI 

COMMON  /BLOCK  t/ITOT  »NI  TER»  AOF  3*  MUP.BUOG  ET  •  SKF  •  PAR 
COMMON  /BLOCK 2/  NAPI5C0I  *CCSTlSOO  ItXMUl  900» 

COMMON  /BLOCKS/  KIT  IN  (  900  »  *PAClN 
C 

DOUBLE  PRECISION  PJ  IS  ON  t  R  A  TIC  «B0*  SUM.  MM  tRMA  X 
C 
C 

C  START  mITH  THE  INITIAL  «IT  ANO  CALCULATE  THE  RATIO  OF  RETURN 

C  PER  DOLLAR  INVESTEO  FOR  EACH  ITEM 

C 


ISN 

0009 

PRICE -PRC  IN 

ISN 

0010 

00  10  I  *  1 • ITOT 

ISN 

0011 

KlTl  I  l-KITINI  II 

ISN 

0012 

NMAX1*KIT!I  MINORS*  1MNAP1  ll«2 

ISN 

0013 

IN  NNAXl*e£*60t  NM4 X  1  *6 0 

ISN 

0015 

IXS*KITIII+1 

ISN 

0014 

sumsoblE(pari«i  po  is  on  t  i  *nm  ax  i  i*  po  isoni  i  *ixs  it 

ISN 

0017 

DO  20  J* 1 *MUP 

ISN 

ooia 

I  XJ*KITt  1 1  +4  J*  1 1  •  N4  P(  11  +  1 

ISN 

0019 

IF4  IXJ*  GE.60I  GO  TO  2  I 

ISN 

0021 

NB-PO ISONI I«IXJ+ 11/ POISONt l» IXJ 1 

ISN 

0022 

GO  TO  20 

ISN 

0023 

21 

MRS  1*0+00 

ISN 

0024 

2o 

SUM *SUM+B04  Jl  +0L0G4  MM  1 

ISN 

0029 

10 

RAT  1011  l-SUM/OBLEIC  0ST1  HI 

V 

c 

r 

FIND  THE  MAXIMUM  RATIO 

ISN 

0026 

V 

100 

DO  25  I*  l*  ITOT 

ISN 

0027 

IF!  IcEM.ll  GO  TO  3j 

ISN 

0029 

IFI  RATIO!  II  .LT.RMAx  1  GO  TO  2S 

ISN 

00  31 

30 

IMA  X*I 

ISN 

0032 

RMAXmRATIOI  II 

ISN 

0033 

25 

CONTINUE 

ISN 

0034 

XPR  1CE-PR ICE* COST  II MA X 1 

V 

c 

TERMINATE  IF  THE  TARGET!  ARC  ACHIEVES 

ISN 

0035 

t 

IFI  xPRICE *GT* SUOGET 1  GO  TO  SC 

c 

f 

IF  NOT*  UPDATE  THE  KIT  VALUE  A  AO  THE  RATIO 

ISN 

0037 

V. 

KZTIIMAXImKITIIMAU  +1 

ISN 

0034 

PRICE *X PR  ICE 

ISN 

0039 

NMA  Xl>Kiri  IMAKI+I  NJ  PS  -  1 1  «N  API  IM  AX  1  +  2 

ISN 

00  40 

IFI NMAX1.GT*60I  N MA  XI » 60 

ISN 

0042 

IXS*K  IT  UMAX  1+1 

ISN 

0043 

SUM  >OBL  Cl  PAR  1  *4  PO  IS  ON  1 1  MA  X  #NP  AX  1 1  «POI  SO  M  IM  AX 

ISN 

0044 

00  40  J* l»MUP 

ISN 

0049 

IXJ+KITI  IMAXl+l  II  +NAP4  IMAXI  +  1 

ISN 

0044 

IFI  IXJ*  GE*60I  GO  TU  49 

ISN 

0044 

MR* POISON! IMAK»ZXJ+ll/POISCNlZMAX (IXJ 1 

ISN 

0049 

40 

SUM*SDM+aoi  Jl  «OLOGI  mm  I 

ISN 

0050 

45 

RATIOI  I  MX  xl  ■  S  UM/OBc  El  CUSTI  INAxI  1 

ISN 

0091 

GO  TO  100 

ISN 

0092 

90 

RATIOI IMA XI *0*0 +00 

ISN  0093 
ISN  0099 
ISN  0094 


C 


IFI  RMAX.GT.  0*  D+00*  GO  To  ICO 

RETURN 

END 


ISM  0002 


C 

c 

c 


subroutine  calci kit •eaors* *dc » 


THIS  SUBROUTINE  COMPUTES  THE  EINORS/Xl  A  NO  E(SCO/X» 


ISM 

0005 

DIMENSION  SOOIfKtT  (9001 

ISM 

0004 

COMMON  PO I SONl 3 00 »o 0 •  >80(3(1 

ISM 

0005 

COMMON  /BLOCK  1/ITOr .HI  TEN » AOPS* MU  I 

ISM 

0006 

COMMON  /BLOCK  2/  NAP  (  5  (  0  ItCCSK  900 

ISN 

OOOT 

COMMON  /Sl/S 

ISM 

0000 

C 

DOUBLE  PRECISION  PJ  IS  ON  t  'j»  SUM  .X  NO 

ISN 

0009 

C 

r 

EPS>9*D-10 

V 

C 

f 

E( NORS/K 1  COMPUTATION 

ISN 

001C 

V 

S  UM  *0*0  +  00 

ISN 

0011 

S  (HOPS*  ll-i.O  0 

ISN 

0012 

00  10  K«1(N0PS 

ISN 

0U1S 

S(K  I*  1*  0*00 

ISN 

001* 

00  20  I*ltITOT 

ISN 

001S 

1XJ-K1T(1  l«(fc-li*N4P(  11*1 

ISN 

0016 

I P(  IX J» GT *6 0 i  CO  T.)  20 

ISN 

0010 

S(KI*S(  KI*POISOM  I.  IX  Jl 

ISN 

0019 

IF(  S(  Kl.LT.  EPS1  CO  TO  15 

ISN 

0021 

20 

CONTINUE 

ISN 

0022 

SUM*SUM«S(KI 

ISN 

0023 

CO  TO  10 

ISN 

0024 

IS 

six  >*a»D»oo 

ISN 

0029 

10 

continue 

ISN 

0026 

XN0PS*DFL0AT| NOPSI«S( AOPSt  11-Sum 

ISN 

002T 

c 

END  RS*X  NORS 

c 

f 

E(  SO  0/  X 1  COMPUTATION 

ISN 

0020 

w 

SOOP-O.  0*00 

ISN 

0029 

DO  90  1*1 • ITU T 

ISN 

0030 

j UP  *(  NO  PS  -  1 1  •  N A  P(  II  *1 

ISN 

0031 

00  60  J-1*JUP 

ISN 

0032 

1  XJ*K  IT  1 1  »*J 

ISN 

0033 

IF(  IXJ.CT .60 i  60  TJ  SO 

ISN 

0035 

60 

SOQ P*500P Xl-D ♦OO-PJ IS3N( It  XXJI 

ISN 

0036 

50 

CONTINUE 

ISN 

0037 

SOO«SOOP 

ISN 

0030 

RETURN 

ISN 

0039 

E  NJ 
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APPENDIX  C 

C.l  Notes  on  Conversion  of  the  D029  Computer  Program  to  an  IBM  Compiler 

C.2  FORTRAN  Listing  of  the  IBM  Version  of  the  D029  Algorithm  with  a 
Solution  Output 


APPENDIX  C.l 


NOTES  ON  THE  CONVERSION  OF  THE  CURRENT  D029  COMPUTER  PROGRAM 

TO  AN  IBM  COMPILER 

IBM  and  CDC  FORTRAN  compilers  are  not  exactly  the  same.  Any 
computer  code  compatible  with  a  CDC  FORTRAN  compiler  is  not  compatible 
with  an  IBM  FORTRAN  compiler  and  vice  versa.  Because  of  this,  the 
D029-CDC  code  has  been  converted  into  an  IBM  compatible  code  in  order 
that  appropriate  comparisons  could  be  made  between  the  current  D029 
program  and  the  Greedy  Algorithm  program. 

CDC  FORTRAN  compilers  accept  variable  names  of  seven  characters, 
but  an  IBM  FORTRAN  compiler  can  accept  variable  names  of  only  six 
characters.  Therefore,  variable  names  with  seven  characters  were 
truncated  to  six  characters.  In  addition,  the  following  statements  are 
not  accepted  by  an  IBM  FORTRAN  compiler. 

IF (IMAIN . EQ . 3HRRR) .  .  . 

and  READ(1 , 107) .  .  . 

IF(E0F(1).NE.0)  GO  TO  200 

The  first  statement  is  converted  by  defining  an  alphanumeric 
variable  RRR  in  the  data  statement  as  follows. 

DATA  RRR/ 'RRR 7 

IF ( IMAIN. EQ. RRR).  .  . 

The  second  statement  checks  for  the  end  of  file  in  the  data  set 
from  a  tape.  The  conversion  is  accomplished  by  combining  the  two  lines 
together  as 

READ (1,107, END-200).  .  . 
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The  word  size  of  a  CDC  computer  is  64  bits,  but  on  an  IBM  370,  a 
word  size  of  32  bits  is  utilized.  This  64  bit  word  size  causes  the 
memory  space  on  a  CDC  to  be  expensive  and  memory  saving  measures  are 
frequently  devised.  Consequently,  special  subroutines  for  word  packing 
and  unpacking  are  included  in  the  current  D029  program. 

It  is  mentioned  on  page  11  of  the  Preliminary  Evaluation  of  D029 
WRSK  Model  [8]  that  the  QPA  for  each  item  is  packed  into  six  bits  or 
1/10  of  a  word.  The  cumulative  probability  of  demand  for  an  item  is 
packed  into  20  bits,  or  1/3  of  a  word.  The  probability  values  for 
quantities  K,  K+QPA,  K+2*QPA,  ...»  K+UE*QPA  are  stored  when  QPA  is 
greater  than  one.  When  a  probability  value  for  a  quantity  between 
K+j*QPA  and  K+(j+l)*QPA  is  needed,  this  probability  is  estimated  by 
interpolation.  On  the  IBM  370/168  at  Oklahoma  State  University,  sufficient 
memory  is  available  such  that  word  packing  and  unpacking  is  not  needed. 
However,  double  precision  is  employed  to  increase  the  accuracy  in 
calculating  the  cumulative  probability  values  for  E(N0RS/X)  and 
E(SDO/X). 

The  IBM  version  of  the  D029  computer  program  is  summarized  as 
follows.  The  input  section  of  the  current  D029  has  been  adapted  as 
before  with  some  necessary  modifications  due  to  machine  differences. 

The  optimization  routine  has  not  been  changed.  Subroutine  ADD  has  also 
been  preserved  with  some  modifications.  Subroutine  TABLE  of  the  IBM 
version  calculates  the  cumulative  Poisson  sums  and  stores  the  values  in 
the  array  Poisson  (500,60).  Poisson  sums  are  retrieved  from  this  array 
when  required.  The  FORTRAN  listing  of  the  IBM  version  of  the  current 
D029  program  is  provided  in  Appendix  C.2. 


C.  1.3 


APPENDIX  C.2 


FORTRAN  LISTING  OF  THE  IBM  VERSION  OF  THE  D029  ALGORITHM 


WITH  A  SOLUTION  OUTPUT 


ISN  0002 
ISN  0003 
ISN  000* 
IS*  0009 


IS*  0006 
ZSN  0007 
ISM  0006 
XSN  0009 
IS*  0010 
ISN  0011 
ISN  0012 
ISN  001* 
ISN  OJU 


c  ••••*••*••***•*• «• •••*••*•«• •• *•••«• •••• 

c  *  • 

C  •  0029  ALGORITHM  * 

C  •  U  S  UK  FORCE  * 

c  *  « 

c  *  • 

c 

DIMENSION  KIT  IS 00 1 
00UB1.E  PRECISION  El  IM  it  500 
COMMON  /0UO/ISUH«1P*SS 
COMMON  /PRINT /MINPUT* MxMu.PlTEA 
C 

C  DETAILED  PRINT  OUT  OPTION.  ALL  DEFAULT  TO  ZERO. 

C 

C  MINPuTbI  MILL  print  out  input  oat  a 

C  HXMt»l  mill  PRINT  JUT  AVERAGE  DEMAND  AM  tM  I »  ,  AND  NFPOI  it  FOR  EACH 

C  ITCH 

C  MITER»1  MILL  PRINT  OUT  CPU  TIME  AND  ITEMS  ADOEO  FOR  EACH  ITERATION 

C 

C  «•* 

call  time 

i  call  input ii  price# 
call  table 

CALL  INKITI  AIT.  PRICE! 

CALL  0PT1M1  KIT.ENOR S.  SOJiPfilCEl 
CALL  OTPUTI  KlTiENORS.  SOO.PFICE* 

IF*  IP  ASS.  CO.  2  I  GO  TO  I 
STOP 

in 


ISN  0002 


ISN  0003 
ISN  000* 
ISN  OOOS 
IS*  0006 
ISN  0007 
ISN  OUOB 
ISN  0U09 


ISN  0010 


BLOCK  OAT  A 

C 

C  SLOCK  OATA  SUBPROGRAM  TO  INITIALIZE  VARIABLES  1*  LABELED  COMMON 

C 

DIMENSION  BINTE'R*  101 

COMMON  /ULK6/B1 NTEK/dLKO/ISUN .IPASS 
COMMON  /NLKS/  VPT.S0OT.  VPGOAL.SOOGCL 
COMMON  /C0UNT/IC3U* T.  lADDi  2001 
COMMON  /PRINT  /HINPUT.  MxMu, PITER 
OOUBLE  PRECISION  VPb JAL. SCOGOL 

OATA  BINTER/S  *3  E6 .  3.01b*  ♦ . 5 F6. 6. 0 Eo .  7. SEN • 9. 0E6 . 1 0.0E6 . 

♦  3»0.0/»ti*  A$  S/  l/»  *PT»S  DJ  T/2*  0./.  ISOM/  (/.IC1UMT  /  0/  • 

.  MINPUT  .MXMU.MITi  R/  J*0/ 

END 


ISN  0002 


ISN  0003 
IS*  000* 
ISN  OOOS 
ISN  O0Oo 


SUBROUTINE  CAlClKITi.vP.  SOO  I 
C 

C  THIS  SUBROUTINE  COMPUTES  THE  E  INORS/ XI  ANO  EISDJ/XI 

C 

01MENSI0N  SI*  SOI  (KIT*  3001 
COMMON  P3IS0NI 300* 1 401 
COMMON  /0LA  0/  ISOM*  IP  A  IS 

COMMON  /BLK 1/ ITOT. R OP. B006ET/BL K2 /XMU*  SOO I.NAP1 9001. COST* 9001 .NFP0 
♦  «  SO  0  I  /BLK 3/  NO  PS  •  SKF  •  I  UC  EFF  .IKK 
C0MMON  /BLX7/S1 

OOUBLE  PRECISION  PJ  IS  JN.  S  1  .VP. SOO  .EPS 
OATA  EPS/  1.0-  12/ 


ISN  0007 
ISN  OOO0 
ISN  0009 


c 


p 

1 


f . 

f  _■ 

f 

►  - 


X 

*  " 

>  ' 

i 


3 

■ 


ISN  0010 
ISN  0011 
ISN  0012 
ISN  0013 
ISN  0014 
ISN  031S 
ISN  0016 
ISN  0010 
ISN  0010 

ISN  0020 
ISN  902i 
ISN  0022 
ISN  0023 
ISN  0024 
ISN  002S 
ISN  3026 
ISN  0027 
ISN  0026 
ISN  0030 
ISN  0031 
ISN  9033 
ISN  0034 
ISN  0036 
ISN  0037 

ISN  0030 
ISN  0030 
ISN  0040 
ISN  304J 
ISN  0042 
ISN  0043 
ISN  0044 
ISN  0046 


c 


c 


11 

10 


S  11  NUN*UO«0O 
S  DO  a  0*3  *00 
00  10  I>1*1T0T 
JU7»I  NOPS-llaNAPI  II  41 
00  11  OaltJUf* 
ujHiri  i  i4j 

IFI  IXJ.GE.NFPM  III  SO  TO  1C 
SOU  -SOO 41*0400- POISON!  ItlXol 
CONTINUE 


•  Pa  0*  04  00 
IUa IUEEFF 
ittl«UK 
00  2S  N 2« I* IU 
NaJKKR.02 
S  11  NI«i,i>400 
00  20  I «1« ITO T 

I  XJ*X  IT1 1 141  N“1  1*N*  P4  1141 
1«  1*0.  GE.NFPOI  III  GO  TO  2C 
S  11  N I «S  UNI  aPOISONl  I*  ISO* 

IP!  SMN1.LE.EPS  »  GJ  r  J  22 

20  CONTINUE 

21  VPa  VP4FL0ATIN1*ISU  N411-S11NM 
IF«  SUN  l.£w.  0*0400*  SJ  TO  30 
GO  TJ  25 


22  S  11  NI-0.00 
GO  TO  2l 
25  CONTINUE 
RETURN 

30  00  35  N2«1*N 
36  5  11  N2  1*0.  0400 
RETURN 
ENO 


ISN  0002 


ISN  000 J 
ISN  0004 
ISN  0005 
ISN  0007 
ISN  0000 
ISN  0009 
ISN  0010 
ISN  0011 
ISN  0012 

ISN  0013 
ISN  0014 
ISN  0016 
ISN  0016 
ISN  0017 


SUd ROUTINE  tine 

c 

C  T  41 S  SJBROUTINE  GX0E3  CPU  TIMES 

C 

C  ELAPSEIII  IS  an  ASSENdLV  LANGUAGE  SUBROUTINE  anICH  TIMES  THE 

c  EXECUTION  OF  a  PORTION  OF  THE  CALLING  PROGRAM 

C 

COMMON  /0LK9/ISUM.1PASS 
0  AT  A  K/  0/ 

1 Fi  k.Eu.OI  GU  TO  30 
CALL  ELAPSE  1  ICPOI 
ISiOM*  is  UN  4ICPU 
T CP U«  FLOAT!  lCPUI/lOJO* 

TSUM» FLOAT!  IS UM 1/10  00. 

■  R1 TE16  >101  TSUM.TC Pu 

10  FORMAT  « IX* 101 2X.*:pu  a*.  FT .3. •  SEC.*  »SX • • INC  RE AMENT  a*, 
4  F7.3.2X*  1(H  *  *•  II 

RETURN 

30  CALL  EL  APSE!  ICPUt 
X«l 

RETURN 

ENO 


i 


i 

i 


ZSN  0002  SUM  ROUT  INC  INPUT! (  PRICE! 

c 

C  REAO  IN  THE  PERTINENr,  )(U  FROM  THE  TAPES  END  CAkOS 
C 


ISN 

0003 

DIMENSION  OPGNI  SOI*  0M0S1  30  I.RPRSI  301  «RMNT<  3  01  •  NSNt  4! 

ISN 

0004 

COM NON  /8LK 1/ XTOT *4UP .ttUDGET /DLK2 /XMul  SOO 1 *  NRPIS 00 1 *COSTf  5001 »MPPO 

41  9il0l/BLK3/N0PS*SKF*t  JEEFF.IXK 

ISN 

OOOS 

COMMON  /HLK  0/  ISUM * IPR  SS 

ISN 

0006 

COMMON  /HLK  12/ITOT l *1  TOT 2 

ISN 

0007 

COMMON  /PRINT /NINP3  T» 4ANU*  PITER 

ISN 

00<M 

INTEGER  OArS.UPAI.J  1»*:  .UPA 

ISN 

0009 

INTEGER  XRU.RRR.S4U*  END 

ISN 

0010 

ORTR  XRU/*  LRU*/*RRR  /*i4RR*  /  *OPGM/  7*1.  13  .  .53*  224.52/ 

ISN 

0011 

r 

ORTR  SRU/ «SRU  •/•£««/•  ENO»/ 

ISN 

0012 

V. 

1 

FORMRTI  •  !•///  101  I*  2R*  4PHRSE  *»H**  Z  ALCULAf  IONS*.  2 X*  101  *4*  II 

ISN 

0013 

101 

FORMRT  1  F5.2*1K.  12  •  1X*F  •>.  2*  IK.  12*  T  25  *1  2*  T2  8  *F  6.  4  1 

ISN 

0014 

102 

FORMAT  1  IX*  F5 .2  « 1  X*  12  .  IX. F  S*  2.  IX.  12*725.12.  T2«»F6*4I 

ISN 

001S 

107 

FORMRTI RR  *»  T1 9*R3 *1 2* IX » 12 «T32* FS *4* IK* FS. * • T 4 6* F0*2 *R3 1 

ISN 

0016 

c 

177 

FORMRTI  IS*2X*4R4*2R*R3*12*  IX  *  12 .2  X*FO.  4  .  IX*  FO*  4*  2R.F0*  2*  2X»  R3  1 

ISN 

0017 

V. 

■  RITE  16*11  IPASS 

ISN 

00  Id 

100 

IF!  IP  ASS*  Em*  1 1 

♦  HER  M  5*1011  PGM I (OR  TS  *PGN*  NOP S*  01  *SxF 

ISN 

0020 

I  FI  1P4S  S*£ll«  1 1 

♦  ■RI TEf 6*102 IPGMI* OR TS .PGM. NOPS*  01 *SxF 

ISN 

0022 

MUP*NQPS*1 

ISN 

0023 

IXXaMUP 

ISN 

0024 

IUE£FF*NOPS 

ISN 

0075 

1TM«0 

ISN 

0026 

MM*  0 

ISN 

0027 

IF  1 0 1.  GT  .OAT  SI  D1-D47S 

ISN 

0029 

PGM  1*0* 

ISN 

0030 

0  0  0686  11*1*01 

ISN 

0031 

5606 

PGM  1*PGMI«0PGM<  JJ  1 

ISN 

0032 

IF  1  IPASS.EQ*  ll  PGM*PGMI 

ISN 

0034 

r 

103 

ITM»ITM4l 

ISN 

0033 

V 

105 

REAO  1  1* 107*ENO>20JI  4SN. IPR IN. RC  .MPRI* TOR. OOR  .C0ST1 ITMI .LRU 

ISN 

0036 

IF  (IPASS. Ed. 21  GO  TO  32l 

ISN 

0036 

c 

MM*  MM41 

ISN 

0039 

V 

I FI  MINPUI.EU*  11 

♦■RiTE«6*l77INM*NSN* IN41N*  AC *I1P  AI  *TDk*D OR  *C 0$T  1 1 TM  I*  tRU 

ISN 

00*1 

321 

NRPI1TNI  a  UPRX 

ISN 

00  42 

r 

BRRaTOR-MOR 

w 

c 

r 

ONLT  LRU  6  SRU  ITEMS  MI  Lt  HE  PROCESSED 

ISN 

00*3 

\» 

I  F<  LRU*LT.XRUI  GO  TO  103 

ISN 

00  #5 

IF  (LRU*EO*XRU!  GO  TO  100 

ISN 

0  J*7 

IFI  LR0.2U.5RU  I  GO  ftO  (  105.  1061*  IP  ASS 

ISN 

00*9 

GO  TO  10S 

ISN 

OJf  0 

106 

TOT  RU*3* 

IS  4 

oasi 

Il«Ol 

ISN 

00  SO 

IF  1  IPASS.E 0*  21  GO  TO  109 

ISN 

0054 

IFI  XMRIN*EM*RRNII14,>1»RC 

ISN 

00  So 

109 

oo  no  1*1*11 

ISN 

0037 

110 

ISN 

0054 

IF  I IPASS.EU* tl  GO  TO  1*9 

ISN 

0060 

1*0 »♦! 

ISN 

0061 

00  US  1*3* OR ys 

ISN 

0062 

IIS 

OMO  S<  I*«MPRI*TOR*OPG<M  I! 

ISN 

0063 

IF  (  IMR  IN*NE«  RRRI  i»  0  TO  19! 

c 

C  LRU  RRR  PROCESSING 

C 


C.  2.4 


ISN 

00  65 

mii  i«nc 

ISA 

0060 

IF  IK.GT.30I  A*  30 

ISA 

0060 

00  120  l«l.K 

ISN 

0069 

120 

RPR  SI  11*0. 

ISN 

0070 

RPRHO-O* 

ISN 

0071 

1F1  K.  GS  .  3 01  GO  TO  132 

ISN 

0073 

A*0  14  1 

ISN 

0074 

00  125  1*1.  K 

ISN 

007S 

126 

RPR  HO  «  RPRM04  FLO*  Ti OPAI 1 *f)RR«  OP  GNf 1 1 

ISN 

0076 

RPRSi01*4C*ll*  RPR10 

ISN 

0077 

j *J  ItRC  42 

ISN 

0070 

IF! J.GT.30I  GO  TO  132 

ISN 

0080 

00  130  IaJ*DATS 

ISN 

0  081 

130 

RPR  SI  II*  FLOAT!  OPAI  I.6RROPGN! 1-ACI 

ISN 

0082 

132 

ROMTI  U*  ONOSill-  4 PR  31  1 1 

ISN 

0083 

00  135  I*2*0AVS 

ISN 

0084 

136 

RUNTl  II*R0MTI  I-  il  4  ONOStll  -  R  PR  SI  1 1 

ISN 

0086 

TOT  R0*R0MT4  11 

ISN 

0086 

00  140  1*2*  OA  VS 

ISN 

0087 

140 

IF  IRUMTi  II  .GT.TOTRUi  T OTR  b«RUMT<  1 1 

ISN 

0089 

GO  TO  166 

ISN 

0090 

146 

00  150  1*1. Il 

ISN 

0091 

15  0 

TOT Ru*rorRU40N0SI II 

■  c 
c 

c 

RR  PROCESSING 

ISN 

0092 

GO  TO  165 

ISN 

0093 

iSS 

00  160  I*  1* 0 ATS 

ISN 

0094 

T OT Rtf*T OT R040N0S!  Il 

ISN 

0096 

160 

IF  1 INA IN .NE« RRR 1  R  UN  T 1  I  1  -  TO  TRU 

ISN 

0097 

165 

xNUIXTNI*  TOT  R<1 

ISN 

0098 

GO  TO  103 

ISN 

0099 

c 

2  00 

1T0T*ITN- 1 

ISN 

0100 

L 

ifiipass.eu.ii  iron*  t tot 

ISN 

0102 

IF!  IPASS.EU.2  1  IT OT 2* IT  OT 

ISN 

0104 

IFIIPASS.EU.1I  PRi:t5*0.  0 

ISN 

0106 

■RITE! 6. 1931  ITOT 

ISN 

0*07 

193 

FOR  NAT!  //.SX.  'TOTAL  NUNSER  OF  I  TE  NS  UNDER 

ISN 

0108 

RETURN 

ISN 

0109 

E  NJ 

ISN 

0002 

SUBROUTINE  3T PUT!  M  T.  3N0RS  .SuO.PR 1CEI 

C 

c 

this  subroutine  print  the  final  resllts  of  ' 

ISN 

0003 

CON  NON  /BLX  1/  IT  OT  »N  OF  tb  JOG  FT  /0LK2  /XNUIS  00 

♦1  6001  /8LK.3/N0PSfSRF  il  UEEFF  »IRK/8l  kS/LOOPZ 

ISN 

0004 

CONNON  /SLR  0/  ISUN  >  IP  ASS 

ISN 

0006 

OINENSION  RIT1S0OI. SNI 41 

ISN 

0006 

* 

DOUBLE  PRECISION  ENORS.SOO 

ISN 

0007 

W 

l  TE  N-  0 

ISN 

0008 

■R1TEI8.100I  PRICE* EN ONS. SCO 

ISN 

0009 

100 

FOR  NAT! // *SX  *  *PRICi 1 A  I* *• * F I 6.2 *6 X»* El N OR 

4*  •»  FiO.  5*  //  1 

ISN 

0010 

■  RITE!  6  *1101  LOOPC 

ISN 

0011 

110 

F  OR  NAT!  /  *SX#  *NUNBc  R  OF  ITERATIONS  IS*.IA 

c 

«•« 

ISN 

0012 

CALL  TIME 

ISN 

0013 

RERlNO  1 

ISN 

0014 

IF  ! IPASS.E 0. 21  STOP 

ISN 

0016 

1PASS-XPA  SS41 

ISN 

0017 

50 

RETURN 

ISN 

0018 

ENO 

CONSIDERATION  ISM«i//l 


C.2.5 


1 


ISM 

0002 

r 

SUONOUTINE  TABLE 

w 

c 

f 

connotation  of  the  Cjnulative  suns  of  poisson  distributor 

c 

POISON!  It  Kl  »  CUMULATIVE  PROBABILITY  OF  DEMANDS 

OF  OUANTITT  K-l 

c 

F  OH  ITE  M  I 

c 

r 

for  k  greater  than  nf^ohi,  pois)n«i*k»  •  i 

1SN 

0003 

W 

COMMON  P0IS0N(300t  1  00  1 

ISN 

0004 

COMMON  /bLK  1/  ITOTtMUP  .BUDGET /BL  K2  /XNUt  *00  I  t  NAP  1 500  1  tCOST  (  500 1 1 NFPO 

♦1  So  01 /BIX  3/ NOPS*  SKF  tIUEEFFt INK 

ISM 

ooos 

COMMON  /BLKO/ ISUM,  IPA  SS 

ISN 

0006 

COMMON  /PRINT  /MINPOTt  NXMUtPITER 

tSN 

ooor 

DOUBLE  PRECISION  PJ  IS  ON  t  T  E  FN  •  MM  ■»  SUN  §  TE  RL 

ISN 

0000 

TERL-S.0-11 

ISN 

0  309 

do  io  i-i.itot 

ISN 

0310 

R«MOfK.U  XMUI  II  1 

tSN 

OOli 

term-oexpi-mnni 

ISN 

0012 

PUI  SON!  I, 1 l-T  ERM 

ISN 

0013 

OU  20  3«2# 100 

ISN 

0014 

1  FI  TERMtlTt  TERLtANOtP  I  ISON  4  1 1 3”  11  tGT»9«D-01l  GO 

TO  30 

ISN 

0316 

TER M— TERn-nnm/OFLOA  T4 3-11 

ISN 

0017 

POISON!  It  JI-POXSONi  It  J-l  I*  TERM 

ISN 

0010 

20 

CONTINUE 

ISN 

0019 

NFPO!  II -100 

ISN 

0020 

60  TO  19 

ISN 

0021 

30 

NFPO!  II -3 

ISN 

0022 

00  40  JX-JtlOO 

ISN 

0023 

40 

P  01  SO  Nl  It  OKI-  It  D  »  03 

ISN 

0024 

10 

IF!  MXNutEtft  11 

t-RI  TE 16 , «  1  tiNFPOtlli  NAP!  11,  X  NO  III 

ISN 

0026 

RET  URN 

ISN 

0027 

E  NU 

ISN 

90  02 

SUN  ROUTINE  JNKlTIKlTt PRICE  1 

c 

c 

INCREASE  KIT  VALUfeS  TO  ENSURE  THAT  THE  PROBABILITY  OF  ANT  ITEM 

c 

p 

artICH  a  ILL  NAVE  N  Aw  F  JF  THE  FLEET  03aN  IS  LES5  THAN  50* 

ISN 

0003 

V 

DIMENSION  KITI3001  tNlTOlfOOl 

ISN 

0004 

COMMON  Pi)  I S ON  1  3 00 1 1  00  1 

ISN 

0005 

COMMON  / 0LK 1/ ITOTtMUP tttODGET/BLK2 /KMUl 5  00* t NAP 13 091 tCOST! 5001 tNFPO 

t!  3  9  01  /BLK  3/N0PS,SKrt[JEEFF tlKK/  BLK9/VPfiR0t SOCZHO 

ISN 

9006 

COMMON  /0LKO/ XSUNtIPA  $S 

ISN 

0007 

COMMON  /BLK 12 /ITOTl tl T0T2 

ISN 

0096 

DOUBLE  PRECISION  PJIS3N  , VP2ER0, SDO/RO  tVPtSOC 

ISN 

0309 

DAI  A  KIT  0/300*07 

ISN 

0019 

IF!  IP  AS  St  NE#  1 1  GO  T9  i 

ISN 

0  312 

DO  1  I*  It  IT  OT 

ISN 

0013 

1 

K  IT!  1 1-0 

ISN 

0014 

2 

CALL  CALC!  K  IT  Of  VPZE  RO  •  SO  02  FO  1 

ISN 

0013 

■RITE*  6t991ll  KI  TO!  1*1  »  lK-lttTOTI 

ISN 

9016 

991 

FORMAT  l//lSXt  •  E1PTY  K  IT  * /!  3a  t2Sl  2x  tl  2 1  1 1 

ISN 

0017 

•  RITE  1 6t  992 1  V  PI ERUt  SCO/RO 1 1 UEEFF 

ISN 

0010 

992 

FORMAT!  /  15  Xt'VP't  FlO.StSXt  •  SUO«  tF  10.  5  >3  X  t*  IUEEFF  •  1 1 4 1 

ISN 

0019 

NEPS-MJP/  2 

ISN 

0020 

DO  10  ITEM-i,  ITOT 

ISN 

0021 

20 

NUM-NEPS-NAP!  ITEM  It  XI T< ITE  FI tl 

ISN 

0022 

IF!  NUN. GT tNFPO!  ITEM  II  O0  TO  10 

ISN 

0024 

IF!  POISON!  ITEM.  MUNI  t*  T..5I  GO  T  0  10 

ISN 

0026 

KIT!  ITEM!  -KIT!  ITE  Ml  tl 

ISN 

9  327 

PRICE -PRICE  tCOST!  ITEM  1 

ISN 

0020 

GO  TO  20 

ISN 

0029 

10 

CONTINUE 

C.2.6 


1 

I 


IS* 

003  u 

Z  ILL  CALCIKIT.VP.  SJOI 

ISM 

0031 

»R1TE  16*911  1  Kirill*  1*1*1  TOT  • 

ISM 

0032 

91 

FORMAT!/// IS  X* *  INI Tt At  A I T* */t  3X *291 2X  *1 21  1 1 

I5N 

0033 

■  RITE  16*  9941  IP  .SCO  •  t  tEEFF  .PRICE 

ISN 

0  034 

99* 

FORM  ATI  /  1SX*  •  V*  •»  F  10.  3*5*.  •  SO}  •  »F  10*  5  .5  X*  *  IUEIFF*  •  14.5*  » 
»  •PRICE • •  F  12.21 

ISN 

0035 

RET  URN 

ISN 

00  Jo 

3 

IFl IT0T1.EU.IT0T2 1  U  TO  2 

ISN 

0030 

1 1*  iron*  i 

ISN 

0  039 

00  «  N-l  1*  1T0T2 

ISN 

00*0 

* 

X  IT  1 M  »•  0 

ISN 

0041 

60  TO  2 

ISN 

00*2 

E  NO 

ISN 

00  32 

SUBROUTINE  0PT1M1KIT.  VP*SDO*PRICE  1 

L 

c 

r 

MAIN  OPTIMIZATION  ROUTINE 

ISN 

0003 

W 

DIMENSION  KXTISOO I.XT 37 181 3001 

ISN 

000* 

DIMENSION  B1NTER1  10> 

ISN 

0005 

COMMON  /atlCO/ ISUM.1PA  SS 

ISN 

0006 

COMMON  /dLK  l/IT0T.*UP.dUDG£T/BLK2/AMU(  S  OO  1  *  NA  PI  S  00 1  *  COST  1  5001 »  NFPO 
*1 SOOI/at*3/NOPS*S«F •! JfcEFF *I*x/Bk  *4/lFt *6.1  NT  * RA TI J/BLKS/LOOPC 

ISN 

ooor 

COMMON  /BLK6/B1NTE*  /  8 1*9/9  Pi  ERO.  SO  Qt  «0 

IS* 

0000 

COMMON  /BLK8/VPT  *  *  P6  0  At  *  SO  06  Ot 

ISN 

0009 

COMMON  /COUNT /XCOUvT.  lADDt  2001 

ISN 

0O10 

COMMON  /PRINT /MXNPoT* MaMU.FITER 

ISN 

0011 

r 

0  00  He  £  PRECISION  VP.vPbOAL.V  FOELT  .VPTEST.  SOO*  S0060L. SOOOEL. SOOTES 
♦  .VPZEAO.SDOZRO.V  pt. 

ISN 

0012 

V 

990 

FORMAT!///  ISX.'KIT  3718*  /  /  1  3A  *251  2X  *1  21  II 

ISN 

0013 

99* 

F  OHMATl  T3**  ITERATION*  .14. 5«*  *VP  -**F10.S.  3*  •  *S00  »«.FlO.S.  3x* 

*  *  PRICE  ■«»  F12*2 *3 X** TOTAL  I T£ PS  AOUEO  -*.  1 3. /I  17.  191  15  *  II 

ISN 

001* 

c 

999 

FORMATl  5*. ‘TARGETS  T3  HE  A TT AI NEO * »/ / *5 X. •  C 0 ST *F 10. 2 *5 X* • NORS 
.GOAL*  *.FiO.  5.  SX»*$D06  JAk*  •  *F  10.  S*  //I 

ISN 

0013 

V 

IFk  A6-0 

ISN 

oo  to 

tOO  PC-3 

ISN 

0017 

c 

BUDGET- 0.0 

c 

c 

*  IT  5718  CALCULATION 

ISN 

oota 

OO  10  ITEM- l. ITOT 

ISN 

0019 

KT37181  ITEMI  —  IFlKlA  Mul  ITEMI.SKF*.  -1 

ISN 

90?0 

IF*  *T S7 18IITEMi.tT.  11  KT57  181  ITEM  1-1 

ISN 

0022 

8U0GET-8U06ET  *COSTl  IT £N  1  -F  LO  ATI  * T  *7161 1  TEM 1  1 

ISN 

0023 

r 

10 

CONTINUE 

ISN 

002* 

w 

CALL  CALC1K7S716.VPG01L.S0CG0LI 

ISN 

0025 

■  RITE 16.9901  1KTS718I  I X 1 . I X- 1. 1  TO T 1 

ISN 

0026 

■  HI TE46 .999 1  BUDGET .« PGOAL  .S0060L 

ISN 

0027 

IF  1  1PASS.E0.  21  VPfi-l»l«VPGQAL 

ISN 

0029 

c 

IFIIPASS.CO.il  BUOtET -BUDGET-. 9 

ISN 

0031 

V 

vPO£lT-l  X PZERO- VPGJ Ae  »/5.0 C 

ISN 

0032 

X  PT  EST- VPZE RO-VPOEl  T 

ISN 

0033 

I  VP-0 

ISN 

003* 

SOJOEL-1  SOOZRO-SOOGOl  1/5.0  C 

ISN 

0035 

sootes-soozro-soooe  l 

ISN 

0036 

V P- VP ZERO 

ISN 

0037 

SOO -303  ZAO 

ISN 

0  030 

ISO  0-1 

ISN 

0339 

I OE  LT  A-  1 

ISN 

00*0 

60  TO  223 

ISN 

00*1 

210 

IFl  VP.LE.  VPTESTI  GO  TO  220 

ISN 

00*3 

1  FI  SDO.fLE .  SOOTE  SI  .0  TO  222 

ISN 

00*5 

LOOPC-LOOPC*! 

c 


C.2.7 


ISN  0  046 
1SN  0047 

ISN  0049 
ISN  0050 
ISN  005i 

ISN  00  52 


219  CN.I  AO  04  KIT.  VP.SOJ  »P  RICEI 
I  FI  MlTER.EU.il 

4  nRI TE( 6 .994  I  LJOPC »9P» SOO.PRI Ct.ICOUNT. « 1 ACOI II . 1*1* XCOUNTI 
DO  216  !■  1.  ICOUNT 
2  16  1  ROM  II  >11 
ICOUNTaO 

•  a. 

I  FI  MlTE  R.EO.il 
♦  CALL  TlMEl 


ISN  005* 
ISN  00 56 
ISN  0057 
ISN  0056 
ISN  0059 
ISN  0061 


iniFLA6.EU.il  GO  TO  ISO 
60  TO  210 

220  4PTEST**PT£$T-9P&El T 

IVPalVP+i 

1FI  S00.6T.S00TESI  .3  TO  223 
222  SOOTES*SOOTES-SOOD£L 


SELECT  PROPER  •RATIO*  A  90  *1  NT  * 


ISN  09£2 
ISN  0065 
ISN  0064 
ISN  0066 


ISN  0067 


ISN  0060 


I  SO  OalSOO  +  l 
223  ION  AN«6aI  VP4IS00 

IF!  I0NAN.6T.36I  IBRAn*36 

60  TO  I  334.  305.  306.  30>.  3  06  .307.  303.504. 305. 306  .3  06.307. 
13 03. 393. 3 04 .3 09.3 06 • 5 37 .302* 303. 3 C3. 3 04 .3 06 .3 07. 

2  3  02.302.3  02.3  02.304  .3  07*30  1.301.3  Cl.  301  .3  01. 3  3  51  .IsRAN 
2 So  60  TO  t 311. 311.311. 312.312 *3 12.31  1.3 11. 312. 312.3  13.3 13* 
13  lx.3  12.3  12.3  13*314.3  14.3  12.  512.3  13.  313.3  14.3  15* 

23  12.313.314.3  14*316 .3  16. 312.  313.3  14.  3  15  .3  16 .3  331  .  ISRAN 
30 1  4PTEST*  0.  04  00 


•RATIO*  ANU  •  tNT*  SET  TABLE 


302  RAT 
60 

303  RAT 
60 

304  RAT 
60 

305  RAT 
60 

3  06  RAT 
60 

307  NAT 
SOO 
60 


IJaio. 

TO  250 
10*2  5. 

TO  2  SO 
10*5  0. 

TO  250 
10* 150. 

TO  2  5 0 
10*5  00. 

TO  250 
10*1000. 
TES-0.0400 
TO  250 


311  IMT-ITOT/44  1 
60  10  2 10 

312  1NT*IT0T/841 
60  TO  2  10 

313  I  NT  *ITOT/ 1641 
60  TO  210 

314  I  NT *1 TOT/  3241 
60  TO  2  10 

315  I  NT*  I  TOT/ 64  41 
6  0  TO  2  10 

316  INT*l 

60  TO  210 
335  RAT  10*50. 

I  FI  4  SU06E  T*  PR  ICE!  .t*C.  40000.1  60  TC  534 
INT*  1 
60  TO  215 

334  IF!  (8U06ET-FRICEI.6E. 20030C. I  60  TO  335 
I  NT *1  TOT/6441 

6  0  TO  215 

335  I  FI  4  UUOGE  T*  PR  IC  E 1  .6  E»  SOOOOCO.I  60  TO  33  6 
INT  -1TQT/  1641 

60  TO  215 

336  INT*IT0T/841 
60  TO  2  13 

350  1FLAG*0 
RETURN 
e  NO 


ISN  0002 


SUBROUTINE  ADO! KIT* VP .SOO.PRICE I 
C 

C  THIS  SUBROUTINE  AOOS  INT  DIFFERENT  ITEMS  TO  THE  KIT  BV 

C  MARGINAL  ANALYSIS 

C 


ZSN 

0003 

DIMENSION  TSl 5001. K ITl 9001 

ISN 

0004 

OIHENSION  BINT-Rt  lOltS  11  S0t*SSltS0>  .S  10  El  C  SOI 

ISN 

OOOS 

COMMON  POISON!  300.1  10  1 

ISM 

000b 

COMMON  /ULK  0/  IS  OM  .  I  PA  S3 

ISN 

0007 

COMMON  /0LK  1/ITOT.MUP  .6U0GET /BLK2VKMU1  5  001.  NAP  IS  001.  COST!  3001  .NFPO 

*1  50  01 /BLK  3/ HOPS  *  SKF  *1 oEEFF .1 KK/ BL  K*/ IFL  AG  >  I  NT  •  KA  T 10 

ISN 

oooa 

COMMON  /BLK 6/81 NT£R 

ISN 

ooov 

C3MM0V  /BLK7/S1  /6LK6/VPT.  VPG3AL.  S3  3GUL 

ISN 

0010 

COMMON  /COUNT /IC0U1T.  IADO!  2001 

ISN 

0011 

DOUBLE  PRECISION  PJ  IS  JN.TS  .OIF.  VP  .A.SDO.SDOP.T. 

*  Sl.SSl  .SlTEMP.SS  IT  SM.  SlSUM.  SS1SUM.S10EL 

ISN 

0012 

DOUBLE  PRECISION  VP31AL.SD0G3L  .EPS  •  V  PT 

ISN 

0013 

DATA  INTER/l/.MFLAG/  3/.EPS/  1.0-12/ 

ISN 

0014 

V 

NOPSl-NOPS*  1 

ISN 

001S 

I  FI  IPASS.  EG.  2  1  BTE*  P»  dXNTE  *1  INTER  1 

ISN 

00  17 

IF  4  I PASS. EU.  11  BTi'MP  *bUDG  ET 

ISN 

0019 

IF!  IPASS.  Ett.  l.ANO.d  INTER!  t  UTERI  .LT.  BUDGET  1  BT  E  MP -b  INTER!  INTER  1 

ISN 

002l 

r 

S  IT  ENP*  0.  DO 

ISN 

0022 

V* 

00  3  0  US1«  1*  NOPS 

ISM 

0023 

30 

S  lTEMP*SlTEMP*Sl!  J>  II 

ISN 

0024 

StSUM-N0PS4SllN0PSl  1-SlTEMF 

ISN 

002S 

C 

00  99  ITm*I.ITOT 

ISN 

Oo2t> 

1XL -KIT! ITM ). 1 

ISN 

0  027 

lXH-KlT! I  TNI*  ! NOPS  —1 1 aNAPI ITNI *2 

ISN 

0020 

IF!  IXN  .UE.  NFPO! IT  Mil  1 XMaNFPO ! I TM 1 

ISN 

0030 

IF!  IXL  «GE.  nFPOIITMII  t XL *NFPO 1 1 TM 1 

ISN 

0032 

DIF*  POISON!  I  TM.  1X4  l-:POI$OA!  ITM.IxLI 

ISN 

0033 

i» 

DO  20  K»l.N0PSl 

ISN 

003* 

K  J»KlTl  ITNI*I  K  -1I.NAPI1TPI  *1 

ISN 

003S 

KJ1  .<J»l 

ISN 

Oo3b 

IF  1KJ  ,gE.  NFPO!  IT  All  «  60  TO  33  3 

I  SN 

0038 

m  POISON! I  TM.KU1.E0. fl.D 01  60  TO  222 

ISN 

0040 

S Si ! K  t«Sl!K  laPOlSOM  ITM.KUII/POISONIITF.kUI 

ISN 

00*1 

GO  TO  2  0 

ISN 

00  42 

222 

SSl!  *1*  1.  00 

ISN 

00*3 

oo  223  Ki-i.iror 

ISN 

00*4 

K  1 1  *K  IT!  K  1 1  +4  K-ll*«APIKl»*l 

ISN 

00  45 

1  FI  KI  •£  0.  IT  Ml  K  X 1  "K  1 1  *  1 

ISN 

00*7 

IF  IKI1.GE. NFPOiKlI 1  GO  TO  223 

ISN 

00*9 

SSllK|a$S 11 Kl aPOXSJNI  KI  .Kill 

ISN 

1050 

IF!  SSl!  Kl  .UE.CPSI  j  <3  TO  32  1 

ISN 

0  052 

223 

CONTINUE 

ISN 

00  S3 

GO  TO  2  0 

ISN 

0054 

321 

S  Si!  K  1*  0.00 

ISN 

00  S3 

GO  TO  20 

ISN 

0056 

333 

$Sl!K  l*$  11  Kl 

ISN 

0057 

c 

20 

CONTINUE 

ISN 

0056 

V 

SS1TEM-0.D0 

ISN 

0059 

DO  AO  ASS  1.  NOPS 

ISN 

0060 

40 

SSlTEMaSSlTEH*SSl!J  SS  ll 

ISN 

006  1 

SSl SOMaNOPSaS $11  NOP  Si 1-SS ITEM 

ISN 

0062 

A»S  SlSUM 

ISN 

00  63 

TSl  ITMlai  |  VP- AlaDMLSl  RAT  10  I*  DIF  1/  OB  LEI  C  OSTl  IT  M  II 

ISN 

006* 

99 

continue 

C 
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’I 


C  A  00  2  OF  THE  IMT*  ITEMS  ml  Th  TmE  HIGHEST  I M  PM  0  vE  ME  H  T  RATIOS 

C  A 00  1  OF  THE  NeXT  HIGHEST  JNT-lNT/4  ITEMS 

C 


ISN 

0065 

INT4*lNT/4 

ISN 

0066 

00  102  Uli  INT 

ISN 

0067 

T  *T  SI  It 

ISN 

0060 

ITN»i 

ISN 

0069 

oo  tot  o«2*iror 

ISN 

0070 

IF(  TSIOl.LT.T »  GO  TO  101 

ISN 

0072 

T  »T  SI  01 

ISN 

0073 

ITN  mj 

ISN 

0074 

101  Z ON  TINJE 

ISN 

0075 

XIT4  ITMI-MITI  ITMI.* 

ISN 

0076 

1  COUNT* 1C0UNT  41 

ISN 

0077 

I  AO 01 ICOU NT  1*  ITN 

ISN 

0078 

IFI  MFLAGaEU.il  C  ALL  MlC  1  KITaVPaSDOt 

ISN 

0080 

PRICE-PRIC£«COSTI  IT  Ml 

ISN 

0081 

T  SI  ITMI  « J.D.D  9 

ISN 

0082 

IFI  PRICE.  GT.BTEMP.  JR.  i  NFLAGaEo.l  aANDa  1  VP .L E . VPG IAL. AND. 

♦  SOO.  L£ a  SOOGOL III  GO  TC  10S 

ISN 

0084 

IFI  INT4.E0.  01  GO  T]  1  02 

ISN 

Q08o 

x  IT  1 ITNI-XITI  ITMI 41 

ISN 

0087 

lCOUMT*  ICOUNT  41 

ISN 

00  88 

1  AO 01  ICOONT  l«ITM 

ISN 

0089 

1FIWFL1G.E0.1I  :au  :*l:»  KIT.VP.S0D1 

ISN 

OOVl 

PR1CC-PRICE4C 0ST1  ITMI 

ISN 

0092 

IFI PN ICE. GT  a  BTE MPa J  R* 1  MFLXG.Eu.l  aANOa  1  VP.lE . VPGJAl. AND. 

♦  SOO.LE.  SOOGOL  1  ll  GO  TG  105 

ISN 

0094 

INT4W  IN T4 •  1 

ISN 

0095 

c 

102  Z  ON  TINGE 

ISN 

00  9o 

CALL  CAlCIKIT.vP* SOO* 

ISN 

0097 

IF  IVP.LE.VPTI  MFL4  G* 1 

ISN 

0099 

c 

RETURN 

ISN 

OlOO 

v» 

105  IFI  II  IPASS.EU.il.  A49.dT£MP.£6.UUUG£TI  .OR. 

♦  INFLAG.  Eu.l  .AND.IVP.Lfc.VPGOAL.AND. 

♦  SD 0*1.6.  SOOGOL  1  11  IFLAG*1 

ISN 

Oi<>2 

call  CALClXITaVPa SO 9  1 

ISN 

0103 

IF  IVP.LE.VPTI  MFLlG«t 

ISN 

0 1  Os 

IFlIFLAG.NE.il  MRITEl6a995l  8TE  HP 

ISN 

0107 

IFIIFLAG.EU.il  «RIT  St  6.994!  PRICE 

ISN 

0109 

PRITEl  6  .9911  KiniXIa  Ik«1.1T0T1 

ISN 

0110 

■  RITE  16.9931  V  P.  SOOa  IGEEFF  .P  RICE  a  I  TN  ,C  OSTI I  TM1 

ISN 

Olll 

s  loen  ii»sii  1 1 

ISN 

0112 

DO  95  U«2  aNOPS 

ISN 

0113 

95  SlOELIJIa  S1IJI-SUJ-II 

ISN 

0114 

f 

■  RITE  16.9961  1  SIDE  Ll  1 1  a  I>  la  NOP  Si 

ISN 

OllS 

V* 

IFlIFLAG.NE.il  INTER.  INTER  »1 

ISN 

0117 

RETURN 

ISN 

0118 

V 

991  FORMAT!  /  1 5  X  a  •  K  IT  •  a  /  •  1  3  X  *2  512  X*  12  1 1  1 

ISN 

0119 

993  FORMAT!  /15Xa»V»*.  FlOaSaSXa  •  S03  •  »F  10*4  aSR  a  •  IUESFF  ••  I  Aa  SX  a 

»  •PRICE* aFl2a 2.5X I* LAST  ITEM  AOOEO*  al 5a SX a •! TEM  COST* • .F 10a 2 a 

♦  //I 

ISN 

0120 

995  FORMAT!  ///  2x  a  •BUDGET  »«*aF  12.21 

ISN 

0121 

996  FORMAT!  5!  4X.6  20.121  1 

ISN 

0122 

994  F  OR  MATl  //  /  2  Xa  ‘PRICE  IF  3PTIMAL  X  IT  *  *  a  F  12.21 

ISN 

0123 

END 
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l 

0029  SYSTEM  C0NTR0L/INP0T  CARD  FORMAT 


I. 


“A"  CARD  :  PROGRAM  CONTROL  CARD 

Description 


Field 


Column 


Type 


flats-"*- 

II. 


III. 


*  TVPe 
I 
k 


-- 

1 

1-5 

15 

Initial  Program  in  Flying 
Hours. 

2 

7-8 

12 

WRSK  Support  Period  in  Days . 

3 

10  -14 

15 

Total  Program  in  Flying  Hours 

4 

16  -17 

12 

Unit  of  Equipment. 

5 

19  -20 

A2 

Run  Sequence  Number. 

6 

22  -23 

12 

Number  of  Kit  Author isations* 

25  -26 

12 

Number  of  Set  Up  Days . 

54 

Al 

A  (Card  Indicator) . 

B"  CARD 

;  KIT 

SUMMARY  CONTROL  CARD 

1 

1  -10 

,,A10 

Kit  Serial  Number. 

54 

Al 

B  (Card  Identifier) . 

C"  CARD 

s  ITEM 

DATA  CARD 

1 

1  -15 

A15 

National  Stock  Number  (Left 
Justified) . 

2 

17 

11 

Note  Code. 

3 

19  -21 

A3 

RR/RRR  Maintenance  Concept 
(Left  Justified) . 

4 

22  -23 

12 

Repair  Cycle  Time  in  Days. 

5 

25  -26 

12 

Quantity  Per  Assembly. 

6 

28  -30 

A3 

Application  Code  . 

7 

32  -36 

F5.4 

Total  OIM  Demand  Rate. 

8 

38  -42 

F5.4 

Deoot  Demand  Rate. 

9 

44 

Al 

Routing  Identifier. 

10 

46  -53 

F8.2 

Unit  Cost . 

11 

54  -56 

A3 

LRU/SRU/EOQ/NOP  Type  Item. 

12 

58  -62 

A5 

Work  Unit  Code  (Left  Justi¬ 
fied)  . 

13 

73  -80 

A8 

Noun .  . 

r  1 


Explanation 

Numeric  input  only#  leading  zeros  as  required,  right 
justified . 

Alphabetic  or  numeric  input  • 

Numeric  Input,  decimal  assigned  by  system,  leading 
zeros  as  required,  right  justified. 
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Phase  1  Data 


1 

1 

25927.00 

1.65423 

RR 

LRU 

2 

1 

3054.00 

1.01932 

RR 

LRU 

■» 

1 

600.00 

1  .03737 

RR 

LRU 

4 

1 

6166.00 

1  •  048  07 

RP  • 

LPU 

5 

1 

11107.00 

1  •  C2b  75 

RR 

LRU 

6 

1 

3162.00 

0.03737 

RR 

LRU 

7 

1 

13900.00 

5  .01276 

RR 

LRU 

8 

1 

2332.00 

p  .04025 

RR 

LRU 

9 

1 

2868.00 

5  .03599 

RR 

LPU 

1C 

1 

3550.00 

5.03599 

RR 

LRU 

11 

1 

8764. CO 

1  .07348 

RR 

LRU 

12 

1 

12366.00 

1 .04025 

RR 

LPU 

13 

1 

1627.00 

1.25553 

RR 

LRU 

19 

1 

2369.00 

1.16422 

RR 

LRU 

1? 

1 

5258.00 

j  .11396 

RR 

LRU 

16 

1 

1325.00 

0.03231 

RR 

LRU 

17 

1 

2975.00 

0.12086 

RR 

LRU 

18 

1 

2520.00 

1.10879 

RR 

LRU 

19 

1 

303.24 

1  .13075 

RR 

LPU 

20 

1 

799.00 

3  .05048 

RR 

LRU 

21 

i 

1424.00 

3.06348  ' 

RP 

LRU 

22 

i 

1056.00 

1.05773 

RR 

LRU 

23 

i 

845. 30 

1.04025 

RR 

LRU 

24 

i 

919.80 

D. 04726 

RR 

LRU 

25 

i 

717.70 

0  .03749 

RR 

LRU 

26 

i 

927.00 

1.09522 

RR 

LRU 

27 

3126.00 

1.10223 

RR 

LRU 

28 

i 

2122.00 

0.06049 

RR 

LRU 

2° 

l 

6891.00  0.21332 

RR 

LRU 

30 

i 

582.00  0.11569 

RR 

LRU 

31 

2 

374  5  .90 

1  .10948 

RR 

L»U 

32 

1 

692.70  0  .05195 

RR 

LRU 

33 

1 

1333. GO  0.07383 

RR 

LRU 

34 

1 

1333.00 

0.08889 

RR 

LRU 

35 

1 

3202. 0u 

3  .22650 

RR 

LRU 

36 

1 

516.00  0  .71875 

RR 

LRU 

37 

1 

1241.00 

3  .29003 

RR 

LRU 

38 

1 

2105.00 

0 .05750 

RR 

LRU 

39 

1 

824.50  0  .07820 

RR 

LRU 

4  P 

1 

1342.00 

0.06014 

PR 

LRU 

41 

2 

3914.00 

1.17480 

RP 

LRU 

42 

2 

891.73 

1  .17227 

RR 

LRU 

43 

2 

933.80 

1  .04692 

RR 

LRU 

44 

2 

2058.00 

0  .19343 

RR 

LRU 

45 

2 

349.1  0 

5  .04653 

RR 

LRU 

46 

2 

1242.00 

3 .06923 

RR 

LRU 

47 

2C 

471.70 

o.38ieo 

RR 

LRU 

*8 

5 

304.30 

1  .27370 

RR 

LRU 

49 

1 

1317.00 

3  .04393 

RR 

LRU 

50 

2 

2083.00 

3 .05865 

RR 

LRU 

51 

2 

104*. CO 

0  .15502 

RR 

LRU 

52 

u 

637.90 

0.13156 

PR 

LRU 

53 

2 

487.90 

3.38318 

RR 

LPU 

54 

2 

9901.00 

D. 32798 

RR 

LRU 

55 

2 

275.30 

D  .06624 

RR 

LRU 

56 

1 

2975.00 

3  .17043 

RR 

LRU 

57 

4 

1130.00 

3.96370 

RR 

LPU 

58 

1 

886.00 

0  *07b24 

RR 

LRU 

59 

5 

320.30 

3.14U87 

RR 

LRU 

6C 

1 

1303.00 

0.04174 

RR 

LRU 

81 

1 

1300.00 

0 .05566 

PR 

LRU 

62 

1 

456.90 

3 .06175 

RR 

LRU 

63 

1 

309.70 

0 .?4335 

RP 

LRU 

64 

1 

945.50 

D .20493 

RR 

LRU 

65 

1 

286.00 

9.12121 

RP 

LRU 

69 

2 

412.  0C 

RR 

LRU 

70 

2 

262.10 

0*06992 

RR 

LRU 

71 

2 

1358.0  C 

0 .84111 

RR 

LRU 

72 

1 

1874.00 

Oil  728  4 

RR 

LRU 

73 

1 

28G.00 

0  .164  2  2 

HR 

LRU 

74 

1 

3658.00 

0  *06382 

RR 

LRU 

7*> 

1 

335. CC 

0  .09579 

RR 

LBU 

76 

2 

2129.00 

0.37237 

RP 

LRU 

77 

1 

2166.00 

0  .14973 

RR 

LRU 

78 

1 

1158.00 

0,12788 

RR 

LRU 

7o 

1 

3550.00 

0.19251 

RR 

LRU 

80 

1 

3290.00 

0443493 

RR 

LRU 

81 

1 

80  2.4  0 

0.0  4715 

RR 

LRU 

62 

1 

in  5.00 

0.25392 

RR 

IRU 

83 

2 

267.80 

0.23000 

RR 

LRU 

64 

1 

1051.00 

0  .54383 

RR 

LRU 

05 

1 

875.50 

0.13524 

RR 

LR'J 

86 

2 

7B6.50 

0 .24978 

RR 

LRU 

87 

1 

1401.00 

0  .14846 

RR 

LRU 

€8 

1 

1585.00 

0.06290 

RR 

LRU 

09 

1 

6180.00 

0  .64628 

RR 

LR'J 

90 

1 

2553.00 

0  .05186 

RR 

LRU 

91 

2 

2292.00 

0.11707 

RR 

LRU 

99 

1 

831.2U 

0  .05681 

RR 

LRU 

93 

1 

641.69 

0  .11638 

RP 
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